EUV and X-ray Spectroscopy of 
the Active Sun 




Claire L. Raftery 
[S chool of Physicsl 
University of Dublin, Trinity College 



A thesis submitted for the degree of 
PhilosophieE Doctor (PhD) 
2009 December 



Declaration 



I, Claire L. Raftery, hereby certify that I am the sole author of this 
thesis and that all the work presented in it, unless otherwise refer- 
enced, is entirely my own. I also declare that this work has not been 
submitted, in whole or in part, to any other university or college for 
any degree or other qualification. 

This thesis work was conducted from October 2006 to December 2009 
under the supervision of Dr. Peter T. Gallagher at Trinity College, 
University of Dublin. 

In submitting this thesis to the University of Dublin I agree that the 
University Library may lend or copy the thesis upon request. 



Name: Claire L. Raftery 
Signature: 



Date: 



Summary 



This thesis strives to improve our understanding of solar activity, 
specifically the behaviour of solar flares and coronal mass ejections. 
An investigation into the hydrodynamic evolution of a confined solar 
flare was carried out using RHESSI, CDS, GOES and TRACE. Evi- 
dence for pre-flare heating, explosive and gentle chromospheric evapo- 
ration and loop draining were observed in the data. The observations 
were compared to a 0-D hydrodynamic model, EBTEL, to aid inter- 
pretation. This led to the conclusion that the flare was not heated 
purely by non-thermal beam heating as previously believed, but also 
required direct heating of the plasma. An observational investigation 
in to the initiation mechanism of a coronal mass ejection and eruptive 
flare was then carried out, again utilising observations from a wide 
range of spacecraft: MESSENGER/SAX, RHESSI, EUVI, Corl and 
Cor2. Observations provided evidence of CME triggering by internal 
tether-cutting and not by breakout reconnection. A comparison of 
the conflned and eruptive flares suggests that while they have differ- 
ent characteristics, timescales and topologies, these two phenomena 
are the result of the same fundamental processes. Finally, an investi- 
gation into the sensitivity of EUV imaging telescopes was carried out. 
This study established a new technique for calculating the sensitivity 
of EUV imagers to plasmas of different temperatures for four different 
types of plasma: coronal hole, quiet sun, active region and solar flare. 
This was carried out for six instruments: Proba-2/SWAP, TRACE, 
SOHO/EIT, STEREO A/EUVl, STEREO B/EUVl and SDO/AIA. 
The importance of considering the multi-thermal nature of these in- 
struments was then put into the context of investigating explosive 
solar activity. 
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solid line), RHESSI (green triangles), and the two components of 
the SAX fit: high temperature (blue diamonds) and low tempera- 
ture (red asterisk), (e) Flare emission measure from GOES (black 
solid line), RHESSI (green triangles), and the two components of 
the SAX fit: high temperature (blue diamonds) and low tempera- 
ture (red asterisk) 11691 



6.6 Observations published in Sui & Holman (2003) that display sim 



ilar behaviour to the observations presented in this chapter. Top 
panel: RHESSI light curves in three energy bands (from top to 
bottom): 3-12, 12-25, and 25-50 keV, scaled by 2.0, 0.5 and 1.0 
respectively. Middle panel: Time histories of the loop height (ob- 
tained from 10-12 keV images) and the coronal source height (ob- 
tained from 10-25 keV images). Bottom panel: Height of the loop 
and the coronal source at different energies at 23:11:00 UT. The 
horizontal bars represent the energy bandwidths of the RHESSI 
images 11711 
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6.7 RHESSI images in different energy bands from |Sui fc Holman 



(2003). The three contours on the solar disk indicate the top of 



loops in the energy bands (contour line shade from light to dark): 
6-8, 10-12, and 16-20 keV. The contours above the limb are (from 
light to dark) 10-12, 12-14, and 14-16 keV bands. The asterisk 
marks the centroid of each source. The crosses mark the two foot- 
points of the X-ray loop 

6.8 Schematic diagram showing the asymmetric eruption of a CME 



that fits well with the observations presented in this chapter. (Tri- 
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pathi et al, 2006) IT75] 



7.1 Schematic diagram from Warmuth (2007) describing the various 



coronal disturbances that are known to occur during a solar erup- 



tion (adapted from Aurass et al. . 2002 ) 
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Chapter 1 
Introduction 



This Chapter introduces the fundamental physics and concepts that are discussed 
in this thesis. This begins with a general introduction to the Sun and the var- 
ious layers of the solar atmosphere and is followed by a detailed discussion of 
the solar corona and solar flares. The physical conditions required for flares are 
presented, including an introduction to concepts such as magnetohydrodynamics, 
magnetic reconnection and chromospheric evaporation. Finally, a discussion of 
coronal mass ejections and their connection to solar flares is presented. 
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1.1 The dynamic Sun 

The Sun has piqued the interest of humankind for thousands of years. Records 
from the Stone Age show that even then, the behaviour of our nearest star was 
noted and utihsed. Newgrange, for example, is a megahthic passage tomb located 
in Co. Meath, Ireland. Built in ~3,200 BC it is beheved to be 200 years older 
than Stonehenge, making it the oldest surviving "solar observatory" in the world. 
While not an observatory in the modern sense, the construction of this place of 
worship required a detailed understanding of the behaviour of the Sun and its 
motion relative to the Earth. At dawn on the days surrounding the winter solstice, 
the tomb is illuminated by light from the Sun. This happens as a result of a near 
perfectly aligned roof box. 

Much has changed since Neolithic times. Today solar observatories utilise 
cutting edge technology to make high quality observations of the Sun. These 
high resolution, high cadence data have revolutionised not only solar physics, 
but stellar astronomy, geophysics and planetary sciences to name but a few. 
With recent advances in space technology, the dynamic nature of the Sun, and 
its atmosphere in particular, is only now beginning to come to light. Explosive 
releases of energy, charged particles and the solar wind all have an impact on the 
Earth. This, of course, is in stark contrast to the image of the Sun most people are 
familiar with. While it is widely understood that the Sun is the source of energy 
and heat for the Earth, few people understand the detrimental impact "space 
weather" can have on Earth. From loss of communications satellites to electricity 
grid failures, even to the interruption of long-haul flights, the full understanding 
and accurate prediction of solar storms is essential to the future well being of this 
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planet's inhabitants. 

Although many details of solar behaviour remain elusive, the general under- 
standing of the Sun and its atmosphere has improved significantly in the past 
number of years. Revolutionary spacecraft both old and new such as the Solar 



Maximum Mission (SMM; Bohlin et al, 1980), the Solar and Heliospheric Ob- 



servatory (SOHO; Domingo et al, 1995) and the Reuven Ramaty High Energy 



Solar Spectroscopic Imager (RHESSI; Lin et al. 2002) are the cornerstones of a 



vast fleet of solar observatories. Between ground- and space-based observatories, 
there are instruments designed to study almost every aspect of solar activity. 
The solar corona is the focus of much investigation and instruments such as the 



Transition Region and Coronal Explorer (TRACE; Handy et al. 1999 ) and the 



Coronal Diagnostic Spectrometer (CDS; Harrison et al, 1995) on board SOHO 



are dedicated to the study of the upper layers of the solar atmosphere. 



The Sun's atmosphere, the layers of which are shown in Figure LI, is defined 
to be the part of the Sun that lies above the visible surface, or photosphere. 
It can be divided into four regions based on their differing physical properties. 
Thermodynamic properties such as temperature and density are highly sensitive 



to height above the photosphere, as shown in Figure 1.2 This can have a signif- 
icant impact on the composition and characteristics of the plasma in each of the 
layers. The density and temperature gradients also affect a parameter known as 
the plasma /3. It is given by the ratio of the gas to magnetic pressure: 



/3 



gas 



B 



E2/87r 



(1.1) 



where Pqas is gas pressure, Pb is magnetic pressure, B is magnetic field strength. 



3 



1.1 The dynamic Sun 




Figure 1.1: Cartoon showing the layers of the Sun from the core to the outer 
atmosphere. 

He is electron density, fcs is Boltzmann's constant and T is temperature. The 
density change though the atmosphere (An ^ 10^^ cm~'^) is greater than that of 
temperature (AT 10^ K) or the square of the magnetic field strength (Ai?^ 
10^ G). Therefore between the photosphere and the corona the (3 value drops 
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Figure 1.2: Temperature (Te) and density (Ne) of the solar atmosphere as a 



function of height ( Gabriel Sz Mason 



1982 



Gallagher 1999). 



1.1.1 The photosphere 

The photosphere is a cool, dense region and is the optically thick layer of the 
atmosphere that is seen when viewing the Sun in optical continuum. The visible 
spectrum of the Sun is in excellent agreement with a black body radiator with an 



effective temperature of ~5,800 K (Figure 1.3) and has a density of ~10^^ cm~^. 
The photosphere is opaque and emits a continuous spectrum which is crossed by 
Fraunhofer lines. Granulation is a feature of the photosphere. This is the term 
used to describe the photospheric manifestation of the large convective motion 
that occurs in the convection zone beneath the photosphere. The continuous 
shifting of plasma by convective motions acts to tangle and stress the magnetic 
field, increasing its non-potential energy. Another interesting feature found in 
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Observed solar tlux per unit wavelength 




3 4 
log(wavelength X [A]) 

Figure 1.3: The observed solar flux through the various layers of the atmosphere. 
The photospheric flux compares very well to a black-body spectrum at a temper- 



ature of 5,762 K (Aschwanden, 2004). 



the photosphere are sunspots. These are large concentrations of magnetic flux 
which appear as dark regions in white light and continuum images. Their dark 
appearance is as a result of the suppression of convection in those areas, result- 
ing in temperatures of 3,000-4,000 K which is much cooler than the surrounding 
photosphere. Sunspots can often be divided into two parts - the central, dark 
umbra, where the magnetic fleld is approximately normal to the photosphere and 
the surrounding, lighter penumbra, where the magnetic flelds are more inclined. 
Figure 1.4| shows a composite image of two interacting sunspots taken with the 
Solar Optical Telescope (SOT) on board Hinode. The top panel shows the chro- 
mospheric Call line. In this image you can see the dark umbra of the northern 
sunspot and chromospheric granulation patterns. The flare can be seen as the 
bright elongated feature stretching laterally between the sunspots. The middle 
frame shows G-band (optical) emission. Both sunspots are clearly visible, along 
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Figure 1.4: Three different views of two interacting sunspots that resulted in a 
solar flare, observed with Hinode/SOT. Top panel shows the chromospheric Call 
line, the middle frame shows G-band (optical) emission and the bottom panel shows 
the magnetogram image of both sunspots with their opposite polarities (black vs. 
white) . 
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Figure 1.5: Ha image of the Sun showing a prominence (top left on the limb) 
and filaments (dark regions on disk). Credit: Jack Newton. 



with granulation patterns. Note there is no indication of a flare occurring in this 
image. The bottom panel shows the magnetogram image of both sunspots of op- 
posite polarities (black vs. white). The vertical nature of the fleld in the umbra 
is clear here as there are no magnetic fleld measurements made in this region. 
The penumbra are showing strong magnetic fleld strengths. Large magnetic fleld 
gradients exist in the mixed polarity region in the center of this image making it 
the location of the solar flare observed in the top panel. 



1.1.2 The chromosphere 



It is clear from Figure 1.4 that the outline of sunspots are not as well deflned 



as height is increased from the photosphere into the chromosphere. At densities 
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Figure 1.6: Image of the structure of the solar chromosphere taken with Hin- 
ode/SOT. 

of ~10^^ cm"'^ and temperatures of ~10'^ K, the chromosphere is dominated by 
absorption lines and continuum emission. The strong Ha line shows bright regions 
called plage in the vicinity of active regions and sunspots. Another interesting 
feature of note in the chromosphere are filaments. These are regions of cool plasma 
suspended above the photosphere. They are seen as long, dark structures on disk 



in Figure |1.5| and as arcade-like features on the limb, where they are known as 
prominences. It is believed that the chromosphere is heated by some combination 
of conduction of heat from the hotter transition region and by the deposition of 
energy by waves. It is believed that acoustic waves are generated by turbulent 
motions in the photosphere which then form shocks as they propagate upwards 
through the chromosphere. With a plasma /3 of ~1, this highly ordered layer has 



very well defined structures, as Figure 1.6 clearly demonstrates (Carlsson et al. 
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Figure 1.7: A SOHO/EIT image of tlie Sun's transition region observed in tlie 
304 A (He ii) passband. 



1997 De Pontieu et ai. 2004) . 



1.1.3 The transition region 

Above the chromosphere, there is a very narrow layer only a few hundred km 



thick known as the transition region (Gallagher et a/. 1998 Mariska, 1993). The 



transition region is a highly dynamic interface between the chromosphere and the 
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the smooth temperature and density profiles found in the chromosphere, there 
are very steep gradients across the transition region, with values changing from 
~10^— 10^ K and ~10^^ — 10^ cm~^ across a height of the order 10^ km high. Con- 
duction (Fc) is highly sensitive to both the temperature (T) and the temperature 
gradient as it scales as = kqT^/'^'WT where kq is the coefficient of thermal con- 



ductivity (see Equation 1.11 for further details). The steep temperature gradient 
across the transition region means conduction is very efficient at transferring heat 
energy from the hot corona downwards, heating the upper layers of the transi- 
tion region as it does so. As a result, the upper transition region emits strongly 
in the ultraviolet (UV) and extreme ultraviolet (EUV) portion of the spectrum. 



The transition region appears brightest in areas called active regions (Figure 1.7). 
These are the EUV manifestations of the dense magnetic field regions observed 
as sunspots in the photosphere. 

1.1.4 The corona 

As we increase in height into the corona, increased temperatures of >10^ K and 
greatly reduced densities of ~10^ cm~'^ compared to the photosphere mean the 
corona emits strongly in EUV and in X-rays, particularly in active regions where 
large intricate loop systems are present. The complexity of these loop systems 
along with large magnetic field gradients make them the most likely place for 



solar fiares to occur ( Conlon et al. , 2008 Gallagher et al. , 2002a McAteer et al. 



2009). The tenuous, ambient emission in the visible corona however, can be very 
difficult to observe due to the stark intensity contrast between it and photo- 
spheric emission. Until the development of coronagraphs, the corona could only 
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Figure 1.8: The visible solar corona up to 6 Rq taken during the 2009 solar 
eclipse. The linear open field of polar coronal holes are clear, along with the closed 



loops of helmet streamers closer to the equator ( Druckmiiller et al. 2009 ) . 



be imaged in visible light during a solar eclipse. Figure 1.8 shows an image of 
the extended corona taken during the 2009 solar eclipse with the moon blocking 
out emission from the solar disk, making it possible to image the visible corona 
in striking detail. The linear features in this image result from emitting plasma 
flowing along magnetic field lines. 



There are two distinct features that can be seen in Figure |1.8[ At the poles 
there are regions of open field called coronal holes. These features exist perma- 
nently at the poles and intermittently closer to the equator. Open magnetic field 
lines stretching between the photosphere and interplanetary space have a large 
pressure gradient along them. This drives plasma out of the solar atmosphere 



into space in what is known as the solar wind (Altschuler et a/., 1972). As a 
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Figure 1.9: Extreme ultraviolet image of a set of flaring loops taken with TRACE. 



result, these regions tend to have a low density (~ 10^ cm ^; Wilhelm, 2006), 
showing up as dark regions in EUV images. 



Closer to the equator in Figure 1.8[ streamers and helmet streamers can be 
seen. These large, trans-equatorial loop systems have coronal plasma trapped 
along the field lines. High in the atmosphere however, the weakening field is 
dragged into a cusp shape by the solar wind. Helmet streamers are often found 
above active regions and prominences. Active regions are areas of enhanced 
magnetic activity. These tend to form in bands north and south of the equator 
that migrate towards the equator as the solar cycle progresses and are frequently 
associated with sunspots. The increased magnetic stresses in these regions result 



in large, impulsive releases of energy known as solar flares (^1.2). Flares are often 
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Babcock 



1961). Panel 



Figure 1.10: The solar dynamo (ail effect, adapted from; 
(a) shows the standard bipolar poloidal field of the Sun. Panel (b) shows differential 
rotation resulting in the wrapping of field lines around the Sun, leading to a toroidal 
topology (a effect). Panel (c) shows the increasing internal pressure inside flux 
tubes results in them rising through the photosphere effect). Panel (d) shows 
that magnetic reconnection with the bipolar field and with neighbouring magnetic 
loop systems results in a gradual cancelling of the toroidal field and a return to a 
poloidal topology. 



associated with coronal mass ejections ( ^1.3 Tousey et al.. 1973). These are the 
ejections of plasma, energy and magnetic field into interplanetary space. The 
remainder of the solar disk is what is known as quiet sun. Despite the name, 
these areas are far from quiet. It is believed that small scale activity is occurring 



constantly in the quiet sun by way of micro- and nano-fiares (Gallagher et al. 



1999 Klimchuk & CargiU, 2001 Parker, 1983 Schmelz et a/. 2009). 
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Active regions, and therefore solar flares, exist as a result of the emergence 
of areas of concentrated magnetic fleld through the photosphere. This occurs 



directly as a result of differential rotation (see Hoyng, 1990, for review). How- 
ever, the details of the solar dynamo remain elusive. The consensus within the 



community is that the solar cycle can be explained by the ail effect (Figure 1.10 



Babcock, 1961). The basis for this is the MHD mean field dynamo equation: 



d{B) 
dt 



V X {{v) X (B) + {v' X B') - r/V X (B)) 



;i.2) 



where B and v are the net magnetic field strength and fiow velocities of the 
large scale mean components and B' and v' refer to the small-scale turbulent 



motions (Charbonneau, 2005). r] is the magnetic diffusivity of the Sun. This 



essentially refers to the viscosity of the fluid. At the tachocline (base of the 



convection zone; Spiegel & Zahn, 1992), two characteristics affect the magnetic 



fleld. Firstly, r] is known to be small (Charbonneau, 2005) and so the cross 



product (advection) terms in Equation 1.2 dominate. Low diffusion means the 



magnetic fleld is "frozen-in" to the plasma and will move with both the large- 
scale, global plasma flows and the small turbulent motions. Secondly, above 
the tachocline, the convection zone rotates differentially (i.e. the Sun no longer 
rotates as a solid body above the tachocline). This forces the magnetic fleld 



to deviate from its initial poloidal state (Figure 1.10 1) and become wrapped 



up, generating a toroidal fleld (Figure l.lOj p). As regions of high magnetic fleld 
density develop, the internal pressure of the magnetic flux tube begins to increase 
until the pressure gradient is sufficient to cause the magnetic fiux tube to rise. The 
flux tube protrudes through the photosphere and can appear as a sunspot (Figure 
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TW^). Re gions where flux emergence is prominent can become active regions. 
As flux emerges through the photosphere, it is twisted and therefore contains 
non-potential energy. In an attempt to return to a force-free state, excess free 
energy can be released in solar flares and coronal mass ejections (Figure 1.10[ i). 



1.2 Solar Flares 



The study of solar flares began more than 150 years ago. Contrary to general 
knowledge, the first recording of the Sun-Earth connection was published by an 



Irishman, Colonel Edward Sabine (1852). Sabine hypothesised that the number 



of sunspots was connected with the level of auroral activity observed at Earth: 
"...it is quite conceivable that affections of the gaseous envelope of the Sun, or 
causes occasioning these affections, may give rise to sensible magnetical effects at 
the surface of our planet, without producing sensible thermic effects." The first 



image of a solar flare was published 7 years later by Carrington (1859), shown in 



Figure 1.11 The magnetic storm associated with this flare (now known to be a 
CME) was recorded all over the world. Articles were published in the Irish Times 
detailing the visibility of the aurora as far as Cork in the south of Ireland at 51° 



(Figure 1.12). Thus began the inquest into the cause of these transient and very 



intense brightenings (flares) and their "affection" to the Earth (CMEs). 

The mid-20th century brought the technology to record these brightenings 
using rockets and balloon flights. The first hard X-ray emission from a flare was 



recorded by Peterson & Winckler (1959). The development of orbiting satellites 



in the latter half of the century resulted in some ground breaking telescopes. 
Observatories of note include Skylab and SMM which were vital in furthering the 
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Figure 1.11: Diagram sketched by Carrington (1859). This is the first recorded 



image of a solar flare, indicated by the white regions marked A and B. 



^Strong et al. 


1999 


Sturrock 


1980) 



The Japanese Solar-A mission, later named Yohkoh (Ogawara et al. . 1991), was 
designed to study solar flares in the keV - GeV range. One of the most noteworthy 
revelations of the Yohkoh mission was the observation of what became known 



as the Masuda Flare (Masuda et al.. 1994). This paper was revolutionary as 



it proclaimed the reconnection region above the solar flare as the location of 
particle acceleration. While Yohkoh was concentrating on the soft and hard X- 
ray emission (SXR and HXR respectively), SOHO was broadening horizons with 
its suite of twelve instruments. While the primary focus of the SOHO mission 
was not the investigation of solar flares, its diverse set of instruments has ensured 
it played its part in the understanding of their behaviour. From measurements 
of the magnetic field out to white light observations of coronal mass ejections. 
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Figure 1.12: Irish Times article published on 8-Sept-1859 describing auroral ac- 
tivity in Cork, Ireland following the 1859 flare (courtesy of Peter T. Gallager). 



this observatory continues to facilitate the study of flares over much of their 
temperature ranges, length scales and wavelengths. While SOHO was unique in 
its vast range of instruments, RHESSl was revolutionary in a different sense. The 
primary goal of this mission is the investigation of particle acceleration and energy 
release during solar flares. RHESSI's unique ability to simultaneously produce 
high-resolution images and spectra in the X- to 7- ray regime is as a result of 
a combination of finely tuned rotating grids and cooled germanium detectors 



(see ^3.5 for further details). This instrument has facilitated the investigation 



of not only particle acceleration, but also the location of reconnection regions 



e.g. Krucker et al, 2009) and the motion of this reconnection region with the 



evolution of flares (e.g. Grigis & Benz, 2005). It also facilitated a statistical study 
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Figure 1.13: Frequency distribution of the thermal energy of microflares (RH), 
nanoflares (TA, TP, EB) and active region brightenings (SS). They have been fitted 



with a power-law function of index 2 (dotted line) . ( Hannah et al. , 2008 ) 



of 25,000 microflares (Cliriste et al. 2008 Hannali et al, 2008). 



The aforementioned investigation by Hannali et al (2008) reveals one of the 



many important reasons for studying solar flares. The frequency distribution 



of flares, shown in Figure 1.13, shows that, for the most part, the frequency 



distribution of flare energy follows a power-law function of index 2.0 ± 0.4. This 
is important as it is believed that this distribution of energies may be what is 
driving the heating of the corona. It is believed that the energy of the many 
nano flares that occur may be sufficient to supply enough energy to the corona to 



19 



1.2 Solar Flares 



Table 1.1: GOES classification scheme for solar flares. 



GOES class Minimum flux 
[Watts m-^] 

X 10=^ 

M 10"^ 

C 10-6 

B 10-^ 

A 10-8 



achieve the temperatures we observe. The critical value is the power-law index. 
If the index is less than 2 then the energy in micro- and nano-fiares is insufficient. 
However if it is >2 then "minor" flares have the potential to be the heating 
mechanism for the solar corona. The magnitude of flares are classified according 
to the maximum fiux observed by the GOES Satellites. A fiare is assigned a class 



based on a logarithmic scaling shown in Table |1.1[ The highest classification is 
"X-class" with a fiux of lO""^ W m"^. For each dex from lO""^ to lO^^ W m"^, 
the classes X, M, C, B and A are assigned. A secondary classification is used to 
indicate the level within each class. E.g. an M3.2 fiare has a maximum GOES 



fiux of 3.2x10-^ W m-2. See Qfor further details. 

Solar fiares are highly complex events that span temperature ranges over four 
orders of magnitude (10^ — 10^ K) and energy ranges over six orders of magnitude 
(~keV - GeV). These impulsive bursts of energy are some of the most powerful 
events in the solar system, releasing up to 10'^^ ergs (10^^ J) in tens of minutes 



(Emslie et ai, 2004). There is generally considered to be two classifications of 
solar fiares: compact and eruptive. In compact fiares little or no loss of material 
occurs, while in eruptive fiares, there is generally an associated coronal mass 
ejection that carries away plasma and magnetic field from the system. During all 
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Figure 1.14: Schematic diagram showing the processes involved in the standard 
model for solar flares that has been adapted from Dennis &: Schwartz (1989). The 



main processes in the evolution of a solar flare are labelled, (a) is the energy 
deposition region, (5) shows the path of propagating electrons towards the dense 
chromosphere (c). Heating of the chromosphere results in chromospheric evapora- 
tion {d) which fills the loops with hot, SXR emitting plasma (e). Direct heating of 
the looptop can result in large pressure gradients that drive conductive fronts (/). 
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types of flares, tfiere exists two phases: tlie impulsive phase during which plasma 
is heated to high temperatures and the decay phase during which the flare cools 



back to equilibrium (Dennis & Schwartz, 1989). 



Figure 1.14 gives an overview of the processes involved in the standard model 
for solar flares. During the impulsive (or rise) phase of a flare, energy is deposited 



into the loop (marked a on Figure 1.14). It is generally accepted that this is 
driven by magnetic reconnection ( ^1.2.1 ) and follows the thick target model of 
Brown ( |1971 ). During magnetic reconnection, energy stored in the loops (e.g. 



as twist) is released and is used to accelerate coronal particles. These particles 



propagate down magnetic field lines towards the chromosphere (Figure 1.14?). 
The sudden increase in density at the chromosphere (Figure 1.14[ :) results in the 
beam particles interacting with ambient particles in the chromosphere through 
coulomb collisions, resulting in the emission of Bremsstrahlung HXR radiation. 
The energy transferred by the beam particles is absorbed by the chromosphere 
and, where possible, radiated away. However, should the rate of energy deposition 
be too great for the chromosphere to efficiently radiate, pressure gradients can 
build up in the plasma. This results in the expansion of the plasma into the loop 



(Figure 1.14i), filling the loop with hot, SXR and EUV emitting plasma (Figure 
1.14^ ) in a process known as chromospheric evaporation. 



Fisher et al.\ ( [l984[ [l985a||bP were among the first to study chromospheric 



evaporation in detail from a theoretical perspective. By running simulations to 
replicate the behaviour of the chromosphere to different fluxes of non-thermal 
beams of electrons, they established a threshold for the flux of non-thermal elec- 
trons required to drive explosive chromospheric evaporation. It was found that 
fluxes of non-thermal particles of less than 10^° ergs cm~^ s~^ do not generate 



22 



1.2 Solar Flares 



sufficient pressure gradients to drive explosive evaporation and the velocities ob- 
served were of less than 20 km s"^. This is what is now classified as gentle 
chromospheric evaporation. Gentle evaporation can also be driven by pressure 
gradients that result from the direct heating of the looptop, driving conduction 



fronts towards the chromosphere (Figure 1.14 f; Antiochos & Sturrock, 1978). 
This is believed to occur early in the decay phase when conduction is most ef- 
ficient due to high plasma temperatures. Explosive chromospheric evaporation 
results from fiuxes of greater than ~3x 10^*^ ergs cm~^ s~^. This drives upfiows of 
hundreds of km s~^ due to the large pressure difference between the heated mate- 
rial and the tenuous corona. Low velocity downflows (tens of km s~^) predicted 



by their models were later observed by Zarro & Canfield (1989). Although the 
velocity of the downfiows are orders of magnitude smaller than the upfiows, the 
components and densities of the chromosphere and corona are such that momen- 



tum within the system is conserved (Canfield et al, 1987 Teriaca et al., 2006). 



The velocities expected from chromospheric evaporation were calculated in 



Fisher et al. (1984). It begins with the equation of motion in one dimension: 



Dv 



:i.3) 



where v is velocity, P is pressure and 



D _ f d d 



;i.4) 



is the so-called convective derivative in one dimension. This takes account of the 
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Figure 1.15: Simulation results from Fisher et al. (1985c) showing temperature 



(T), number density (n), pressure (P) and velocity (v) as a function of column 
depth measured from the loop apex (N) at significant times in the simulations 
(i.e. 2 and 4 seconds after the electron beam is imposed on the system). The 
flux of non-thermal electrons at energies less than 20 keV (F20) are shown for the 
explosive (right) and gentle (left) chromospheric evaporation profiles and preflare 
profiles (dashed curves). Positive velocities correspond to downflows (redshifts) 
and negative velocities correspond to upflows (blueshifts) . 



rate of change of the velocity as it moves in space (z) and time (t). Since the mass 
density, p can be expressed as the product of the mean mass per hydrogen nucleus, 



m and number density of hydrogen nucleii, n, Equation 1.3 can be rewritten as 



' dv dv\ 1 dP 
dt dz ) n dz 



;i.5) 
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If we ignore the time derivative, assuming a constant velocity evolution at any 
given height in the loop. Assuming that the coronal protons and electrons have 
the same temperature, the pressure can be defined as Ptot = Pion + Peiectron and 
since it can be assumed that Pion = Peiectron = nksT in the corona, we can write 
Ptot = ^nksT (Antiochos & Sturrock, 1976 Kivelson & Russell, 1995 Krall et al. 



1998). Therefore, Equation 1.5 reduces to 



dv J ^1 dn 
mv— = —zkbT- 



dz 



n dz 



;i.6) 



Integrating between the chromosphere {ch) and the front of the expanding mate- 
rial in the corona (co), we obtain: 



miVco - Vch 



2/cbT [ln{nch) - ln{nco)] 



;i.7) 



If we assume that once the fiare occurs, velocities in the chromosphere are negli- 
gible, rearranging gives us 



4fcr^^ / rich 
m \ rirn 



1/2 



1 /2 

This can be written in terms of the sound speed Cg = {•ykBT/m) where 7 is 
the heat capacity ratio. 



In 

5 / KUrn 



1/2 



;i.9) 



The plasma velocity in the corona works out to be approximately twice the coronal 
sound speed. Changing the density ratio (rich/nco) from e.g. 10^ to 10^ has 
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little effect on this with Vco changing from ~ 2.35cs to ~ 2.88cs. [Milligan et al. 



(2006a) recorded upflow velocities of 230 ±38 km s ^ and simultaneous downflow 



velocities of 36 ± 16 and 43 ± 22 km s ^ at chromospheric and transition region 



temperatures respectively. Raftery et al. (2009) observed maximum up flows of 



95 ± 10 km s^^ and simultaneous downflows of 23 ± 10 km s^^. Considering the 
model adopts a constant velocity approach, thus ignoring any "start-up" time 
required to accelerate the plasma from rest, this result is quite reasonable. 

The rate at which the loop fills as a result of chromospheric evaporation was 



found to be closely correlated to the HXR flux (Neupert, 1968). This effect 



named after its discoverer relates the flux [F) of hard and soft X-rays as follows: 

dFsxR 



dt 



oc Fhxr (1.10) 



This means that the derivative of an SXR light curve can be used to approximate 
the size and duration of an associated HXR burst. This is believed to stem 
from the time difference between the heating of chromospheric plasma during 
chromospheric evaporation and the filling of the loop with evaporated plasma. 
It takes a finite time for the heated plasma to rise into the loop and reach a 
sufficient density to begin to emit in SXRs. This suggests that the energetic 
electrons responsible for the HXR flux are also responsible for chromospheric 
heating of the flare loop. In reality, multiple-loop systems and other heating 



mechanisms result in slight deviations from this relationship (e.g. Dennis & Zarro 



1993). However, it can be a very useful technique for estimating the duration of 
a HXR burst. The flux of HXRs are often determined using observations from 
the RHESSI spacecraft. This spacecraft however, is subject to regular night-time 
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passes, or Earth occultations due to its low Earth orbit. As the satellite orbits 
the Earth, it can pass into the night-side of the Sun-Earth line. RHESSl also 
passes through the South Atlantic Anomaly (SAA). This region is where the 
Earth's inner Van Allen belt is at its closest. Since the Van Allen belt is aligned 
with the planet's magnetic axis and not its rotational axis, the belt is closest 
to the Earth over the south Atlantic ocean. As a satellite passes through this 
region, it will be exposed to strong radiation, contaminating any observations it 
is taking. The frequency and duration of the Earth occultations and SAA passes 
vary throughout the year. However, with a 90 minute orbit and up to 20 minutes 
of interference, they can still have a dramatic effect on the number of events that 
are observed unhindered. However, SXR observations from the GOES satellite 
in 0.5-4 A and 1-8 A range are taken every 3 seconds with no interruptions. 
Therefore, when HXR observations from RHESSl are contaminated, the HXR 
burst can be approximated by the derivative of the GOES SXR lightcurve. 

The peak of the SXR flux often occurs simultaneously with the peak of emis- 
sion measure. This is not surprising, since the emission will have its highest 
intensity when the loop is at its densest. This occurs just after the time of max- 
imum temperature. Like the Neupert effect, this delay IS cLS cL result of the loop 
filling time. Around the time of maximum temperature/emission measure con- 
duction is found to be the dominant heat transfer process. Spitzer conductivity 



(Spitzer, 1956) is widely accepted as the form for conduction in solar plasmas: 

Fc = KqT^/^VT (1.11) 
where k,q is the Spitzer coefficient. Spitzer conductivity is not the most suitable 
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Temperature (K) 

Figure 1.16: The radiative loss function {Er{T)) as a function of temperature 
for coronal abundances (upper curve) and photospheric abundances (lower curve) 



(Dere et al. 2009). 



form for solar flares, considering it was devised for non-magnetised plasma in 
equilibrium. However, in the absence of a more appropriate relation, it is widely 
accepted within the solar physics community. During the very early decay phase 
while temperatures are still high, conduction is very efficient in redistributing 
heat throughout a system, although it does not actually remove heat from the 
system. If for example, the corona is studied as an isolated system, it can be 
heated by the conduction of heat into the corona and cooled by the removal of 
heat to the transition region. During the early decay phase while temperatures 
are high and the temperature gradient is steep, conduction is very efficient. It is 
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the primary mechanism for heat loss in the corona for the first ~10^ seconds. As 
the loop thermalises and temperatures begin to fall, its efficiency is reduced. As 
the temperature approaches ~ 10^ K, the efficiency of heat loss by radiation is 



maximised (Figure 1.16). The radiative loss rate is given by: 



dt 



;i.l2) 



where A(T) is the optically thin radiative loss function, as shown in Figure 1.16 
making radiation an effective cooling mechanism for ~ 10^ seconds in the mid- 



to late decay phase (Culhane et al. 1970 Raftery et al. 2009) 



These processes describe the standard model for solar flares. While the stan- 
dard model is widely accepted, it is also known to include some serious flaws, such 
as the number problem. It has been shown that the number of electrons required 
to produce the Bremsstrahlung emission by the thick target model is approxi- 



mately 10^^ — 10^^ electrons per second (Holman et a/., 2003). With 10^ cm~^ as 
an average electron density in the corona, a cubic volume of 10^'^ cm^ is required 
to be evacuated into the chromosphere every second. Considering reconnection 



is believed to take place across length scales of the order of meters (see {1.2.1 
for further details), and an average coronal loop is ~ 10^° cm long, this proves 
to be a signiflcant problem in the standard thick-target model. The volume of 
the corona surrounding the entire loop would have to be accelerated and replen- 



ished every second. The concept of return currents (e.g. Benz, 2008) has been 



presented as a possible mechanism of returning accelerated electrons back to the 



corona. 



Fletcher & Hudson (2008) have also presented a possible solution to the 



number problem by moving the site of particle acceleration to the chromosphere. 
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where the number density of electrons is significantly higher. 
1.2.1 Magnet ohydrodynamics 

Magnetic reconnection is widely believed to be the driving force behind solar 
flares. The interaction of solar magnetic field with plasma can be understood 
using the equations of magnetohydrodynamics (MHD). These set of equations 
describe the behaviour of the electric field E, the magnetic field B, current density 
j and plasma velocity v. 

Gauss' law relates the distribution of electric charge Pq to the resulting electric 
field by way of the permittivity of free space, eo: 

V-E=^ (1.13) 
Co 

Faraday's law of electromagnetic induction states that changing a magnetic 
field in time will induce an electric field: 



Gauss' law for magnetism states that there are no magnetic monopoles: 



Ampere's law states that either a current or a time varying electric field will 
produce a magnetic field. 



V X E = - 



(9B 

~dt 



(1.14) 



V-B = 



(1.15) 



V X B = 



1 aE 



(1.16) 
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where /io is the permeability of free space and is related to the speed of light 
by = l//ioeo- In addition to Maxwell's equations (Equations 1.13 to 1.16), 



Ohm's law states that the current density is related to both the electric field of 
the plasma and the motion of the plasma at velocity v relative to the magnetic 
field. 

j = (T(E + vxB) (1.17) 

where a is electrical conductivity. 

MHD also incorporates the equations of fluid dynamics for a plasma with 
density p and pressure P. The equation of motion for a parcel of fluid states that 
the rate of change of fluid velocity is governed by the pressure gradients: 



p— = -VP-p^? + (jxB) (1.18) 



The mass continuity equation states that the rate at which mass enters a system 
is equal to the rate at which mass leaves the system: 



^ + pV-v = (1.19) 



The energy equation describes the rate of change of energy {E) in a plasma as 
a result of various heat sources and sinks, such as the heating rate (Q), the 
divergence of the conductive flux (Fc) and the radiative loss rate {E^). 



— = Q-^.F,-Er (1.20) 
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Assuming non-relativistic velocities, Equation |1 . 16| can be approximated as: 



V X B = /ioj 



;i.21) 



Substituting Equations 1.17 and 1.21 into Equation 1.14 gives: 



(9B / 1 

— = Vx (vxB V X B 

at \ ajiQ 



;i.22) 



Utilising the identity: 



V X (V X B) = V (V • B) - (V ■ V) B, 



1.23) 



and recalling Equation 1.15, Equation 1.22 can be rewritten as what is known as 



the induction equation: 



~dt 



V X (v X B) + r^V^B 



;i.24) 



where magnetic diffusivity is given hj r] = l/afiQ. The two terms on the right 



hand side of Equation |1.24| describe advection and diffusion respectively. The 
ratio of these terms is known as the magnetic Reynolds number, Rm- 



Rr. 



V X (v X B) vL 



r?V2B 



7] 



;i.25) 



For a perfectly conducting plasma, i.e. a — )■ oo and ?7 — )■ 0, the change in the 
magnetic field is completely dominated by advection and the field is carried along 
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by the flowing plasma. Thus Equation 1.24 can be approximated by: 

— ^ V X (v X B) (1.26) 

and the fields are said to be "frozen-in" to the plasma. Although the plasma is 
not perfectly conducting, in most astrophysical plasmas, Rm ^ 1- For example, 
in the solar corona, rj ^ 1 m^ s^^, L ^ 10^ m, and v ~ 10'^ m s~^ gives a value of 
Rm ~ 10^. The primary exception to this is in the presence of large magnetic field 
gradients. Under these circumstances, the changing field can be approximated 
as: 

— ^ r/V^B (1.27) 

Here, Rm ^ 1 and diffusion dominates over advection, allowing the magnetic field 
to slip through the plasma. This is one of the conditions required for magnetic 
reconnection. A simple schematic of the processes involved in magnetic reconnec- 



tion is shown in Figure 1.17 When regions of opposite polarity magnetic fiux are 
in close proximity to each other, large magnetic field gradients will be established 
as the value of the field goes from positive in one region, to zero at the neutral 
line to negative at the opposing fiux. The large gradients in the magnetic field 
result in Rm ^ 1. This allows the magnetic field to diffuse through the plasma 
and reconnect, resulting in an energetically more favourable topology. In doing 
so, the plasma is ejected perpendicular to the in-flowing flelds which creates a 
drop in pressure. This in turn pulls more plasma and magnetic fleld into the 
diffusion region, repeating the process. 



Since Equation 1.24 can be approximated as Equation 1.27 in the diffusion 
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Figure 1.17: Schematic diagram of magnetic reconnection. Field lines approach 
the diffusion region where they reconnect. The resulting topology is energetically 
more stable. 



region, the timescales can be estimated as 



:i.28) 



The timescales for reconnection are known to be on the order of seconds. There- 



fore, to satisfy Equation 1.28 , there is a requirement for very short length scales 
(on the order of meters). The reconnection model presented by Sweet and Parker 



(Parker, 1957 Sweet, 1958) utilised a thin current sheet with length of the order 



of a coronal loop along which reconnection can take place. With L ~ 1 Mm, the 
reconnection rate is of the order of 10^° seconds, or roughly a billion years! Since 



reconnection on the Sun requires timescales on the order of seconds, Petschek 



(1964) proposed a system where the diffusion region is of the order of meters, a 



fraction of a scale length. This means that direct observational evidence of mag- 
netic reconnection is not yet possible. The low density of the diffusion region, 
combined with its small length scale is beyond the capability of current instru- 
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ments. However, indirect measurements of this phenomenon are widespread in 



literature, from solar jets (Bain & Fletcher, 2009) to magnetospheres (Slavin 



et al. 2009) and in laboratory plasmas (e.g. Cothran et al, 2003). 



1.2.2 Hydrodynamic modelling 



While the use of MHD is essential to the understanding of the behaviour of field 
and plasma, to quantitatively solve the equations of MHD in 3 dimensions is a 
very complex task that requires a huge amount of computing power. As such, 
approximations are often made to these equations in 0- and 1- dimensional codes. 
Even the computation load of 1-D codes is not trivial, however the use of 0- 
D hydrodynamic simulations are very quick (~seconds on a personal computer 
versus hours to weeks on a cluster for 1-D). The term "0-D" stems from the fact 
that these models sacrifice spatial resolution in favour of efficiency by assuming 
that any energy in a coronal loop undergoing magnetic reconnection is uniformly 
distributed in the corona. This is a reasonable representation since temperature, 
density and pressure are approximately uniform along the magnetic field, with 
the exception of the steep gradients in the transition region. It has been shown 
that these codes are very robust in their responses and compare very well to more 



detailed 1-D codes (see e.g. Klimchuk et al, 2008, for details of comparison). 
In general, hydrodynamic models begin with the 1-D time dependent version 



of the energy equation, given in Equation 1.20 



:i.29) 



This equation states that the rate of change in energy of a system (E) is balanced 
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by the heating rate (Q), and the loss rates by conduction radiation (two right 
hand terms which are defined in Equations 



1.11 



and 1.12) 



1.2.2.1 The Cargill model 



The simple model presented by Cargill (1993, 1994) is an effective method for 



making a first approximation for the cooling timescale of a fiare. Beginning with 



the basic energy equation |1.29| and assuming velocities generated by evaporation 
are much smaller than the sound speed, we can substitute for the different terms. 
Energy is given hj 3/2nkBT, conduction is approximated from the Spitzer formula 
LTT|as2/7/toT^/VL^ 



given in Equation 



3 d [nksT] 2 T^' 

2^^ = Q-7-o^-nA(T) 



;i.30) 



The temperature evolution of a system can be obtained from Equation 1.30 



Following work laid out in Antiochos & Sturrock (1978) and Antiochos (1980) 



the assumption is made that at any given time, cooling is done by only a single 
mechanism, i.e. in the early decay phase radiation is ignored and in the late decay 
phase, conduction is ignored. For example, the conductive cooling time can be 
calculated by setting Q and the radiative loss rate to 0. This leaves 



3d[nkBT] _ 2 

2 dt ""7"°^^ 



;i.31) 



Rearranging this equation to give 



3nk 



To 



-2 /to 
7 L2 



dt 



1.32) 



to 
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we can integrate to get the temperature evolution for the conductive phase of a 
flare: 



T{t) = To 1 



t 



Tr. 



-2/7 



;i.33) 



The temperature evolution for the radiation phase is obtained in a similar fashion 
and is given by 



Tit) = To 1 



3 t 



2r, 



;i.34) 



In Equations [L33] and [L34J the parameters Tc and Tr correspond to the conductive 
and radiative cooling timescales: 



and 



3nkBT 



;i.35) 



1.36) 



' n^KiT) 
respectively. 

The time and temperature at which the cooling mechanism dominance changes, 
and respectively, can then be expressed as: 



Co 



' CO 



7/12 



1 



;i.37) 



and 



Tr 



-1/6 



;i.38) 



T(t) is shown in the top panel of Figure 1.18 with T*(r*) highlighted as the break 



between conductive and radiative cooling. Note that this occurs in conjunction 



with the peak density (middle panel. Figure 1.18). The bottom panel shows 
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Figure 1.18: Sample plots from Cargill (1994) showing the cooling profile of a 
flare of density 5 x 10^ cm~^, initial temperature of 10^ K and a half length of 
2.5 X 10^ cm. The top two panels show the temperature and density evolution of 
the flare. The cutoff time and temperature {t^:,T^:) are shown in the top panel. 
The bottom panel shows the instantaneous radiative and conductive cooling times. 
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the ratio of Tr/Tc- The radiative cooling time is significantly longer than the 
conductive cooling time during the first ~ 1000 seconds of the flare. After this 
time, the radiative cooling time becomes shorter than the conductive timescale, 
making radiation the dominant cooling mechanism. 

One assumption of the Cargill model is that of independent cooling mech- 
anisms. It is assumed that while conduction is efficiently removing heat from 
the system, radiation is negligible and is therefore ignored, and vice versa. This 
is a reasonable thing to assume at the beginning and end of the decay phase. 
However, there is a time when both conduction and radiation have significant 
contribution to the removal of energy from the system. Therefore, during this 
time, a single loss mechanism will underestimate the energy removed from the 
system and potentially result in a higher temperature or longer cooling time. 

1.2.2.2 Enthalpy Based Thermal Evolution of Loops (EBTEL) 



The Enthalpy Based Thermal Evolution of Loops model (EBTEL; Klimchuk 



et ai, 2008) accounts for both conductive and radiative losses throughout the 
lifetime of the flare, thus eliminating some of the limitations of the Cargill model. 
EBTEL takes explicit account of the important role of enthalpy in the energetics 
of evolving loops. The basic assumption of the EBTEL model is that n, P and T 
in the corona can be represented by spatial averages since they generally vary by 
less than a factor of 3 through the corona. The base of the corona is defined to 
be the point at which conduction switches from a heating term in the transition 
region to a cooling term in the corona. This is based on the categorisation of the 
enthalpy flux. As the heating rate increases during, say, a flare, the chromosphere 
is unable to radiate all the absorbed energy, and so plasma rises into the corona. 
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Therefore, an excess in heat flux can be associated with the impulsive phase of 
a flare. Conversely, when the heating rate decreases, the coronal temperature 
begins to fall, creating pressure gradients in the loop and causing plasma to flow 
towards the footpoints. Thus a deficit of heat flux is associated with plasma 
cooling. 

By incorporating the enthalpy H of the system, we take account of the work 
done on and by the system: 



where Q is the combined direct and non-thermal heating rate, s is distance along 
the loop, Kq is the Spitzer conductivity coefficient, T is the temperature, is the 
electron density and A(T) is the radiative loss function. Considering the enthalpy 
can be written the sum of the internal energy [E) and the work done {PV), we 
can write 




(1.39) 




(1.40) 



Substituting for the thermal and kinetic energies: 



^ 3„ 1 , 



(1.41) 



we can write 



sap 1 d (pv^) 
2~dt ^ 2 dt 



Q- 



d{Pv) 3dPv ld{pv^)v 



ds 2 ds 2 ds 




(1.42) 
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Assuming a subsonic flow, the kinetic terms can be neglected giving 



3dP bdPv d 



2 dt ^ 2 ds ds 

EBTEL makes the assumption that the energy equation can be solved inde- 
pendently for the corona and the transition region. If we designate the subscript 
"0" to define the base of the corona (see above for definition) and assume a coro- 
nal half length of L and a transition region half length of /, integrating Equation 



1.43 1 over these regions gives: 

3 dP - 5 

-L— ^QL+ -Povo + F,, - Fr^ (1.44) 
3 dP - 5 

-/— ^Ql- -Povo - - F,,, (1.45) 

for the corona and transition region respectively where the overbar refers to spa- 
tially averaged values along the appropriate region. Fr refers to the radiative 
cooling rate per unit area and F^ to the conductive flux. Note the difference 
in signs for the conduction term. This implies that heat conducted out of the 
corona is conducted into the transition region. It is assumed that the enthalpy 
flux is negligible at the base of the transition region (i.e. the majority of the 
heat flux energy is distributed through the upper layers, heating each consecu- 
tive layer). Since IdPtr/dt <^ LdPc/dt, it is also assumed that pressure in the 



transition region is constant and we can write Equation |1.45| as 

5 



^PoVo = -F^-Fr,^ (1.46) 



41 



1.2 Solar Flares 



Equation 1.46 directly describes the directionality of the heat flux. When \Fcg\ > 



Fr^^ then transition region radiation is insufficient to remove the heat and it is 
conducted into the corona. If, however \Fcq\ < F^.^^, then the transition region is 
radiating efficiently enough to draw heat from the corona. This, combined with 



Equation 1.44 leads to 



dP 
~dt 



Q-- {Fr^ + 



1.47) 



One important assumption made when solving Equation 1.47 is 



C'rad 



(1.48) 



at all times, where Crad is a constant of the system. This is a reasonable approxi- 
mation during the decay phase of the flare since the long decay timescales mean 
the system does not deviate far from static equilibrium at any given time. How- 
ever, during the impulsive phase this is not necessarily valid. The quickly chang- 
ing density of the system can invalidate this approximation. Despite comparisons 
to 1-D hydro models that suggest that this ratio does not have a significant effect 



on the resulting parameters (Klimchuk et a/., 2008), an investigation carried out 



by Adamakis et al. (2008) revealed that the resulting output parameters are in 



fact sensitive to this ratio. 



Since F^^ = nlK{T)L and from Equation 1.48, F^^^ can be written in terms of 
Fr^, the coronal pressure can be written as: 



dP 
~dt 



nlk{f)L 

Q (1 + Crad) 



(1.49) 
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The conservation of mass requires that the enthalpy flux through the transition 
region be approximately constant. Since the total mass contained within a length 
L will change with the evaporation of material into the corona, this must equal 
the electron flux from the transition region to the corona (i.e. through the base 
of the corona), we can write: 

=nv^Jo (1.50) 
dt ^ ^ 



Combining this with Equation 1.46 and considering P = 2nkBT, we can describe 
the rate of change of density as 

dfl 1 

Finally, the temperature evolution can be obtained from the ideal gas law as: 



T(if-^f] (1.52) 



dt \P dt n dt 
The evolution of the average coronal pressure, density and temperature can there- 



fore be obtained as a function of time from Equations 1.49 1.51 and 1.52 



Figure |1.19| shows the time evolution of the three parameters: temperature, 
density and pressure of a nanofiare heated impulsively for 500 seconds. The com- 
parison of the EBTEL simulation (solid line) to a 1-D simulation (dashed line) 
shows very good agreement between the two models, despite the significantly sim- 
pler approach taken with EBTEL. It is also interesting to note that unlike Figure 



1.18 there are no discontinuities in the EBTEL evolution. This is a significant 



improvement on the Cargill model. The efficiency of combining losses by con- 
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Figure 1.19: Sample plots from Klimchuk et al. (2008) showing tiie coronal av- 



eraged temperature, density and pressure of an impulsively heated nanoflare. The 
solid lines are the EBTEL simulation while the dashed lines are from a more com- 



plex 1-D simulation for comparison (Antiochos &; Klimchuk, 1991). 
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duction and radiation throughout the flare are evident from the sharp reduction 
in temperature in the first ~600 seconds of the decay phase. Note the change in 
slope of the temperature evolution around the time of the density maximum. As 
the temperature falls, conduction is no longer as significant. Simultaneously, the 
density peaks and radiation becomes the dominant loss mechanism. 

The basic model scenario is such that the corona is heated primarily by con- 
duction fronts (direct heating). However, it is also possible to use EBTEL to 
approximate the combined effects of direct heating and heating by a non-thermal 
electron beam (non-thermal heating). The inclusion of the non-thermal particle 
heating is not thorough but is sufficient to give a reasonable estimate of the effect 
of a beam. It is assumed that any accelerated particles originate from within the 
system, specifically the corona. Thus, the number density of the entire loop does 
not change. It is also assumed that the accelerated particles can stream freely 
and not interact with any particles until they reach the chromosphere. The final, 
and most concerning approximation, is that all energy found in the beam is used 
to evaporate chromospheric plasma upwards into the loop. This is a concern 
because it is well known that propagating particles result in the expansion of 



plasma downwards as well as upwards (e.g. Fisher et al, 1984 Milligan et al. 



2006a ; Raftery et al. , 2009 ) . Chromospheric emissions are another result of accel- 



erated particles interacting with the chromosphere, albeit a tiny fraction (10 ^) 



of their energy (Dennis, 2007). 

The EBTEL model can accommodate any combination of direct and non- 
thermal heating functions, though it cannot accommodate purely non-thermal 
heating. Following the procedure above, the average temperature, density and 
pressure in the corona as a function of time are calculated. A secondary, but 
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important result is the evolution of the conductive and radiative losses in the 
coronal portion of the loop, allowing the user assess to the efficiency of the cool- 
ing mechanisms for the duration of the event. Other model options, including 
the differential emission measure of the loop may also be calculated. However, 
the inclusion of a non-thermal beam results in ambiguous transition region dif- 
ferential emission measure values, as the dependence on deposition depth are not 
considered. 



1.3 Eruptive flares and CMEs 

Solar flares are well known to be associated with coronal mass ejections (CMEs). 
CMEs are ejections of material, magnetic field and energy from the solar corona. 
They are generally bulbous structures threaded with magnetic field that grow 
radially as they propagate away from the Sun. CMEs occur across many differ- 
ent length- and time-scales. CMEs appear on ever increasing lengthscales as the 
propagate away from the Sun. They are most frequently observed using white 
light emission by imaging the Thompson scattered light of the K-corona (pho- 
tospheric emission scattered off free electrons inside the CME). They generally 
have a three part structure: a bright leading edge, a dark sparse cavity and a 



bright dense core. The three components are highlighted in Figure 1.20 



The association of CMEs with solar flares is widely known but not well un- 



derstood. Gosling (1993) postulated that geomagnetic storms are produced by 
CMEs and not, as previously believed, by solar flares. This declaration, dubbed 
"The Solar Flare Myth" , led to the misunderstanding that solar flares were not an 
important aspect of solar physics research as they had no effect on life on Earth. 
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Figure 1.20: A coronagraph image of a typical CME. Tlie briglit front, dark 
cavity and briglit core are highlighted. 



This belief divided the community and was contested on numerous occasions 
(e.g. 



Hudson et al. 



1995 



Reames 



1995 



Svestka 2001). The significance of solar 



flares has since been restored and was summarised nicely by Svestka (2001): "It 



is misleading to claim that flares are not important in solar-terrestrial relations. 
Although they do not cause the CME phenomenon that propagates from the Sun 
eventually hitting the Earth, they are excellent indicators of coronal storms and 
actually indicate the strongest, fastest and most important storms." Since then, 
the ideas connecting flares and CMEs have left the "cause and effect" paradigm, 
and it is now widely believed that eruptive flares and CMEs both result from the 



same driving mechanism (Zhang et al.. 2001) 



The combined work of Carmichael (1964), Sturrock (1966), Hirayama (1974) 



and Kopp & Pneuman (1976), known as the CSHKP model, began the unification 



of the fiare-CME theories. This model began with the proposition that low lying 
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Figure 1.21: The CSHKP model, adapted from Magara et al. (1996) 



coronal loops and open overlying field lines could create an upside-down Y-type 



magnetic topology (Figure 1.21, left panel). Following this, shearing motions may 



trigger a tearing instability near the neutral line located above the loop (within 
the diffusion region). This can result in magnetic reconnection and the acceler- 
ation of particles which propagate to lower altitudes and result in a solar flare 
beneath the reconnection region. The sling shot effect of the reconnected fields 
can result in shocks perpendicular to the incoming field. Accelerated particles 
can cause HXR footpoints and chromospheric evaporation of hot material into 
the loop to create SXR emitting loops beneath the reconnection region as with 
confined flares. Thin target emission is believed to be the cause of short lived 
HXR sources that form at the top of the loop/base of the reconnection region 



and may be the signature of the reconnection site (e.g. Krucker et al. , 2009 ). The 



reconnection and reorganisation of the field results in the next set of field lines 
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being drawn into the diffusion region to be reconnected. Thus a run-away process 
begins. With each set of reconnecting field lines, the CME is freed some more: 
as more field reconnects there is less overlying field preventing the propagation of 
the CME and it can therefore rise faster. This in turn results in faster reconnec- 
tion and an explosive burst of acceleration occurs. This acceleration burst has 
been found to be closely associated with the hard X-ray profile of the underlying 



eruptive flare (Temmer et a/., 2008). As successive field lines are reconnected, 



the post flare loops appear to rise as they are heated into the passbands of our 
instruments. Thus, as the CME rises, the post flare arcade is also observed to 
rise beneath the CME. 

1.3.1 CME initiation 

There have been many models put forward to try to simulate the behaviour of 
eruptive flares and CMEs. A successful model must be able to reproduce the 
observations as closely as possible. However, with ever evolving and improving 
observations, it is difficult to simulate every aspect of the system. This is espe- 
cially difficult considering the large energy and length scales involved. There have 



been many reviews on this topic. For example Forbes et al. (2006) contains an 



overview of the entire research field, including CME initiation, propagation, struc- 



ture, modelling and shock formation. The review by Klimchuk (2001) contains 



a description of CME initiation models using illustrations analogous to simple 
mechanical systems such as springs, pulleys and bombs. While this is instructive, 
it can also be misleading in how simplified the models appear to be. The author 



believes that Moore & Sterling (2006) strikes a good balance between conceptual 
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progression and physical description. Moore & Sterling (2006) describe three in- 
dependent methods of initiating a CME: internal tether cutting, external tether 
cutting and ideal MHD instability. 

In the internal tether-cutting case the initial topology is as shown in the top 



panel of Figure 1.22 This 2-D sketch implies a central sheared core (represented 
by the innermost loop of the quadrupole), tethered by a central arcade. The 
presence of neighbouring arcades is also a possibility but not a necessity for this 
model. Before eruption, the central arcade is in force free equilibrium and no 
current sheet exists between it and the overlying field. However, a current sheet 
does exist between the legs of the arcade as a result of their slow shearing due to 
photospheric motions. When this current sheet becomes sufficiently thin enough 



to allow reconnection across it (as in Figure 1.22 middle panel), a "run-away 
tether-cutting" process begins. The field lines above the reconnection site are 
now no longer tethered to the photosphere and begin to erupt upwards while the 
field lines beneath the reconnection site become the site of a solar fiare. As the 
reconnection progresses, the plasmoid above the reconnection site slowly erupts 
upwards and compresses the null point, forming a current sheet with the overlying 
field. This will result in explosive breakout reconnection above the plasmoid 
which reconnects field lines into the neighbouring arcades, both heating the side 
arcades and removing field blocking the path of the plasmoid. This run-away 
process results in the launch of the CME. In this scenario, the CME is triggered 
by internal reconnection, with breakout reconnection occurring later in the event. 
Therefore it is expected that signatures of internal reconnection e.g. fiaring of 
the central arcade, a SXR source at the reconnection site etc. will occur before 
those of the side lobe restructuring. 
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Figure 1.22: Schematic diagram describing internal tether cutting, from Moore 



& Sterling (2006). The filament is shown as the innermost loop. It is tethered by 



a central arcade and has neighbouring sidelobe arcades on each side. The entire 
system is then covered by overlying field lines. 
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Figure 1.23: Schematic diagram describing external tetlier cutting, from Moore 



& Sterling (2006), as in Figure 1.22 
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External tether-cutting, or tlie "breakout" model by Antiochos (1998), begins 
with a quadrupolar magnetic topology: a central arcade between two side lobe 



arcades (Figure 1.23). The structure inside the central arcade is such that no 
current sheet exists between the arcade legs (or if it does it is too thick to support 
reconnection) but a current sheet does exist between the top of the arcade and 



the overlying field, as shown in the top panel of Figure |1.23[ This can result from 
e.g. further emergence of the central arcade, which works to compress the null 
point between the central arcade and the overlying field without the generation 
of a current sheet between the arcade legs. Reconnection above the arcade shifts 
the force balance so that the central arcade begins to rise, stretching the field and 
drawing the legs of the central arcade together to create a second current sheet 



beneath the filament (Figure 1.23 middle panel). This results in run-away tether- 
cutting reconnection as in the internal tether-cutting case. Unlike the internal 
tether-cutting case, breakout reconnection begins first. Therefore evidence of 
heating or reconnection in the side lobes would be expected before evidence of 
the same in the central arcade. 



The third case we consider is the catastrophe model (e.g. Forbes & Isenberg 



1991 Forbes & Priest, 1995; Isenberg et ai, 1993). This model differs from 



the first two in that it is not triggered by magnetic reconnection. Instead, the 
continued shearing and twisting of the central arcade gradually evolves the field 
until it is forced out of force-free magnetostatic equilibrium. The field seeks a 
new equilibrium by erupting upwards, generating two current sheets, one between 
the stretched fields of the arcade legs and a second between the top of the arcade 



and the overlying field (middle panel Figure 1.24). It has been shown that this 



can occur without the use of magnetic reconnection (e.g. Chen & Shibata, 2000 
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Isenberg et al. 1993 ; Roussev et al. 2003 ). Following the formation of the current 



sheets, magnetic reconnection can take place and run-away tether-cutting drives 



the launch of the CME, as before (Figure 1.24 bottom panel). In this case, 
one would expect to see a rising flux rope before any indication that magnetic 
reconnection had taken place. 

While current theoretical models have made signiflcant progress from the ideas 
of Carmichael, Sturrock, Hirayama, Kopp and Pneuman there remains signiflcant 
work to be done in this fleld. The development of new and improved technology 
and instruments, e.g. the SWAP instrument on board Proba-2 with its extended 



fleld of view (see ^3.2.1 for further details) continuously pushes the envelope on 
the performance of current theories. 



1.4 Outline of thesis 

The work presented in this thesis attempts to improve the understanding of the 
connection between solar flares and coronal mass ejections. To date, the rela- 
tionship between these two phenomena has been fraught with complications and 
competing theories. This thesis attempts to better understand the behaviour of 
flares, beginning with the evolution of a conflned flare. Extending this study 
to investigate the evolution of an eruptive flare attempts to categorise both the 
similarities and differences between these two flare categories. 

Chapter |2] presents the atomic processes involved in producing the emissions 
observed from the Sun which are employed when discussing the methods of ob- 
servation in Chapter |3| As an extension to the discussion of instrumentation, a 
study of the calibration of EUV imaging telescopes is presented in Chapter |4j As 
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part of the Proba-2 team, the analysis of the SWAP instrument's sensitivity to 
temperature was investigated for coronal hole, quiet sun, active regions and flares. 
This was extended to investigate the corresponding responses of other EUV im- 
agers, namely TRACE, SOHO/EIT, STEREO A/EUVl, STEREO B/EUVI and 
SDO/AIA. 

Chapter |5] explores the evolution of a conflned flare. The physical mechanisms 
involved in the flaring process are investigated by comparing observations from a 
wide range of spacecraft to a 0-D hydrodynamic model, EBTEL. The combination 
of observations and theory allow for a more extensive investigation than otherwise 
possible. The investigation into the physics of flares is then expanded in Chapter 
[6] to include an "eruptive" flare. Using the knowledge and techniques developed 
in Chapter |5j the relationship between flares and CMEs are investigated. This is 
done through a study of the hydrodynamic evolution of the flare and the kinematic 
evolution of both the flare and the CME. 

Chapter [7] discusses the implications of the work presented in this thesis both 
from observational and theoretical perspectives and suggests improvements and 
developments for future work. 
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Chapter 2 
Atomic physics 



The basis of all solar observations, from radio to 'y-ray emission, requires a de- 
tailed understanding of how that particular emission is generated. This thesis 
focusses on the EUV and X-ray regimes. To better understand the production 
of EUV and X-ray line and continuum emission, the use of an atomic physics 
package, CHIANTI, is used. In this chapter, we describe how CHIANTI models 
emission line and continuum radiation in the range of interest, the assumptions 
made in the calculations and the capabilities of the package. 
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2.1 Introduction 

The outer atmosphere of the Sun is a hot, tenuous plasma. Light elements such as 
hydrogen and helium are completely ionised while heavier elements are at least 
partially ionised, depending on the temperature. The main processes involved 



in the ionisation of atoms, shown in Figure |2.1[ are collisional ionisation and 
excitation autoionisation while atoms recombine by radiative recombination and 
dielectronic recombination. In equilibrium, the ionisation fraction is the number 
density of a particular ion relative to the number density of a particular element, 
determined from the balance between the ionisation and recombination processes 



of that particular element. The ionisation fraction for Fe is shown in Figure 2.2 



calculated from the ionisation fractions of Mazzotta et al. (1998). 

The solar corona emits strongly in the EUV and X-ray part of the spectrum 
through both emission lines and continuum. In this regime, the spectrum contains 
a multitude of emission lines, two-photon, free-free and free-bound continuum. 
To fully understand the origin and the conditions of the plasma emitting such 
radiation, it is necessary to model the spectrum in this regime. This is not a 
trivial task. In order to do so, parameters such as the energy levels, transitions, 
radiative transfer probabilities and excitation rates must be well understood for 
each individual line. In this thesis, we make use of the CHIANTl]^ atomic physics 



package (Dere et al., 1997). The CHIANTI database is used extensively by the 
astrophysical and solar communities to analyse emission line spectra from astro- 
physical sources. It was established in 1996 has, over the years, been maintained 



^CHIANTI is a collaborative project involving NRL (USA), RAL (UK), and the following 
Universities: College London (UK), of Cambridge (UK), George Mason (USA), and of Florence 
(Italy). 
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Figure 2.1: Schematic diagram showing the main ionisation and recombination 



processes in the solar atmosphere. Adapted from Aschwanden (2004). 



by Ken Dere, Helen Mason, Brunella Monsignori-Fossi, Enrico Landi, Massimo 
Landini, Peter Young, Giulio Del Zanna. It contains the parameters listed above 
for the known emission lines and includes sample differential emission measure 
functions for use in the simulation of astrophysical spectra. 



2.2 Emission lines 



CHIANTI calculates the flux of an emission line by representing the flux in terms 



of a series measurable and theoretical parameters. Following Mariska ( 1993 ) and 



Dere et al. ( 1997), an electron in an excited state can spontaneously decay through 



a bound-bound process with probability Aj^i. In a species X of ionisation state 
m, an electron that transitions between an upper state j and a lower state i will 
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Figure 2.2: lonisation fraction of iron calculated using CHIANTI and data from 



Mazzotta et al. (1998). 



produce a photon of energy = AE = Ej — Ei = hu = hc/X (Figure 2.3). 



(2.1) 



The volume emissivity, (ej^i) of a plasma with upper level population, Uj, and 
lower level population density, rii, is given by: 



3 



ej^i = huj^iAjiTij [ergs cm s 



(2.2) 
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Figure 2.3: Two level diagram showing an electron decaying from an upper level 
with energy Ej to a lower level of energy Ei and emitting a photon of energy AE 
in the process. 

where Aj^i is the Einstein coefficient for spontaneous radiative emission giving the 
probability per unit time that the electron in the excited state will spontaneously 
decay to the lower state. For a volume of optically thin plasma, AV, which is a 
good approximation for the outer layers of the solar atmosphere, the flux observed 
at Earth at a distance R from the Sun is proportional to the number of emitting 
ions in the line of sight to the observer and the fraction of those ions that are in 
a given energy state producing the emission line. We can begin by writing the 
flux as: 



Practically, the volume element dV is deflned by the spatial resolution of the 
instrument. Assuming the lines observed are optically thin, emission from all 
material along R will be accounted for. 

The number density of the upper level can be determined from: 




[ergs cm ^ s ^] 



(2.3) 



(2.4) 



n. 



'ion ^el 
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where nj/rtion is the population of the upper level relative to the total number 
density of the ion and is a function of temperature and density, riion/nei is the 
relative abundance of the ion and is a function of temperature, ngxjnu = ^ei is the 
element abundance relative to hydrogen and nn/ue is the hydrogen abundance 
relative to the electron number density. We can therefore rewrite flux of a line at 
Earth as: 

The number density of the excited state, rij must be populated by balancing the 
excitation processes with de-excitation processes. The energy state of an emit- 
ting ion can be changed by a range of different process. For example, the energy 
of the ion will change when an electron is excited into a higher level by e.g. a 
collision. The energy will again change if that electron decays back to its original 
state. These processes are generally faster than the processes that effect the ioni- 
sation state, as shown in Table 2.1| ( [Mariska 1993 p 18). This makes separating 
excitation/de-excitation and ionisation/recombination calculations possible. The 



excitation and de-excitation processes are shown in Figure 2.4 Each diagram in 



this figure corresponds to a term in Equation 2.6 An electron can be excited into 
a higher energy state by one of four processes: a collision with a free electron, a 
collision with a free proton or simulation by radiation. The electron can decay by 
one of four processes: electron and proton collision, stimulation by a photon or 
by spontaneous radiative decay. This is represented mathematically in Equation 



2.6 and schematically in Figure 2.4 
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Figure 2.4: Schematic diagrams of the excitation and decay processes with each 
figure corresponding to a term in Equation 



2.6 



Adapted from ( Aschwanden , 2004 ) . 



e~ collision p"*" collision stimulated spontaneous 



(2.6) 



The C^^ terms, given in cm^ s~^, are the collision rate coefficients for electrons 
and protons (where a = e and p respectively) for the 6*^ - c*^ transition. Ri, c 
in refers to the stimulated absorption coefficient for the h^^ - c*^ transition. 
The radiative transfer probabilities, A^^c for the 6*^ - c*'^ transition is measured 
in . These parameters act to populate and depopulate the excited level j. 
Initially, CHIANTI did not consider the proton excitation rates. However, they 
were later introduced to take account of the fine structure transitions in highly 



ionised systems (Young et al, 2003). 



63 



2.2 Emission lines 



Table 2.1: Transition rates of important atomic processes calculated for Civ at 
a temperature of 10^ K and a density of 10^*^ cm~^ (from 



Mariska 



1993). 



Process 


Rate 


Characteristic time 








(De-)excitation processes 


Collisional excitation 




2 X 10^2 


Collisional de-excitation 




2 X 10-3 


Spontaneous radiative decay 




4 X 10-9 


ionisation / recombination processes 


Collisional ionisation 




107 


Autoionisation 


^e'^ionQauto 




Total ionisation rate 




107 


Radiative recombination 


^e^ion^rad 


88 


Dielectric recombination 






Total recombination rate 
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In what is known as the coronal approximation, it is assumed that the pop- 
ulation of excited states occurs primarily by collisional excitation by electrons 
from the ground state and the de-population of excited states occurs primarily 
by spontaneous emission. Since the majority of electrons are in the lower (ground) 
state, we can say ^ ■ Ui = rii 



. Thus, Equation 2.6 can be approximated as: 



k 



(2.7) 



With rij/riion = UeCfj/ Aj^k, the emissivity can now be written as: 



{21 



The electron collision rate coefficient Cij for Maxwellian distributed electron 
velocities is given by: 



(^i,j{v)f{v)vdv 



(2.9) 
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where 



fiv) 



An 



m 



27ikBT 



3/2 



2 ; l/2mv'^ 



(2.10) 



(jjj is the electron excitation cross section by collisions and is commonly ex- 
pressed in terms of the collision strength Qij{E), the incident electron energy 
E = l/2mt>^ (measured in Rydbergs), the Bohr radius and the statistical 
weight of the i*^ level, Ui: 



a. 



UiE 



Ui is introduced to adhere to the principle of detailed balance, ensuring the net 
exchange between any two levels will be balanced i.e. the number of excitations 
caused by electrons in range dEi is balanced by collisional de-excitations by 
electrons in range dE2 so Ei = E2 + ^Eij. This is only valid for energy levels 
in thermal equilibrium and is not valid for the coronal approximation. However, 



CHIANTI does not make the assumption of Equation 2/7 and so the principle of 
detailed balance is used for the calculation of collisional excitation rates as: 



Cj^i — CijCxp 



00 



AEi 



(2.12) 



Combining Equations 2.9, 2.10 and 2.11, the collision rate coefficient can be 
written as: 

(2.13) 



8.63 xlO-'^ P ^ f~AE\ ,^ 



UiksT^/^ J^E. 



Equation 2.5 can now be written as 



F{X) 



hu 8.63 X 10 ^fij uh riei 



2 



ne riH Jav VT Uel 



-exp 



-AE 
kT 



dV 



(2.14) 
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In order to solve this equation for a given emission line, CHIANTI acquires 
the parameters from a number of sources. Energy level information is obtained, 
where possible, from the National Institute of Standards and Technology (NIST) 



database of observed energy levels (Martin et ai, 1995). This has been supple 



mented by theoretical estimates where the energy levels are not known (Dere 



et a/., 1997). These have been calculated using the UCL SSTRUCT program 



(Eissner et a/., 1974). Einstein coefficients for radiative transitions are, for the 



most part, obtained from literature or calculated using the SSTRUCT code. Col- 



lision rate coefficients (Cij) are scaled according to Burgess & TuUy (1992) and 
the de-excitation rates are obtained from the principle of detailed balance. The 
elemental abundance and ionisation fraction of the plasma are user defined. Coro- 



nal abundances and the ionisation fraction of Mazzotta et al. (1998) are used 
throughout this thesis. 



2.2.1 Contribution functions and emission measures 



The temperature sensitive components of Equation 2.14 can be extracted in what 



is known as the contribution function, G{n, T), samples of which can be found in 



Figure 2.5 The contribution function is given by: 



riei 



kT 



(2.15) 



and defining the emission measure (EM) to be the amount of emitting plasma in 
a given volume dV 



EM = I nidV 



(2.16) 
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or column depth dh 

EM = J nldh (2.17) 
the line flux can then be written as: 

F{X) = c X EM X G{T) (2.18) 

where c takes into account the physical constants, the hydrogen abundance [nn / Ue 
and the elemental abundance {nei/nn)- The contribution function provides in- 
formation regarding the formation temperature of a given line. In this thesis, we 
utilise the G{T) of emission lines to convert between the intensity of a line and 
the EM of the line at a given time and temperature (Chapter [5, Raftery et al. 



2009). This is done by inversion of Equation 2.18 Estimation of density can be 



made from: 



1^ 

This inversion requires the estimation of V , the volume of the emitting source. 
This of course can be a signiflcant problem. Even if the source can be imaged 
using e.g. a spectrometer raster, the two dimensional nature of the observation 
immediately places uncertainty on the volume. However, with little or no other 
options, this is sometimes the only method by which the density can be calculated. 

The emission measure of a plasma is suitable, assuming that the spectral line 
is emitted over a homogeneous, isothermal volume. For cases of multi-thermal 
plasma it is useful to deflne the differential emission measure (DEM) for a volume 
as 

DEMiT) = nl [dVjdT) (2.20) 
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Figure 2.5: Contribution functions G{n,T), obtained from CHIANTI for Hei 
(584.33 A), Ov (629.73 A), Mgx (524.94 A), Fexvi (360.75 A) and Fexix 
(592.23 A). These were calculated using a density of rie = 5 x 109 cm ^, coro- 



nal elemental abundances and ionisation fractions from Mazzotta et al. (1998) 



or for a column of depth as 

DEM{T) = n\ {dT/dhy^ (2.21) 

The DEM relates the amount of material in the temperature interval dT in a 
volume and so can give information about the structure of the atmosphere. For 
example, the temperature gradient dT/dh in the column DEM can be used to de- 
termine the temperature gradient between the transition region and the corona. 
It also takes account of the large volume that is the chromosphere and the tem- 



peratures found within it, as shown in Figure 2.6 Likewise, the thinness of the 



transition region results in a reduced DEM over transition region temperatures. 
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Figure 2.6: Quiet Sun DEM function of the solar atmosphere calculated using 
values from CHIANTI. 

The large volume of the corona means the DEM is increased from transition re- 
gion values, although the low densities found in the corona suppress the DEM 
somewhat. During a solar flare however, the coronal temperatures and densities 
are raised considerably resulting in a DEM peak at high temperatures. CHI- 
ANTI supphes sample differential emission measure files. These have been calcu- 
lated from a combination of the work by Pottasch (1964) and UV/EUV spectra. 
CHIANTI does not currently include volume emission measures and so column 
emission measure must be considered throughout while using the database. 
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Figure 2.7: Continuum emission from free-free, free-bound and two-photon pro- 
cesses in the wavelength range 1 - 500 A calculated using CHIANTI. 

2.3 Continuum emission 



There are three primary mechanisms by which continuum emission is formed: 
free-free, free-bound and two-photon. The relative intensity of these processes 



between 1 and 500 A are shown in Figure 2/7 In the SXR range, free-bound 
continuum has a significant contribution, however it is clear that free-free contin- 
uum dominates throughout the spectrum. Therefore, we shall focus on free-free 
continuum emission here. 



2.3.1 Free-free continuum 

In a hot coronal plasma, free electrons and ions can suffer multiple interactions. 
The most frequent of these is when a free electron is scattered in the Coulomb 
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Figure 2.8: Schematic diagram of how free-free Bremsstrahlung continuum emis- 
sion is formed: as the path of an electron is bent by interaction with the electric 



field of a stationary particle, it emits a photon. (Aschwanden 2004). 



field of an ion (Z). The scattering is such that the electron remains free after the 
interaction. This is what gives rise to the name "free-free" continuum, otherwise 
known as Bremsstrahlung or braking radiation. During the scattering process, 
the electron loses some of its energy which is released as a photon with energy 
= Ep. — Ef, Thick target Bremsstrahlung occurs when electrons are 

/ '^i7icia ^scat O O 

accelerated to high energies in a collisionless plasma and become collisionally 
stopped when they interact with a thermal plasma such as the chromosphere. 



Following Aschwanden (2004), the Bremsstrahlung total power Pj(f, z/) radi- 



ated by a single electron in a plasma of rii ions in a frequency range du is 



Pi(v, v)dv = nivar{v, v)dv 



(2.22) 



where cxj. is the radiating cross section in cm? ergs Hz ^ and is given by 



, , 16 Z'^e^ /■''™"^ dh Ti , ^ 



(2.23) 
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where 6, the distance from the electron to the atom is shown schematically in 



Figure 2.8 and g(i',T) is the Gaunt factor which is approximately unity in the 



corona. The volume emissivity for electrons can be written as 



^^ = ^ / P^iv,u)f{v)dv (2.24) 



where the Maxwell-Boltzmann distribution f{v) is given by 



f{v) = Att ( ] v^exp ( 4^ ) (2.25) 

■'^ ' ^ 27ikBT J ^ V 2A;bT ' ^ ' 



and the flux of emission at Earth can be written as 



F oc y e^dV (2.26) 
In the corona, g{v,T) ^ 1, n = rii = Ug, and since we are predominantly dealing 



with protons, Z ^ I. Therefore, we can combine Equations 2.22 to 2.26 to write 



the Bremsstrahlung flux in terms of the photon energy eph = hv 

F(e,.) oc exp j T-'/'ri'dV (2.27) 



Remembering that emission measure EM = J n^dV we can write Equation 



2.27 



F(e,,) oc exp (^-^^ T-'/'EM (2.28) 

This function demonstrates the sensitivity of the Bremsstrahlung spectrum to 
both temperature and emission measure. These are the two significant factors 
taken into account when fitting the thermal component of the RHESSI and SAX 
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Figure 2.9: Photon spectrum of a large solar flare extending from SXR out 
to 7-rays. The spectrum shows the different processes that dominate the dif- 
ferent regimes. The SXR spectrum is dominated by thermal electrons, HXR 
(Bremsstrahlung) radiation by non-thermal electrons and 7-rays by energetic elec- 
trons, ions and pions. The electron-positron annihilation line is shown at 511 keV 



along with the neutron capture line at 2.2 MeV. (Aschwanden, 2004) . 



spectra in Chapters [5] and |6} Figure 2.9 shows an example of a solar energy 
spectrum where the thermal component is dominated by free-free continuum. 
The amplitude of the thermal spectrum will change with emission measure while 
the width of the spectrum will change with temperature. 

The CHIANTI code is extremely useful in identifying the emission lines and 
continuum found in the solar spectrum in the EUV and X-ray range. A com- 
parison of theoretical line intensities to observed line intensities can reveal the 
physical conditions required for their production, information that would oth- 
erwise be unattainable. This can be especially helpful when determining the 
contributions from blended lines observed by an instrument that cannot resolve 
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them, such as those in Chapter [5j The database has been constructed in such 
a way as to facihtate the easy distribution and updating of files. It is freely 
accessible through the internelj^ and through a package in the Interactive Data 



Language and SolarSoftware package (SSWIDL; Freeland & Handy, 1998). 



^ http : / / www . chiant i. rl .ac.uk/ 
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Chapter 3 



Instrumentation 



This chapter presents details of the instruments used in this thesis and how the 
observations are obtained, prepared and analysed. Beginning with CDS, a discus- 
sion of how data are obtained and analysed is presented. Following this, EUV 
imaging telescopes are discussed, in particular SWAP, TRACE and EUVI. Ex- 
tending our interest into the heliosphere, the techniques employed to analyse the 
data from STEREO Corl and Cor2 coronagraphs are presented. Going from the 
EUV and white light emission into X-rays, we begin with the GOES Satellites, 
focussing on the interpretation of observations for diagnostic purposes. Finally, 
the capabilities of RHESSI and the Mercury MESSENGER Solar array for X- 
rays are laid out. The imaging techniques used for RHESSI and the spectroscopic 
capabilities of both instruments are reviewed, focussing in particular on their per- 
formance in the thermal energy range. 
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3.1 The Coronal Diagnostic Spectrometer 

The Solar and Heliospheric Observatory (SOHO) was launched in 1995 as a joint 
endeavour between ESA and NASA and has been providing uninterrupted obser- 
vations of the Sun for 13 years (with the exception of a 3 month period, during 
which communication with the space craft was lost). Orbiting the first Lagrange 
point - LI - at 1.5 million km from the Earth, it is always on the day-side of the 
Earth. This thesis takes advantage of observations from one SOHO instrument 



in particular - the Coronal Diagnostic Spectrometer (CDS; Harrison et al, 1995). 
CDS is an EUV spectrometer designed to study the solar corona. Initially, the 
instrument was not used for studying solar flares for fear of damaging the detec- 
tors. It has however become an integral part of flare studies. This thesis utilises 
five emission lines at a wide range of temperatures, from He I (~0.03 MK) to 
Fexix at ~8 MK. These data can provide users with temperature and emission 
measure diagnostics, along with Doppler velocities. 

CDS is a Wolter II grazing incidence telescope whose optical layout is shown 



in Figure pTi] Light passing through a scan mirror is focussed onto one of six slits 
(2" X 2", 4" X 4", 8" X 50", 2" x 240", 4" x 240", 90" x 240"). The beam is passed into 
one of two telescope apertures by reflection at grazing incidence. Light entering 
the Grazing Incidence Spectrometer (CIS) hits a grating in grazing incidence and 
the spectrum is dispersed onto one of four 1-D detectors. Light entering the 
the Normal Incidence Spectrometer (NIS) passes through two toroidal gratings 
in normal incidence. A small out-of-plane tilt in each of the diffraction gratings 
(mounted side by side) means that the spectra formed by both gratings can 



be dispersed onto a single 2-D CCD detector as shown in Figure 3.2 (note in 
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Figure 3.1: CDS optical layout (Foley, 2003) 



this figure the spectra are displaced. In practice the spectra are dispersed in 
the central regions of the CCD). Only the central portion of the CCD is used 
for observations. In each of the four corners, a 10 pixel xlO pixel box remains 
unexposed. These regions are used for debiasing and will be explained later (see 



^3.1.1.1). NIS is a rastering spectrometer that builds images by stepping a slit 



from right to left using a scan mirror. A schematic diagram of such an image is 



shown in Figure 3.3 Here, each column represents the row of pixels along the 
slit from which a spectrum can be obtained. As the slit moves from right to 
left, a 3-D structure is built in solar x, Y and A. This enables the observer to 
extract the spectrum for an individual, or group of pixels. This type of image is 
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Figure 3.2: A quiet Sun spectrum as recorded by CDS by passing light through 
two gratings (NISI above and NIS2 below) and imaging the dispersed spectrum 
onto a detector. 
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Figure 3.3: Schematic diagram of an CDS raster. The slit is located along solar 
Y and stepped from right to left along solar X. At each AXAY, the spectrum is 



recorded onto a CCD (Haugen, 1996) 
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achieved at the expense of temporal resolution. For CDS, the effective cadence 
of a single image (i.e. the time it takes to step through the x dimension) is 
~11 minutes. Alternatively, if timing is an important aspect of the investigation, 
the slit position can remain fixed and the x dimension can be exchanged for 
temporal resolution. The instrument has two gratings with wavelength ranges 



of 310 - 380 A for NISI and 520 - 630 A for N1S2, shown in Figure |3j with 
the brightest lines identified. Due to telemetry limitations, only a portion of 
the whole spectrum can be down-linked. The selection of lines to study, the 
slit used to generate the spectrum and the exposure times are three of the main 
components that make up a CDS study. 

In this thesis, we utilise the CDS observing study flare_ar. This focusses 
on five emission lines spanning a broad range of temperatures. The He I line at 
548.5 Ais an optically thick chromospheric line and has a formation temperature 
of 3 X 10^ K. A 0.25 MK O v line is observed at 629.8 A. Coronal lines Mg x (625 A) 
and Fexvi (350.89 A) are formed at 1.2 and 2.5 MK respectively. Finally, an 
Fexix flare line that has a formation temperature of 8 MK is observed at 592.3 A. 
Each raster, taken with the 4" NIS slit utilising only 180" in the Y direction, 
consists of 45 slit positions, each 15 s long. This give an effective cadence of 
~11 minutes and a 180"xl80" field of view. 

3.1.1 CDS data analysis 

Each CDS file contains spectral and spatial data of a single raster, along with 
engineering information. This information is stored in a FITS file (Flexible Im- 



age Transport System; Wells et ai, 1981). FITS is designed specifically for sci 
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Figure 3.4: A "Sun as a star" spectrum observed with NISI and NIS2 across the 
full spectral range of both gratings. The brightest lines are identified within both 
spectral ranges. 

entific data and is widely used for storing photometric and spatial calibration 
information, together with image origin metadata. In order to retrieve spatial 
and spectral data from the fits structure, routines from SSW are used. SSWIDL 
consists of a "tree" of routines developed by the solar physics community, with 
"branches" for the preparation and analysis of different instruments. 

3.1.1.1 Data prepciration 

The CDS branch of the SSW tree consists of a series of routines to retrieve and 
prepare data. The following steps were applied to many of the other instruments 
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Figure 3.5: The "salt and pepper" effect caused by cosmic rays on a SOHO/EIT 
image. 

used for this thesis, especially the EUV imagers. The preparation steps included: 

• Cosmic ray cleaning: NIS, and all space-based CCD detectors are subject 
to contamination by cosmic rays passing through the instrument. When this 
happens, the cosmic ray can leave a high intensity streak across images. 



Figure |3.5| shows the "salt and pepper" damage these high energy particles 
can have on an image. Using cds_CLEAN, any pixels above the Sex level are 
replaced using bilinear interpolation of its neighbours. 

Flat fielding: Each CCD pixel has a different response to incoming pho- 
tons. The difference in gain and dark current can be removed by dividing 
data by pre-fiight fiat-field exposures. This should result in a uniform sen- 
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sitivity across the pixel. For NIS, the flat flelding is done automatically 
using VDS-CALIB when NIS.CALIB is called. 

Debiasing: A DC bias potential is applied to the CCD before it is read 
out. This ensures that the readout is always positive. This bias is removed 
using VDS_DEBIAS (again called automatically by nis_CALIb). In each of 
the four corners of the CCD, a small (10 pixxlO pix) region remains unex- 
posed. Subtracting the average value of each of these dark "frames" from 
the respective quadrant removes the DC bias. 

Calibration: The conversion of raw data numbers into physical units re- 
quires consideration of the reflectivity of the primary, secondary and scan 
mirrors, the efficiency of the grating, the quantum efficiency of the detec- 
tor and the wavelength calibration. For CDS, these effects are considered 
in NIS_CALIB along with the flat-flelding and debias techniques mentioned 
above. This technique, applied to EUV imagers, is described in detail in 
Chapter |4} 

Slant and tilt correction: The CDS wavelength axis is not perfectly 
aligned with the detector axis. This is as a result in an offset between the 
gratings and detector. This results in the dispersed spectrum appearing 
slanted relative to the detector (this effect can be seen clearly in Figure 



3.2). This tilt, along with an error along the Y axis (due to an error in the 



slit assembly) is removed using the nis_ROTATE routine. 
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3.1.2 CDS line fitting 

Once the data has been prepared and cleaned, an investigation of individual lines 
can take place. As mentioned before, spectral data is available for individual pix- 
els. However, in many cases (especially for non-fiaring conditions), the intensity 
of a line in a single pixel is too low to be accurately modelled. Therefore, it is 
often necessary to sum the emission from multiple pixels in order to achieve a 



well defined spectral line, as shown in Figure 3.6 



Before the three month loss of contact, the CDS emission lines were modelled 
using Gaussian distributions. However, damage to the instrument resulted line 



profiles with extended wings (see comparison in Figure 3.7). In order to anal- 



yse these lines, a single Gaussian was insufficient. Thompson (1999) modelled 



post recovery line profiles, approximated by one or more broadened Gaussian 
functions: 



G{\) = exp 



W{X) 



A-Ar 



(3.1) 



(3.2) 



(^) +1 

where G represents the Gaussian term, W represents the wings, A is the wave- 
length, Aq is the rest wavelength, a is the Gaussian width and a' is the FWHM. 
The combined profile can be expressed as 



5(A) = / [(I -a)G(A) + aVr(A)] 



(3.3) 



where / is the amplitude of the line profile and a is the amplitude of the red or 
blue wing. The best fit of a particular function to a line was established by 
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Figure 3.6: CDS Fexvi image (top left) and spectra for a single flaring pixel 
(spectrum A), single quiet Sun pixel (spectrum B) and multiple quiet Sun pixels 
(spectrum C). 



minimisation. 



X 



E 



(3.4) 



where af is the RMS error in lobs- 



If a spectral window is known to contain more than one line, i.e. the line of 



84 



3.1 The Coronal Diagnostic Spectrometer 



(7) 
C 
0) 



0) 

> 



0) 



0,10 



0.0 




360 



055 370 
Wavelength (A) 

NIS-2 



375 



1 00.00 



0,00 



.00 



0.10 



0.0" 




575 



580 5a5 

Wavelength {A^ 



590 



595 



Figure 3.7: Comparison of line profiles from before and after loss-of-contact for 
NISI and NIS2. The faint lines illustrate the lines profiles prior to loss-of-contact 
on 25 June 1998. The bolder lines show the broadened profiles after recovery in 



October 1998 (Thompson, 1999). 
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Figure 3.8: An example of a blended Fexvi line, where red corresponds to the 
Fexvi line, the green line is a blend, the yellow histogram is the total fit and the 
yellow horizontal line is the background level. The black histogram corresponds to 
observed data and the vertical lines are the centroid position. 

interest is blended or contaminated, a second, and sometimes third broadened 
Gaussian function is required to account for and remove all line blends (e.g. 



Figure 3.8). The relative intensities of the blended lines changes throughout 
the course of a flare. As the temperature rises, "flare" lines, i.e. those formed 
at high temperatures, become more dominant. Often it is the case that at the 
peak of a flare, the lines of interest completely dominate the spectral window. 
However, as the flare cools, the blends become more significant and can influence 
the parameters such as width and centroid of the line of interest. 



3.1.3 CDS Doppler shifts 

A useful diagnostic tool that emission line spectroscopy provides is the ability 
to analyse the velocity of the plasma. Moving plasma will emit Doppler shifted 
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radiation. Measuring the shift in wavelength can give the magnitude of the ve- 
locity and measuring the direction of the shift i.e. towards the blue or red end 
of the spectrum, will give the direction. Caution must be taken when analysing 
Doppler shifts as only the line of sight component of the velocity will register any 
shift. If the plasma motion is not parallel to the line of sight, then a correction 
factor is required. The Doppler shift is given by: 



V A — Ao 
c 



Ao 



(3.5) 



where v is the plasma velocity, c is the speed of light, Aq is the rest wavelength 



and A is the shifted wavelength. Figure [3^ gives a clear example of Doppler shifts 
for the footpoints of a loop. The deviation from the rest wavelength (solid line) is 
shown, along with the Gaussian profiles and the shifted wavelength (dotted line). 



CDS has been used to resolve velocities to an uncertainty of 4.7 km s ^ ( Gallagher 



et al.: 1999) for He I, although a more conservative uncertainty of 10 km s ^ is 



often used (e.g. Brekke et al. 1997; Raftery et al. , 2009). 



3.1.3.1 Rest wavelength calculation 

The calculation of the rest wavelength for individual emission lines is essential 
for the accurate measurement of the Doppler shifts. For the cooler lines (He I, 
Ov, Mgx Fexvi) the rest wavelength was obtained from a region of quiet Sun, 
since it contains plasma emitting at the temperatures at which the above lines 
are formed. The rest wavelengths for CDS chromospheric and transition regions 



Brekke et al., 


1997 


Doschek et al. 


1976 



Feldman et al. ^ 1982). Brekke et al. (1997) measured the redshift of the Ov 629 A 
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Figure 3.9: Fexix map with footpoints highlighted by black boxes. The inte- 
grated spectrum of each footpoint are shown with the centroid and rest wavelengths 
marked by dotted and solid lines respectively. The red component shows the fit of 
the FeXix lines, the green and blue components are blended lines' fits, the yellow 
horizontal line is the background and the yellow histogram shows the overall fit. 



line using the Solar Ultraviolet Measurements of Emitted Radiation (SUMER; 



Wilhelm et al., 1995). They found a maximum shift of 16±3 km s ^ for the Ov 



line and 6±1.5 km s^^ for Mgx. Since no absolute Doppler shifts can be obtained 
with CDS, we adopt the upper limit of 16 km s~^ as the error on our observations. 
Obtaining a rest wavelength for the Fe xix line was more difficult since no Fe xix 
exists in the quiet sun. The orientation of the loop was such that any plasma 
flowing along the loop would have little or no line of sight velocity component. 
Therefore, it was assumed that the net velocity of plasma at the top of the loop 
would be zero as the components flowing upwards from each leg will cancel out, 
resulting in the "rest" wavelength for Fexix. 
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Figure 3.10: A schematic diagram of SWAP's optical components. 

3.2 EUV imagers 

3.2.1 Sun Watcher using Active Pixel System detector 
and Image Processing (SWAP) 

The launch of the ESA Proba-2 mission in November 2009 has brought with it the 
commissioning of the Sun Watcher using Active Pixel System detector and Image 



Processing (SWAP; Berghmans et al. , 2006 ): a full disk EUV imager with a single 



filter centred on the 171 A Feix/x line. With a plate scale of 3.1" per pixel, the 
spatial resolution of the instrument is comparable to SOHO/EIT. However, the 
extra wide field-of-view (54' x 54') makes this unique instrument ideal for tracking 
coronal mass ejections from the low corona out to more than 3i?0. 



swap's optical layout, shown in Figure 3.10 includes a set of aluminium 
foil filters. The front filter is supported by a polyimide membrane and the rear 
filter by a nickel grid. The filters, combined with EUV multilayer coatings on 
the mirrors tune the sensitivity of the instrument to between 171 A and ~183 A. 
With first light not expected until January 2010, we cannot include details on 
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image processing here, however the telescope was used in a calibration effort for 
EUV imagers described in Chapter |4] 

3.2.2 Transition Region and Coronal Explorer (TRACE) 



The Transition Region and Coronal Explorer (TRACE; Handy et al., 1999) is 
a space based imaging telescope, designed to take high cadence, high resolution 
images of the different layers of the solar atmosphere. The large aperture, small 
field of view (8.5' x 8.5') and small plate scale (0.5" pixels) make the resolution 
of TRACE the highest of any current EUV imagers, though at the cost of cov- 
erage area. Similar to SWAP, TRACE'S multilayer mirror coating and filters 
tune the sensitivity of the instrument to one of eight passbands in the UV and 
EUV. TRACE has seven filters, three of which are in the EUV range (171, 195 
and 284 A), four in the UV (1216, 1550, 1600, 1700 A) and a broad whitelight 
passband. TRACE data is available online and, like CDS, can be cleaned and pre- 
pared using SSWIDL routines. The basic techniques are similar to those described 



in ^3.1.1.1 



3.2.3 The Extreme Ultraviolet Imagers (EUVI) 



The Extreme Ultraviolet Imagers (EUVI; Howard et al. 2008) on board the twin 
Solar Terrestrial Relations Observatories (STEREO; Kaiser et al.. 2008), Ahead 
(A) and Behind (B) have provided the solar community with two full disk EUV 
imagers whose spatial resolution are almost 40% better than EIT. EUVI has 
four passbands; 304 A, 171 A, 195 A, and 284 A. Again, the multilayer coatings 
and filters tune the instrument sensitivity to these passbands. The cadence of 
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EUVI is significantly better than EIT (2 mins vs. 10 mins) and the 43'x43' field 
of view makes the EUVI instruments optimal for studying large scale coronal 



disturbances. EUVI (along with the coronagraphs listed in [3.3) are prepared 
and cleaned using the secchi_prep routine, in a similar manner to §3.1.1.1 



3.3 Coronagraphs 



The Cor coronagraphs from the Sun Earth Connection Coronal and Heliospheric 



Investigation (SECCHI; Howard et al. . 2008 ) suite of instruments on board STEREO 



B were used for work in this thesis. A coronagraph is an instrument that sim- 
ulates a solar eclipse by blocking radiation from the photosphere. This allows 
the faint emission from the extended corona to be integrated over the line of 
sight through the optically thin corona. The instrument is designed not only to 
eliminate photospheric emissions, but also to reduce stray light from aperture 
diffraction. 



3.3.1 Corl 



The Corl telescope (Howard et al. . 2008) is an internally-occulted coronagraph 



based on the classic design by Lyot (1939) that has a field of view from 1.4 - 



4 Rq. a schematic diagram of the telescope's optics are shown in Figure 3.11 



Light enters the telescope and is focussed by the objective lens onto the occulting 
spot. This spot absorbs light from the center of the field of view, removing much 
of the emission from the central source, in this case, the solar disk. This has the 
effect of forming diffraction rings around the occulting spot. A Lyot stop is used 
to remove these rings. Additional stray light is removed by placing baffles at 
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Figure 3.11: A schematic diagram of the Corl coronagraph. 



various points between the front aperture and the Lyot stop. A hnear polariser 
placed between the Lyot stop and the final lens is used to suppress scattered light 
(F corona emission) and to extract the polarised brightness signal from coronal 
emissions. 

3.3.2 Cor2 

Cor2 is an externally occulted coronagraph with a field of view of 2.5 to 15 Rq. 
Unlike Corl, the solar disk is occulted at the entrance aperture of the telescope. 
This has a vignetting effect on the image formed at the objective lens (this is why 
an internal occulter was chosen for Corl). The external occulter is cantilevered 
with a single pylon to reduce diffraction. Internal baffles between the occulter 
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and the objective lens help to reject stray light. Further to this, a heat rejection 
mirror guides any remaining photospheric light out of the light box, preventing 
the reflection off the back end of the occulter. An internal occulter rejects the 
bright fringe of diffracted light that forms around the primary occulter. The light 
path is then guided through a Lyot stop and a fllter wheel, onto the detector. 



See Howard et al. (2008) for further details 



3.4 Geostationary Operational Environmental Satel- 
lites (GOES) 

The GOES fleet consist of a series of satellites, the flrst of which was launched 
on 16-October-1975. This study makes use of data from the twelfth of the series. 
These spacecraft are used for "now-casting", providing data for severe storm 
evaluation, information on cloud cover, winds, ocean currents, fog distribution, 
storm circulation and snow melt on Earth. The instrument also takes continuous 
observations of the integrated soft X-ray flux of the Sun in the wavelength ranges 



0.5-4 A and 1-8 A (Figure 3.12). This data is frequently used for Space Weather 



forecasting. The classiflcation of a flare is obtained from the magnitude of soft 



X-ray flux observed by GOES (see ^1.2). This satellite is integral to the past 



present and future of solar physics, as it is the only instrument that has been 
studying the Sun for such an extensive period of time. 

Observations from GOES take the form of a flux measurement, taken every 



three seconds. Two channels are used: 0.5-4 A and 1-8 A. Following Thomas 



et al. (1985), the flux of a given channel i can be expressed as: 
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GOES Xray Flux (5 minute data 



Begin: 2007 Dec 29 0000 UTC 



-1-1 
-1-1 




Updated 2007 Dec 31 23:31:07 UTC 



Jan 1 

NOAA/SWPC Boulder, CO USA 



Figure 3.12: GOES 0.5-4 A and 1-8 A flux for the 29-Dec-2007 to l-Jan-2008. 
Note the extremely quiet conditions on the 29-Dec-2007 are suddenly disturbed by 
the rotation of a flare-producing active region onto the solar limb. Note the large 
C8.0 class flare at 00:45 UT on 31-Dec-2009 is the focus of Chapter |6| (Image 
taken from www.solarmonitor.org) 



Fi = EM / G{X)idX/G = EM x k 



(3.6) 



where 



6. = / G{X)dX/G 



(3.7) 



contains the temperature dependent components of flux - the transfer function, 



G(A), shown in Figure 3.13 and the wavelength range, dX - and is normalised by 
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Figure 3.13: GOES 0.5-4 A and 1-8 A transfer functions (White et al. 2005) 



the wavelength averaged transfer function, G. 

The ratio of two channels - 0.5-4 A (subscript 4) and 1-8 A (subscript 8) 
gives the temperature dependent value R{T), where 



R{T) 



(3. 



and the response of RiT) is shown in Figure 3.14 Once the effective temperature 



has been found from the ratio of the fluxes, fog can be determined and the emission 



measure can be calculated by inverting Eqn. 3.6 
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Figure 3.14: Flux ratio, R{T) response as a function of temperature (Thomas 



et al. 



1985) 



3.5 Reuven Ramaty High Energy Solar Spec- 
troscopic Imager (RHESSI) 



RHESSI (Lin et al.. 2002) is a NASA Small Explorer Mission designed to investi- 



gate particle acceleration and energy release in solar flares through imaging and 
spectroscopy of high energy continuum and lines. The instrument consists of nine 
bi-grid rotating modulation collimators (RMCs) in front of a spectrometer and 
nine cryogenically cooled germanium detectors, each one behind an RMC The 
sensitivity ranges from 3/6 keV up to 17 MeV depending on the attenuator state 



see Table 3.1), covering soft X-rays to 7-ray emission. 



The spacecraft, which rotates once every 4 seconds, is pointed and oriented 
by the Solar Aspect System (SAS) and the Roll Angle System (RAS). The SAS 
uses white light images of the solar disk to ensure spacecraft pointing is correct to 
within 1.5". This is done through accurate measurement of pitch and yaw angles. 
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Table 3.1: RHESSI attenuator states and the lowest corresponding energies that 
can be observed. 



Attenuator state Attenuator thickness Lowest reliable energy 
AO None 3 keV 

Al Thin 3 keV 

A2 Thick 6 keV 

A3 Thick + thin 6 keV 



The RAS uses background stars to position the spacecraft once every rotation, 
defining the spacecraft orientation in relation to the star field. 

RHESSI has dual analysis capabilities: imaging and spectroscopy. The details 
of each are very different and equally complex. They are detailed below. 



3.5.1 RHESSI imaging 

Imaging high energy photons is very challenging. This is because the wavelength 
of the light is shorter than the spacing between atoms used in detectors. As such, 
a more complex method of imaging must be employed to successfully capture 
the high energy emission from solar flares. This problem is overcome by using 
the fact that radiation creates a shadow when an object is obstructing the light 
path. This technique, collimator-based Fourier-transform imaging is detailed in 



Hurford et al. (2002). 



At a given instant in time, an off-axis source will either be occulted by the 
grid or not. As the spacecraft rotates, the grids will pass infront of the source, 
modulating the detected signal observed by the detectors as the grids cast a 
shadow on the detector. The modulation pattern will vary depending on location, 
intensity and size of the source. The detectors, which have no spatial resolution. 
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Figure 3.15: Schematic of one RHESSI RMC setup. (Hurford et al. 2002). 



record the arrival time and energy of the incident photons, recording the count 
flux function of rotation angle. 

The nine grids each have a different angular resolution, given by: 



9 = tan 



P 
2L 



(3.9) 



The pitch P is the distance between the edge of one slat and the corresponding 
edge of the neighbouring slat and L is the length of the collimator (see Figure 



3.15 for clarification). The grid pitches range from 34 fira to 2.75 mm in steps of 
•s/3. Thus, the grids provide resolution from 2.3" to 180". 

In order to reconstruct an image using RHESSI data, one must select an 
image reconstruction algorithm. Back projection is the most basic routine. A 
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"probability" (or dirty) map is created of the origin of each photon and the prob- 
ability map for all photons are then summed. In order to remove any artefacts, 
such as side lobes, a more complex algorithm should be used. The simplest and 
most efficient of these is the "Clean" algorithm, which has been employed for all 
RHESSI imaging in this thesis. Clean assumes a source can be represented by 



the superposition of multiple point sources. Figure |3.16| shows the various steps 
in the processing of a "Clean" image. Clean removes the sources in the dirty map 
{I{x,y)) which it believes are not real, leaving only the real sources {0{X,Y)) 
convolved with the instrument point spread function {P{x,y)). 

/ + 00 
0{X,Y)P{x-X,y~Y)dXdY (3.10) 
■oo 

The algorithm begins by first computing the dirty map /„ (or residual map for 



n 7^ 0) as with the back projection method (Figure 3.1Q\). The maximum value 



of the dirty map, is ascertained (Figure 3.16p ). The point spread function 



is applied to the point source (Figure 3.16p ) and removed from J„, leaving the 



residual map (Figure 3.16^ ). This process is repeated until the residual map 



is at the noise level or a negative maximum is reached. The final Clean map is 
the convolution of the maxima with the point spread function, added to the final 
residual map (Figure [3l6p ). 



3.5.2 RHESSI spectroscopy 

The second major function of the RHESSI instrument is for spectroscopy of high 
energy solar events. RHESSI has the capability to observe continuum from the 
soft X-ray regime to 7-rays, along with some strong line features (such as the 
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Figure 3.16: Steps involved in creating a CLEAN image, (a) shows a "dirty" 
or pre-processed map. (b) shows the residuals following the removal of a CLEAN 
component (c) which is the convolution of the maximum intensity pixel in the map 
(d) with the instrument point spread function. The resulting CLEANed map is 
shown in (e) along with scaled residual and component maps (f, g) and the scaled 
CLEAN components (h). 

Fe/Ni complex in the thermal range, the electron/proton annihilation line at 
511 keV and 7-ray lines in the 1-10 MeV range). 

Photons entering the Ge detectors lose their energy via photoelectric absorp- 
tion, Compton scattering and electron/hole pair production. The high voltage 
across the detector attracts the electron or hole to the relevant electrode, creating 
a pulse. However, not every photon behaves in exactly this manner. Some effects 
that may modify the spectrum are: 

• Absorption by the mylar blankets, cryostat windows and grids. 



Compton scattering out of detectors. 
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• Compton scattering off the Earth's atmosphere. 

• Noise in the electronics. 

• Detector degradation due to radiation damage. 

• Low energy cut off due to electronics. 

In order to accurately interpret the spectrum, it is necessary to understand how 
detectors interpret photons of different energies. This is known as the detector 
response matrix (DRM). The diagonal components represent the efficiency of 
the instrument for detecting photons at their proper energy and the off diagonal 
elements represent the changing of photon energy from their true energy to a 
different energy. 

The detected count spectrum can then be converted to a photon spectrum 
using 

CT = BG + DRM*PH (3.11) 

where CT is the count rate, BG is the background count rate and PH is the 
photon rate. 

The photon spectrum can then be modelled using the Object SPectral Exec- 



utive package (OSPEX; Schwartz et al. 2002) using a variety of photon spectra 



components. Alternatively, and more applicable to the studies carried out in this 



thesis. Equation 3.11 can be inverted and the model photon components can be 
converted to count flux and compared to the raw count flux observations. Since 
this is highly dependent on temperature and emission measure, the changes in 
these parameters can be observed more clearly in count space than in photon 
space 
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3.6 MESSENGER Solar Array for X-Rays 

The Mercury Surface, Space, Environment, Geochemistry and Ranging Instru- 



ment (MESSENGER; Santo et al, 2001) is a satelhte that is designed to study 



the planet Mercury. On board is the Solar array for X-rays (SAX; Schlemm et al. 



2007). SAX, a Silicon-PIN detector, is a component of a larger instrument de- 
signed to measure X-rays on the surface of Mercury. SAX is a calibration device 
that measures the incident SXR flux from the Sun. The high spectral resolution 
in the range 1-10 keV makes it ideal for measuring the temperature and emission 
measure of flares. 

The OSPEX routines in SSW have been adapted to accommodate SAX 
spectra. Therefore, the analysis tools used for RHESSI can also be applied to 
SAX spectra. As with RHESSI, the count spectrum is compared to photon models 



that have been convolved with the SAX DRM flle (Equation 3.11). 

The SAX instrument is sensitive to only thermal emission (i.e. it observes no 
hard X-rays). For isothermal plasma. Maxwell Boltzmann distributions are the 
only available functions with which to flt the spectra. For multi-thermal plasma, 
there are a number of options. All of these are based on differential emission 
measure functions. OSPEX has two built in functions: an exponential DEM flt 
and a power law DEM flt. Both flt a function (either exponential or power law) 
between the DEM peak at low temperatures (photosphere/chromosphere) and 
at high temperatures (corona). These models ignore the dip that occurs in the 
middle of the DEM curve at the transition region. This approach is reasonable 
for RHESSI observations in which high temperature emission dominates. How- 
ever, for SAX data, there is also a signiflcant contribution from low temperature 
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plasma. Therefore, the accurate measurement of the relative contributions of 
both the high and low temperature plasma is important. As such, the observed 
plasma is approximated by two isothermal functions: one corresponding to the 
low temperature component of the DEM and one to the high temperature com- 
ponent, thus modelling both the low and high temperature plasma observed by 
SAX simultaneously. 
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Chapter 4 

Temperature response of EUV 
imagers 

This chapter discusses the temperature response of six extreme ultraviolet imaging 
telescopes in the 171 A passband. The intent is to investigate the response of the 
instruments to various types of plasma, including coronal holes, quiet sun, active 
regions and flaring emission. Beginning with a differential emission measure dis- 
tribution for each type of emission, synthetic spectra at single temperatures were 
generated in the range 10^ to 10^ K. The product of these spectra with the instru- 
ments' wavelength response functions resulted in a series of isothermal transmitted 
spectra. The instruments' response to temperature was then obtained by integrat- 
ing the transmitted spectra over wavelength to give the expected counts per second 
per pixel at a single temperature and emission measure. Despite variations in 
the instruments' wavelength responses, the resulting temperature responses were 
found to agree very well. The most significant deviation was the offset between 
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temperature responses. This is a direct result of the efficiency and size of each 
instrument. However, comparing the responses for different types of emission re- 
vealed significant differences, especially for the coronal hole and fiaring scenarios. 
The peak temperature of the passband when observing coronal holes is 7.5 x 10^ K 
versus 1 x 10^ K for quiet sun and active regions and the response to fiaring 
plasma proved to be multi-peaked as a result of the significant contribution from 
free-free continuum at high temperatures. This method will enable scientists to 
directly compare the sensitivity of a range of instruments to different solar con- 
ditions. This work has been submitted to the Astronomy & Astrophysics journal 
(Raftery, Bloomfield, Gallagher, Seaton & Berghmans, 2010, submitted). Con- 
siderable assistance was given to the author by Dr. Shaun Bloomfield and Dr. 
Peter Gallagher. Access to BESSY data was facilitated by Dr. David Berghmans 
and knowledge of the Proba-2 spacecraft was imparted by both Dr. Berghmans 
and Dr. Daniel Seaton. 



105 



4.1 Introduction 



4.1 Introduction 



The upper layers of the atmosphere emit strongly in the EUV portion of the 
spectrum, from the chromosphere through to the corona. This part of the solar 
atmosphere contains features across many different size scales, from oscillations 



of coronal loops (on the order of arc seconds; De Moortel et al. . 2002 ) to EIT 



"waves" (on the order of R©; Thompson et al., 1998). Active regions are the 
strongest features visible in EUV emission on the solar disk. As the origin of 
solar flares and CMEs, the accurate interpretation of EUV emission from active 
regions is of particular interest to the author. The EUV part of the spectrum 
is dominated by emission lines formed across a wide range of temperatures from 
~10^ to more than 10^ K. With the strong dependence of temperature on height 
through the atmosphere, this suggests that EUV emission can originate from 
anywhere between the upper chromosphere and the corona. Density is also known 



to scale with height in the atmosphere (Figure 1.2). Since emission measure (EM) 



scales with density (Equation 2.16) it is clear that not only does the temperature 



of a feature give some indication of the height at which it is being observed, but 
we can also draw conclusions regarding the EM and density of the region, as 



Chapters 5 and 6 showed ( Raftery et al. 2009 2010b ) . 



Imaging of the Sun using full disk EUV images is a very useful method of 
studying solar features across many different size scales at a wide range of tem- 
peratures. Normal incidence telescopes make use of constructive interference to 
amplify the reflectivity of their mirrors. The use of varying layers of Molybde- 
num and Silicon tune the wavelength sensitivity to different portions of the EUV 
spectrum. At ~ 10 A wide, the wavelength sensitivity is broader than that of a 
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spectrometer but generally sufficiently narrow to focus on one or two emission 
lines under "typical" solar conditions. However, the solar spectrum can vary sig- 
nificantly when studying various features, even within the narrow band of sensi- 
tivity. For example, the emission coming from a coronal hole is significantly cooler 
then that coming from a solar flare, or even an active region. Coronal holes have 
low temperatures and densities (~ 10^ K, ~ 10'' cm~'^ respectively) and appear 
as dark regions in EUV images. A solar flare, on the other hand has emission 
at very high temperatures and high electron densities (~ 10'' K, ~ 10^^ cm^'^ 
respectively). The appearance of a cool line in a coronal hole spectrum or a hot 
line in a flare spectrum can have a remarkable effect on the temperature of the 
emission detected by the EUV instrument. This can lead to the misinterpretation 
of "bright" features in EUV images. 

In an attempt to avoid such confusion, investigations into the characterisation 
of EUV imagers' temperature responses have been carried out. A characterisation 
of each instrument has been carried out by their respective team and published in 
the relevant instrument paper (see ^4.2.1 for papers). Typically this is attained 
for a single type of solar feature (usually quiet-sun) using a constant emission 
measure (normally 10'^^ cm~^). This clearly does not account for the order of 
magnitude difference in EM between e.g. the transition region and the corona. 

With the level of understanding of atomic processes ever evolving, the re- 



sponses presented in instrument papers require continuous re-analysis. Del Zanna 



& Mason (2003), for example, highlighted the importance of including FeVlll in 



the calculation of the TRACE 171 and 195 A passbands. The Fevill line had not 



been included in the calculation by Handy et al. (1999) and had an signiflcant 
effect on the 195 A channel in particular. This, however was only carried out 
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for quiet-sun emission at a constant pressure. Phillips et al. (2005) continued 
the development of the TRACE temperature response by incorporating the high 
temperature contribution to the spectrum by free-free continuum, as predicted 



by Feldman et al. (1999). Although they expanded their study to include emis- 
sion from quiet-sun, active regions and flares, the response curves published by 



Phillips et al. (2005) have been normalised, making it difficult to compare to 



other curves. Brooks & Warren (2006) again improved the temperature response 
of the TRACE and EIT EUV filters using spectroscopic observations taken with 
CDS. Recalculating the ionisation fractions using the CDS data facilitated a more 
accurate calculation of the instruments' temperature response curves compared 
to those in the instrument papers. Again, however, this was only completed for 
quiet-sun emission. It is clear that the studies carried out in the past do not facil- 
itate the direct comparison of the temperature response of multiple instruments 
to many different types of plasma. 

In this chapter, version 5.2 of the CHIANTI atomic physics package is used 
to investigate the temperature response of six solar instruments in the commonly 
termed "171 A" wavelength band. These responses are calculated for four classes 
of solar plasmas: coronal holes, quiet sun, active regions and solar fiares. The 
variation of EM with height is taken into account through the use of a differential 
emission measure curve (DEM). The calculation of isothermal spectra is detailed 



in 5:4.2 The instruments under investigation and their corresponding wavelength 



responses are described in ^4.2.1 The nature of the transmitted spectra are 



presented in ^4.2.2 The resulting temperature response curves for the four solar 



conditions and for the six instruments are described in detail in ^4.3 Finally, the 



conclusion and implications of this study are discussed in ^4.4 
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Figure 4.1: The expected field of view of SWAP compared to AIA. Tlie image 
shown is an EIT composite of multiple wavelength bands. 

4.2 Method 



4.2.1 Instruments and wavelength responses 

There are six instruments under investigation in this chapter. EIT on board 
SOHO, TRACE, the twin STEREO/EUVI imagers, the soon to be launched At- 
mospheric Imaging Assembly on board the Solar Dynamics Observatory (SDO/AIA; 



Title, 2009) and the recently launched SWAP instrument on the Proba-2 ESA 
satellite. As a member of the SWAP instrument team, this telescope will be the 
focus of my investigation. 

SWAP is a single filter imaging telescope with the unique capability to off- 
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Instrument 


Plate scale 


Nominal cadence 


FOV 


Aperture area 




[arcsec/pixel] 


[Minutes] 


[arcmin] 


[cm2] 


SWAP 


3.1 


1 


542 


8.55 


TRACE 


0.5 


<1 


8.52 


706 


EIT 


2.6 


10 


452 


13 


EUVI 


1.6 


2 


432 


3.01 


AIA 


0.6 


<1 


412 


113.1 



Table 4.1: Specifications for tlie EUV instruments used in this investigation. 



point out to 3 Rq. It was launched successfully in November 2009 with first light 
expected on 26 January 2010. At 54' x 54', its field of view already surpasses 



that of all the other instruments (see Figure 4.1 and Table ?? for comparison to 
instruments). With its remarkable ability to move the field-of-view to exclude 
the solar disk, this technology demonstration will enable users to image the solar 
corona in EUV to greater distances then ever before using the first active pixel 
sensors (APS) flown on a satellite. The high cadence of the instrument, at 1 
image per minute, will allow investigations of fast moving phenomena such as 
CMEs and EIT "waves" to be carried out in great detail. This will be especially 
important in the coming years as the STEREO spacecraft move further from the 
Earth. As their angle increases, their cadence will be reduced to due telemetry 
limitations. In addition, their field of view will be close to 180° from each other, 
making co-ordinated studies with each other and Earth orbiting instruments ex- 
tremely difficult. SWAP, with its similar capabilities to EUVI, will be an excellent 
replacement for these instruments. 

Launched in 2006, the EUVI instruments on STEREO A and B are also full 
disk imagers, though with a slightly smaller field of view compared to SWAP. 
These twin spacecraft have image cadence of ~2 minutes and have extremely 
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high resolution with a plate scale of 1.6"/pixel. The instruments have four EUV 
fihers (171, 195, 284 and 304 A) and have provided users with the first EUV 
stereo images of the Sun. 

EIT, launched in 1995 is often used as a "third eye" for the EUVI instruments. 
With similar field of view and the same filters as EUVI, it can often provide a 
complementary or comparative view of an event. For example, in Chapter [6] 



(Raftery et a/., 2010b), the EUVI instrument was primarily used for diagnostic 
purposes as the flare was mostly occulted in the EIT field-of-view. However, the 
EIT lightcurve provided a very useful tool for proving that a "double peak" ob- 
served in the GOES lightcurve was as a result of occultation and not a secondary 
event. EIT has significantly reduced cadence and resolution compared to the 
newer instruments. However, as a pioneer in the field of solar EUV imaging, it 
was once revolutionary in its capabilities. 

TRACE is another instrument that transformed solar physics using its seven 
passbands, 3 of which are in the EUV range (171, 195 and 284 A). Unlike the other 
EUV instruments, TRACE does not have a full disk field-of-view. Instead, it has 
extremely high resolution. With a cadence of less than one minute and a plate 
scale of 0.5"/pixel, this instrument revealed small scale features that had never 
been seen before. TRACE also deviated in the type of detector used. Unlike 
the other instruments listed here, TRACE used a front-illuminated CCD. The 
structures in front of the photosensitive portion of the front-illuminated device 
act as a deadlayer, reducing its sensitivity. Therefore, the overall efficiency of the 
TRACE instrument is expected to be lower than other instruments, despite its 
large aperture. 

The AIA instrument has significant heritage from TRACE, although it has 



111 



4.2 Method 




Figure 4.2: Images of the Sun in all four instruments currently taking measure- 
ments (i.e. excluding AIA and SWAP). Top left shows EUVI A in 284 A. Bottom 
left shows EUVI B in 195 A. Top right shows EIT in 304 A and bottom right 
shows TRACE in 171 A. In the centre there are four images corresponding to the 
100" X 100" regions marked on the main images. 



reverted to the use of a full disk, back-illuminated structure. This instrument, due 
for launch in early 2010, will be another in a long line of revolutionary telescopes. 
With resolution comparable to TRACE, nine EUV passbands (94, 131, 171, 193, 
211, 304, 335 A), full disk resolution and a cadence of less than a minute, AIA will 
hopefully reveal even more new science of the solar corona. A sample TRACE 



image is shown in the bottom right of Figure 4.2 AIA will have similar resolution 
to this but for the entire disk. 



Figure 4.2 highlights the different resolutions of the instruments in question. 
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This figure shows images of a recent period of activity on 22-Nov-2009 from 
the four telescopes currently taking data. In the center is a small 100" x 100" 
region to more clearly show the resolution differences. In addition, the different 
wavelength bands in which these images are taken emphasises the differences in 
solar emission at different temperatures. At 2 MK, the 284 A image is the hottest 
image and it is clear that the structures are significantly more diffuse that the 
other images. The 195 A image is showing more structure as is the somewhat 
cooler, with a peak temperature of just above 1 MK. Slightly cooler again, the 
171 A image displays a similar level of structure to 195 A, although the details 
are far superior thanks to TRACE'S high resolution. Finally, the EIT 304 A at 
chromospheric temperatures shows a very different Sun. With little extended 
emission and showing prominences on the limb, it is clear that the resolution of 
EIT is quite poor. The structure one would expect to observe in the chromosphere 
is completely blurred out. 

The wavelength responses R{X) of five of the imagers (all but SWAP) were 
obtained from SSW routines or from the instrument teams. The SWAP response 
was calculated from observations taken at the Berliner Elektronenspeicherring- 
Gesellschaft fr Synchrotronstrahlung (BESSY; Berghmans et al. private commu- 
nication). This experiment consisted of passing a calibrated synchrotron beam 
through the telescope's optics at wavelengths between 165 A and 195 A. Compar- 
ing the transmitted signal {hrans) to the known input signal {Ibeam), the number 
of counts, or data numbers (DN) were obtained for every input photon. 

i?(A) = ^-p^ [DN phot-']. (4.1) 

-'beam 
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To ensure the response was compatible with the synthetic spectra generated in 
the following section, the wavelength responses of the six instruments were, where 
necessary, corrected for the telescope aperture size (cm^) and for the solid angle 
sky coverage of one pixel (str pix~^). 



The wavelength responses are shown in Figure 4.3, The most striking feature 
is the offset between the instruments' response magnitude. As expected, the 
TRACE response is the lowest of all, mainly due to the low efficiency of its 
detector. Despite being built at the same location, AIA has the highest response. 
The shape of the AIA response is also significantly different to the others. This is a 
testament to both its large aperture size and the developments in technology since 
TRACE was built. The EUVI and EIT instruments have very similar responses 
despite their significantly different aperture sizes. This again, can be attributed 
to the developments in the efficiency of modern detectors. SWAP falls between 
TRACE and EUVI/EIT. A combination of its small aperture size and the nature 
of an APS detector reduces its overall sensitivity. Since an APS chip requires 
individual electronics for each pixel, the overall effective area of the detector is 
reduced. 

Moving through increasing wavelength, it is interesting to note that all instru- 
ments have a gradual increase in their R{X) over the first 5±1 A of the function. 
It is not surprising that AIA and TRACE have similar slopes in this region, as do 
SWAP and EUVI (and EIT, to a certain extent), since these pairs of instruments 
were developed by the same institutes. From 169 A the instruments' responses be- 
gin to deviate. SWAP, EUVI B, EUVI A and EIT have a very sharp rise between 
170 and 171 A, peaking at (1.2, 7.0, 8.5, 11.3) xlO^^^ DN phot"^ cm^ str pix^^ 
respectively. Each of these four responses "plateau" as they remain at values close 
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the peak 171 A value for 6±2 A. The TRACE and AIA responses, however, have a 
significantly different shape. These responses have a much more gradual rise from 
169 A to 172 A, peaking at 2.5x10"^^ and 2.5x10"^ DN phot"^ cm^ str pix'^ 
respectively. Neither TRACE nor AIA have as broad a plateau as the other four 
instruments. The response of AIA falls off at wavelengths immediately after its 
peak at 172 A. It does however, experience a 1 A wide step at 175 A before 
decreasing to an almost flat response above 180 A. The TRACE wavelength re- 
sponse remains close to the peak value for ~3 A and then decreases slowly to an 
almost constant response above 188 A. The decay of the EIT wavelength response 
follows a similar rate to that of TRACE. EUVI and SWAP however, experience a 
secondary increase in their respective response functions between ~184 to 192 A. 



4.2.2 Synthetic spectrum calculation 

The method presented here obtains the temperature response of an EUV imager 
between the temperatures of 2 x 10^ and 10'' K for a range of solar features. In 
order to accommodate both the change in temperature (T) and EM with height 
for these different features, the T and EM values were taken from a DEM curve. 
The amount of emitting plasma in a column dh can be written as: 

EM = I nldh, (4.2) 



As discussed in ^2.2.1, DEM can be written as DEM = J nl{dT/dh)~^. There- 



fore we can redeflne the EM to be: 



EM = / DEMdT (4.3) 

J AT 
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Figure 4.3: Wavelength response of AIA, EIT, EUVI A, EUVI B, SWAP and 
TRACE. 



The EM can therefore be extracted from the DEM curve for a temperature bin 
logioAT = 0.1 wide. 



Figure 4.4 shows the DEM (top) and corresponding EM (bottom) curves used 



in this chapter that have been calculated using CHIANTI. The coronal hole, 
quiet sun and active region curves for both EM and DEM are very similar up 
to a temperature of ~6xl0^ K. This roughly corresponds to the temperature 
at which the coronal hole DEM curve peaks. Following the split, the coronal 
hole curves start to decrease. This means that there is little or no plasma above 
~lxlO^ K present in a coronal hole. The quiet sun curves continue to rise to 
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10* 10^ 10® 10"" 10' 

Temperature [MK] 



Figure 4.4: Differential emission measure and emission measure as a function of 
temperature for coronal hole, quiet sun, active region and flare. 

a peak at ~1.6xl0^ K, at which temperature, they turn over. While the active 
region curves follow the same trend as the coronal hole and active region at 
temperatures less than ~1 MK, it is clear that the active region is approximately 
two orders of magnitude greater. At temperatures greater than 1 MK, the active 
region curves remain approximately constant to almost 10 MK. In all of the 
coronal hole, quiet sun and active region curves, a dip exists at approximately 
1.5x10^ K. The small height of the transition region (at this temperature) has 
a significant effect on the emission levels that originate here. Notice however, 
that the dip is shifted to closer to 1 MK in the flare EM and DEM curves. This 
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is around the temperature of the active region curves' peak. This suggests that 
during a solar flare, the plasma within an active region loop is heated to very 
high temperatures (~10^ K), thus leaving the active region loops devoid of 1 MK 
plasma. 

At each temperature (T) in the range of interest, an EM value was extracted 
from the appropriate EM curve in Figure |4.4 Along with the unique T and its 
corresponding EM value, a generic density value was supplied. The intensity of 
an emission line depends on the population of the upper levels of the atomic 
transition, which itself depends on the plasma density. Understandably, this 
density will change depending on the feature being studied. However, it can be 
taken as a constant within the study of an individual feature. The densities used 
were: 



Coronal hole: 10^ cm'^ (Wilhelm, 2006) 



Quiet sun: 6 x 10^ cm (Young, 2005) 



Active region: 5 x 10 cm (Gallagher et al, 2001) 



Flare: 10^^ cm'^ (Raftery et al. 2009) 



The abundances were taken to be those of Feldman et al. ( 1992 ) and the ionisation 



fractions of Mazzotta et al. ( 1998 ) were used. 



Utilising the parameters listed above, along with a unique T and correspond- 
ing EM (Figure |4.5^, asterisk), an isothermal spectrum was calculated using 



the CHIANTI routines ch_synthetic and make_chianti_spec (Figure 4.5 o 



Each isothermal spectrum was then multiplied by the wavelength response of each 



instrument (Figure 4.5 to give the expected throughput of the instrument: the 
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EM = 2.6 X 10""- cm 
T = 10° K 




180 185 
Wavelength [A] 




10= 10° 
Temperature [K] 



Figure 4.5: (a) Quiet sun EM curve with a single point highlighted at T = 10^ K 
and EM = 2.6 x lO'^^ cm~^. (b) Synthetic isothermal spectrum calculated for the 
conditions selected in panel (a), (c) The SWAP R{X) function, (d) Transmitted 
spectrum from the product of (b) and (c). (e) The SWAP temperature response, 
calculated by integrating the transmitted spectrum over wavelength. The asterisk 
corresponds to the conditions set in panel (a). 
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transmitted spectrum (Figure 4.5 i). The transmitted spectrum, is the number 
of DNs a single pixel is expected to detect in one second at a given wavelength. 



A sample of transmitted spectra are shown in Figure 4.6, taken from throughout 



the temperature range for a flare DEM. Figure 4^ a, shows the transmitted spec- 
trum at a temperature of logio(Te) = 5. It is clear that with no emission lines, 
the continuum is the dominant source of emission at this temperature. Figure 



4.63 shows the transmitted spectrum at logio(Te) = 5.4 where the spectrum is 



dominated by an Ov line at 172.2 A. The Fevill lines between 167 and 169 A 
along with the O and Fe lines at wavelengths above 184 A are all suppressed by 



the instrument response function. Moving up to logio(Te) = 6 in Figure 4.6:, the 
spectrum is now dominated by the characteristic Feix line at 171 A and has sig- 
nificant contribution from a number of Fex lines at slightly higher wavelengths. 
Notice that at this temperature, emission lines completely dominate the spectrum 
and the continuum, now an order of magnitude less than in Figure |4.6^ , has a 
negligible contribution to the spectrum. At a higher temperature, logio(Te) = 6.2 



in Figure 4.6 i, the contribution of the Feix line is not as significant as the contri- 
bution from the series of Fex and Fexi lines, particularly at 174.5 A. Increasing 



the temperature further to logio(Te) = 6.4, as in Figure 4.63, the continuum level 
begins to rise again, although the line emission still accounts for ~93% of the 
total emission. Nixv lines between 174 and 188 A dominate the spectrum at 



this temperature. Figure 4.6" shows the spectrum at logio(Te) = 7.2 where the 
continuum accounts for 41.6% of the emission. Emission lines at this temperature 
consist mainly of highly-ionised species of Fe, specifically Fexxill at 173 A and 
Fexxiv at 192 A. 



It is clear from Figure 4.6 that the relative intensity of the spectrum changes 
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Figure 4.6: Sample transmitted spectra over the temperature range under inves- 
tigation calculated using a flare DEM and the SWAP wavelength response function 
for logio(Te) = 5.0, 5.4, 6.0, 6.2, 6.4 and 7.2. The ion names of the strongest lines 
are also listed. 



significantly with temperature. Therefore integrating the transmitted spectrum 
over all wavelengths results in the expected number of counts a single pixel is 
expected to detect per second at a single temperature, or the instrument's tem- 



perature response, R{T) (Figure 4.53, asterisk). Repeating this process for all 
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10" 10^ 10° 10' 
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Figure 4.7: Instrument response to temperature for 6 EUV imagers in the 171 A 
passband for (a) coronal hole, (b) quiet sun, (c) active region and (d) flare. The 



numbers in the bottom panel refer to the discussion points in ^4.3.1 



temperatures in the range of interest results in the temperature response curve 
for a single instrument and solar feature. 

4.3 Results 



4.3.1 Inter-instrument comparison 



The temperature responses for each of the four solar conditions are shown in Fig- 



ure 4.7 for all six instruments. At first, it appears that aside from the offset in 
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Figure 4.8: The ratio of the AIA and SWAP temperature responses for coronal 
hole (solid), quiet sun (dotted), active region (dashed) and flare (dot-dashed), i, ii 



and in refer to the discussion point in ^4.3.1 



their magnitude, the responses of all six instruments are in good agreement. Con- 
sidering the significant deviations noted in the respective wavelength responses in 



^4.2.1, the correlation between the various temperature responses are remarkable. 



There are some discrepancies worth noting, marked in Figure 4.7 i and Figure 4.8 



with dashed lines. Figure [478] shows the ratio of the AIA and SWAP temperature 
responses since all of the discrepancies involve the AIA temperature response. 

Below the peak temperature, between (2-4) x 10^ K, there is an increase in 
the relative sensitivity of AIA and TRACE compared to SWAP, KIT and EUVI, 



marked as i on Figures 4.71 and 4.8 Since the peak of both AIA and TRACE'S 
wavelength responses occur from ~172 A, they are more sensitive to the Ov 



lines that are formed at 172 A between 0.2 and 0.5 MK (see Figure |4.6p ). This 
also has the effect of shifting the peak response of AIA to 7.5x10^ K compared 
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to 10 K for all other instruments (marked ii on Figures 4.7i and 4.8). The 
other instruments peak at 1 MK because the spectrum is dominated by Fe ix at 
171 A and Fex at 174 A, both with a formation temperature of ~1 MK (Figure 



4.6| :). Since the relative wavelength response of EIT, EUVI and SWAP is high at 



these wavelengths these instruments are sensitive to plasma at 1 MK. TRACE'S 
wavelength response does not peak until 172 A. Therefore it is not sensitive to 
the 171 A line. Unlike AIA however, TRACE'S wavelength is not shifted. This 
is because of the 3 A wide plateau around its peak wavelength response. This 
increases the relative sensitivity of TRACE to Fex lines formed at 174 A that 
AIA does not see due to its steep fall off. 



The third significant deviation is marked number in on Figures 4.71 and 4.8 
At this temperature, 2.5x10^ K, there is significant contribution to the spectrum 



from Nixv lines formed between 172 and 177 A(Figure 4.63). The steep drop off 
of the AIA wavelength response between these wavelengths significantly reduces 
the sensitivity of the instrument to these hot lines. Therefore, AIA's relative 
sensitivity to emission at this temperature is reduced compared to the other 
instruments. 

4.3.2 Solar conditions comparison 



Figure 4.7i shows the temperature response to imaging a coronal hole. Figure 



47p to quiet sun plasma. Figure |47p to an active region and Figure |4.7| d to a 
solar flare. It is apparent that the responses of the instruments to these different 
features are significantly different, as expected. The response to imaging a coronal 
hole has a primary peak at 7.5x10^ K, a temperature only slightly higher than 
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the EM peak in Figure 4A This is not surprising as a comparison between 
Figure |4.6| 3 and c shows. Since the EM in the corona never rises, the EM at 
the temperature of formation of O V dominates over the EM at the formation 
temperature of Feix/x. Therefore the relative contribution of the O lines is 
greater than Fe. Thus, the Ov lines will dominate the transmission spectrum 
and lower the overall temperature. The step that appears at ~10^ K is due to 



emission from continuum at this temperature (Figure 4.6 1). A comparison to the 
quiet sun temperature response shows that the magnitude of the response to a 
coronal hole is smaller. Therefore, longer exposures are required to image these 
regions. Not only has the quiet sun response an overall higher magnitude but 
the emission at 3 MK is almost eight orders of magnitude greater than that of a 
coronal hole. This compares to a factor of 10 at their peaks. This is even more 
exaggerated in the active region case. The "shoulder" that appears at 3 MK in 
the quiet sun is barely noticeable in the active region response as it is dominated 
by contribution from emission at ~7 MK. Despite this however, the main peak at 
1 MK is more than two orders of magnitude greater than the high temperature 
peak. Therefore, it can still be assumed that the majority of emission observed in 
an active region is being emitted at 1 MK. The active region peak response is 10 
times greater then the quiet sun peak. This is not surprising given the distribution 
of plasma in the active region EM curve. The continuum at these temperatures 
has a considerable effect on the emission levels. Finally, the temperature response 
to a flare is the most dramatic of all. A second peak that occurs at 10^ K is now 
the dominant feature of the response and is more than an order of magnitude 
greater than the response at 10^ K. This is due to a combination of strong, highly 
ionised emission lines (Fexxill, Fexxiv) forming at >10 MK and the rising levels 
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Figure 4.9: Normalised R{T) functions calculated by the instrument response 
teams (black lines) for TRACE (a), EIT (b), EUVI A (c) and EUVI B (d). These 
are compared to the R{T) functions calculated in this chapter for each instrument 
for coronal hole (pink lines), quiet sun (red lines), active region (green lines) and 
flares (blue lines). 



of continuum at these temperatures (Figure 4.6| F). 



4.3.3 Comparison to instrument teams' R{T) functions 



As discussed in §6.1 the temperature response functions for these instruments 
have previously been examined by the instrument teams. A comparison to these 
response functions was carried out and is presented in Figure 4^ Each panel 
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contains the R{T) from the instrument team (black line) and the R{T) curves cal- 
culated in this chapter for the four solar conditions. Beginning with the TRACE 



plots in panel (a), it is clear that the R{T) calculated by Brooks & Warren (2006) 
is in reasonable agreement with our flare R{T) at temperatures less than ~10^ K. 
At the peak temperature of 1 MK, the Brooks response is most similar to our 
active region response, although it appears to be shifted slightly towards cooler 
temperatures in conjunction with our coronal hole R{T). At high temperatures, 
the correlation breaks down, with the Brooks response almost three orders of 
magnitude less than our flare response. This is a direct result of the explicit 
treatment of the four differing plasma types. The Brooks R{T) attempts to in- 
corporate all types of conditions into one single response function and in doing so, 
is removing the significance of the high temperature peak that may be produced 
during a solar flare. 

The EIT teams temperature response function (panel b) shows that the re- 
sponse at 10^ K is less than that of a flare, as is the case with TRACE. This is 
most likely to do with the fact that when this analysis was conducted, the con- 



tinuum levels were not as well understood (Feldman et a/., 1999). Like TRACE, 
there appears to be a broadening of the response across our coronal hole and 
active region peaks at 1 MK. The EIT response to temperatures above 2x10^ K 
is not available. 

The responses of the two EUVl instruments shown in panels (c) and (d) 
are almost identical and so we will discuss them simultaneously. Like EIT, the 
EUVI R{T) function at 10^ K calculated by the instrument team is not as large 
as that of our flare and the response at 1 MK is again, broadened across our 
coronal hole and active region curves. Similar to TRACE, the high temperature 
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continuum is greatly underestimated due to combining the responses of multiple 
solar conditions into one R{T) function. 

4.4 Conclusions 

This chapter studies the temperature response of the SWAP EUV imager, along 
with five similar instruments in the 171 A passband. The temperature responses 
were calculated for four different solar features: coronal holes, quiet sun, active 
regions and flares. Besides an expected offset in the magnitudes of the tempera- 
ture responses, the six instruments were found to agree remarkably well, despite 
significant deviations in their wavelength response functions. Any differences 
that occurred in the temperature responses were as a result of the shape of the 
wavelength responses. 

The ability to compare the responses of individual instruments is important 
to understanding both similarities and differences in co-ordinated observations. 
Presently there exist comparisons between the TRACE and EIT instruments. 
This is understandable since these two instruments have been at the forefront of 
EUV imaging until very recently. However, with the introduction of four new 
imagers in almost as many years there is a clear need to expand the instru- 
ment comparison. While the individual response functions are available in their 
respective instrument papers, different units, methods and initial conditions in 
their calculation make comparison difficult. In addition, the instruments' tem- 
perature response functions are, in general, presented for only one solar feature 
- usually quiet sun. This chapter has presented a solution to this problem by 
combining an in depth study of multiple instruments and different solar condi- 
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tions. A comparison to the existing temperature response functions reveals that 
our results are in very good agreement with the response curves of the instru- 
ment teams at temperatures below ~3xl0^ K. Above this temperature, the high 
temperature continuum contribution that is found, in particular during a solar 
flare, is highlighted though our method of isolating the individual plasma types. 

The facility to analyse the temperature response to different types of plasma is 
as important as the inter-comparison of the instruments themselves. For quiet sun 
and active regions, the interpretation of emitting plasma is reasonably straightfor- 
ward and well known in the community. Coronal hole and flare emission however, 
cannot be interpreted in the same way. The dominance of continuum and cool O 
emission lines in coronal holes creates a shift in the maximum temperature sensi- 
tivity towards lower temperatures. Likewise, the signiflcant free-free continuum 
and hot Fe lines in solar flares results in a multi-peaked temperature response, 
with emission above 10 MK dominating. 

This chapter presents a new, more appropriate method of calculating the 
temperature response of EUV imagers to different types of solar conditions. This 
method takes account of the variation in EM with temperature for coronal holes, 
quiet sun, active regions and flaring conditions. It can be applied to multiple 
instruments, thus facilitating their direct comparison. In addition, the method 
can be easily adapted to investigate other passbands besides 171 A and new 
instruments in the future. This work will facilitate a more precise understanding 
and interpretation of EUV images. 
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Chapter 5 

Multi-wavelength observations 
and modelling of a solar flare 

This chapter discusses observations of a confined, C-class solar flare in light of 
the standard flare model. The evolution of temperature, emission measure, energy 
loss and footpoint velocity were investigated from both observational and theoret- 
ical perspectives. These properties were derived by analysing the systematic cool- 
ing of emission through the passbands of RHESSI (>5 MK), GOES (5-30 MK), 
TRACE 171 A (1 MK), and CDS (r^ 0.03-8 MK). Comparing these observations 
to a 0-D hydrodynamic model, EBTEL, enabled the author to diagnose gentle and 
explosive chromospheric evaporation and analyse the heating and cooling mech- 
anisms involved in the event. This is the first eoctensive study of the evolution 
of a canonical solar flare using both spectroscopic and broad-band instruments, in 
conjunction with a 0-D hydrodynamic model. While our results are in broad agree- 
ment with the standard flare model, the simulations suggest that both conductive 
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and non-thermal beam heating play important roles during the impulsive phase of 
at least this event. This work was published by Raftery, Gallagher, Milligan & 
Klimchuk, Astronomy & Astrophysics, 4.94, 1127, 2009. In the analysis of this 
event, significant assistance was given by Dr. Ryan Milliagan for the preparation 
and interpretation of the CDS and RHESSI observations. The theoretical model, 
EBTEL, was written by Dr. James Klimchuk and considerable assistance in the 
use of this model was offered by him. Guidance and advice on this project was 
given by Dr. Peter Gallagher. 
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5.1 Introduction 

The temporal evolution of most solar flares can be divided into two distinct 
phases. During the impulsive phase, temperatures rise to >10 MK via direct 
heating below the reconnection site in the corona and/or chromospheric evap- 



oration due to accelerated particles (Kopp & Pneuman 1976). Chromospheric 



evaporation, discussed in detail in §1.2 can be classifled in one of two ways: 



explosive or gentle (Fisher et ai, 1985c Milligan et al., 2006a|[b ). Explosive 



evaporation occurs when the flux of non-thermal particles impacting the chro- 
mosphere is greater than a critical value of approximately 3x10^° ergs cm~^ s~^. 
In this case, the chromosphere cannot dissipate the absorbed energy efliciently 
enough. The plasma is forced to expand into the corona as hot upflows of hun- 
dreds of km s~^ and simultaneously into the chromosphere as cooler downflows 
of tens of km s~^. Beam driven gentle evaporation occurs when the non-thermal 
flux is less than ~10^'' ergs cm~^ s~^. Under these circumstances, the chromo- 
sphere is efficient in radiating the absorbed energy and plasma expands slowly 
(tens of km s~^) upwards into the loop. Gentle evaporation can also be driven 
by a downward heat flux from the corona in what is known as conduction driven 
gentle evaporation: conduction fronts that propagate from the hot looptop to 
the footpoints heat the chromospheric plasma and drive the gradual expansion of 
plasma upwards into the loop. 

Once the energy release has ceased, the hot plasma returns to its equilibrium 
state during the decay phase. The cooling process begins with thermal conduction 
as the dominant energy loss mechanism due to the high temperatures present 



see Equation 1.11 ). As the temperature decreases and the radiative loss function 
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begins to increase (Figure 1.16) and radiative cooling becomes the primary energy 



loss mechanism (Culhane et al, 1970). Finally, the "evaporated" material drains 
back towards the solar surface, returning the system to equilibrium. 

There have been a wealth of studies that focus on the hydrodynamic mod- 



elling of these heating and cooling mechanisms (e.g. Antiochos & Sturrock, 1978 



Fisher et al. 


1985c 


Doschek et al. 


1983 


Cargilll 11993 Klimchuk & Cargill 


2001 



Reeves & Warren 


2002 


Bradshaw & Cargill, |2005 ; Klimchuk 


2006 


Warren & 



Winebarger, 2007 Sarkar & Walsh, 2008). For example, Reale (2007) conducted 



an analysis of the details of stellar flares using the Palermo-Harvard theoreti- 



cal model (Betta et al, 1997 Peres et al, 1982). This paper fully describes 
the cooling timescales and plasma parameters of flares in terms of their phases, 
including an investigation of the thermal heating function. However, these, and 
most other theoretical results were not compared to observations. The majority of 
investigations that make this comparison concentrate on broad-band instruments 



and utilise very simple models. For example Culhane et al. (1994) compared 



Yohkoh observations to an over-simplified power-law cooling curve. [Aschwanden 



& Alexander (2001 ) compared broad-band observations to a model similar to that 



of Cargill (1994) that considers a purely conductive cooling phase followed by a 



purely radiative cooling phase. Vrsnak et al. (2006) conducted a similar study, 
again concentrating on broad-band observations and a simple, independent cool- 



ing mechanism model. Teriaca et al. (2006) conducted multi- wavelength analysis 
of a C-class flare, incorporating observations from CDS, RHESSI, TRACE and 
ground based detectors. The cooling timescales were obtained using the simple, 
independent-cooling Cargill model. 

In this chapter, we attempt to improve upon previous studies by comparing 
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high resolution observations over a wide range of temperatures to a theoreti- 
cal model. Observations of a GOES C-class solar flare were made with several 
instruments, including TRACE, RHESSl, GOES-12 and CDS. There are many 
advantages to using spectroscopic data in conjunction with broad-band observa- 
tions. The identification of emission lines are for the most part, well documented 
and therefore individual lines can be isolated for analysis. Also, material as cool 
as ~30,000 K can be observed simultaneously with emission at ~8 MK. Fur- 
thermore, it is possible to carry out velocity, temperature and emission measure 
diagnostics over a wide range of temperatures for the duration of the flare, sig- 
nificantly improving the scope of the analysis undertaken. These observations 
were compared to a highly efficient 0-D hydrodynamic model, EBTEL. The com- 
bination of this extensive data set and the new modelling techniques enabled the 
author to carry out a comprehensive investigation of the heating and cooling of 
a confined solar fiare. 



5.2 Observations and data analysis 

This investigation concentrates on a GOES C3.0 fiare that occurred in active re- 
gion NOAA 9878 on 26-March-2002. Located close to disk centre (-92", 297") the 
event began at ~15:00 UT. The CDS observing study that was used (flare_ar) 
focused on five emission lines spanning a broad range of temperatures. The rest 



wavelengths and peak temperatures of the lines are given in Table 5.1, where the 



quoted temperatures refer to the maximum of the contribution function. The 



contribution functions for the five CDS lines are shown in Figure 2.5 Each CDS 



raster consists of 45 slit positions, each 15 s long, resulting in an effective cadence 
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Table 5.1: Rest wavelengths and temperature of emission lines and bandpasses 
used in this study. 



Ion 


Ao [A] 


Temperature 


[MK] 


He I 


584.45 


0.03 




Ov 


629.80 


0.25 




Mgx 


625.00 


1.2 




Fexvi 


360.89 


2.5 




Fexix 


592.30 


8.0 




TRACE 


171 


1.0 




Instrument 


Range 


Temperature 


[MK] 


GOES 


0.5-4A & 1-8A 


5 - 30 




RHESSI 


3 keV-17 MeV 


>5 





of ~11 minutes. The slit itself is 4"xl80", resulting in a 180"xl80" field of view. 



Figure 5.1 shows the evolution of the fiare in multiple wavelengths. The top 
two rows show the looptop source observed in RHESSI 6-12 and 12-25 keV energy 
bands, with the Ov and Fexix (15:09 UT) contours overplotted. The next four 
rows of this figure show emission observed in Fexix, Fexvi, TRACE 171 A and 
O V. At ~15:00 UT, before the main impulsive phase of the fiare began, evidence 



of low level Fexix loop emission was observed (first Fexix image in Figure 5.1). 
By ~15:09 UT, the footpoints were seen in O v while the Fexix loop top emission 
continued to brighten. At 15:16:40 UT, when RHESSI emerged from eclipse, a 
thermal looptop source was observed in both 6-12 and 12-25 keV energy bands. 
By 15:21 UT the loop was emitting predominantly at ~8 MK and a bright "knot" 
was seen at the top of the loop. The area of the knot is marked with a dotted box 
on the Fe xix images. Such features have been observed in the past, although they 



have not been readily explained (e.g. Doschek & Warren, 2005). By ~15:32 UT 



the loop was emitting mainly in Fexvi at 2.5 MK and had cooled to <1 MK, 
into the TRACE and Ov passbands by ~15:44 UT. 
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Figure 5.1: RHESSI 6-12 and 12-25 keV images are shown in the top two panels 
with O V and Fe xix contours overplotted in red and yellow respectively. Fe xix 
(8MK), Fexvi (2.5 MK), TRACE/171 A (1.0 MK) and Ov (0.25 MK) images are 
shown in the next four rows. The dotted box represents the loop "apex" used for 
comparison to EBTEL. The bottom panel shows the GOES 1-8 A, RHESSI 3-6, 
6-12 and 12-25 keV lightcurves. RHESSI was in eclipse until 15:15 UT and passed 
through the South Atlantic Anomaly between 15:35 UT and 15:58 UT. The vertical 
dotted lines (and corresponding arrows) on the GOES plot represent the start and 
end times of the CDS rasters above. 
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The lightcurves plotted at the bottom of Figure show the GOES 1-8 A 
and RHESSl 3-6, 6-12 and 12-25 keV lightcurves. Although RHESSl was in 
eclipse for the majority of the impulsive phase (up to ~15:15 UT), the observed 
continued rise of the 6-12 keV lightcurve after emergence from night implies that 
the peak of the soft X-rays was observed. However, while a hard X-ray component 
(12-25 keV) was observed, we believe that the HXR peak occurred before this 
time. 

5.2.1 Temperature and emission measure 



The RHESSl spectrum, shown in Figure 5^ was analysed for one minute be- 
tween 15:16:30 UT and 15:17:30 UT. Since it is believed that this was after the 
the HXR peak, this may mean that the maximum temperature was not in-fact 



observed. Following previous studies (e.g. Saint-Hilaire &; Benz, 2002), the data 



were fitted with an isothermal model at low energies and a thick-target model 



up to ~30 keV (see ^3.5.2 and Equation 2.27). The thick target component 
yielded a low energy cutoff of 17 keV and a power law index of 8.2. The isother- 
mal fit to lower energies resulted in temperature and emission measure values 
of ~13 MK and ~lxlO^® cm~^ respectively. A non-thermal electron flux of 
~7xl0^ ergs cm^^ s^^ was also calculated by approximating the footpoint area 
from He I and O v observations. Since it is probable that the HXR peak occurred 
before this time, this non-thermal electron flux must be assumed to be a lower 
limit to the peak flux. The filter ratio of the two GOES passbands produced 



the temperature and emission measure of the event as described in [:3.4, giving a 



peak temperature of 10 MK and an emission measure of 4x10 cm . 
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4 5 6 7 8 9 10 

Photon Energy [keV] 

Figure 5.2: The top panel shows the RHESSI photon spectrum between 
15:16:30 UT and 15:17:30 UT. The data points (crosses) were fit with a thermal 
Maxwell distribution (dotted line) and a thick target, non-thermal model (dashed 
line). Combining these model parameters with the background function (dot-dash 
line) resulted in the overall fit to the spectrum (solid line). The residuals of the 
spectrum fit are shown in the bottom panel. 



Four of the CDS emission lines (excluding He I as it is optically thick) and 
TRACE 171 A images were integrated over the area described by the bright 
"knot" mentioned above. Although the loop is believed to consist of multiple 
magnetic strands, it was assumed that the majority of strands within this region 
were heated almost simultaneously during the impulsive burst. A small number 
of strands can be heated before or after this time, producing a multi-thermal 



plasma, however from studying Figure 5.1, this small region was assumed to be 



approximately isothermal at any one time. The lightcurves of this region are 
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Lightcurves 26-March-2002 
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Start Time (26-Mar-02 14:50:00) 



Figure 5.3: Lightcurves from GOES, Fexix, Fexvi, Mgx, TRACE 171 A and 
O V. Overplotted on the data points are sphne interpolations (black lines) and the 
lightcurves predicted by Equation |5.1| (red lines) . The vertical lines represent the 
times adopted as the peak of the lightcurves. 
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shown in Figure |5.3[ In order to ascertain the thermal evolution of the flare, 
a temperature was assigned to each lightcurve. This was achieved through the 
fitting of a spline interpolation to the data points and comparing them to forward 
modelling results of various emission lines and of the TRACE 171 A passband 



(Figure 5.3). The emission measure and temperature evolution of the event as 



predicted by EBTEL (top two panels of Figure 5.5) were used to calculate the 



value of the contribution functions G{Ti)ion shown in Figure 2.5 and emission mea- 



sure EMi at each timestep i. Their product resulted in the predicted lightcurve 



for a particular ion at a given time, as in Equation 5.1 



h = G(T,)_ X EM, (5.1) 



As Figure 5^ shows, the peak of the predicted lightcurves agree very well 
with the observed lightcurves. The predicted lightcurves, in general, reach their 
maximum intensity at the time {t) the contribution function is maximised (i.e. 
t{G(T)max) = t{Imax))- Therefore, the temperature of the peak of the contribu- 
tion functions were were assigned to the time of the maximum intensity. The 
poor cadence of the CDS instrument was such that the time of maximum inten- 
sity could only be estimated by a spline fit to the datapoints (black lines. Figure 



5.3D . 

The predicted lightcurves for Fexix and Fexvi are in very good agreement 
with the observed lightcurves. This is most likely due to the high intensity of 
these lines during a fiare. The predicted lightcurve for Mgx however, requires 
some explanation. It is clear that there are two peaks in intensity for MgX, one 
each side of the spline peak. By 15:30 UT, the time of the first predicted peak. 
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the temperature (shown in the top panel of Figure 5.5) had already reached its 
maximum and had begun to pass back into the band of sensitivity of the Mgx 
contribution function. At this time, the levels of EM were raised to a signifi- 
cantly higher level compared to during the impulsive phase, when the temperature 
initially rose through the Mgx G{T). The combination of the increasing G{T) 
and elevated EM resulted in a broad peak in the predicted lightcurve around 
15:30 UT. A dip is then seen in the predicted lightcurve at ~15:40 UT. This is 
as a result of the falling EM at this time and the leveling off of the contribution 
function at ~2xl0® K. As the contribution function rises again to its peak at 
1 MK, the lightcurve too rises to a sharp peak. Although we have ascertained the 
theoretical time of the 1 MK peak, the poor cadence of the instrument restricts 
us from being as exact with the observation. 

Similar to Mg X, the O v lightcurve peak was difficult to establish. This may 
be due to the high levels of O v in the surrounding areas, along with the low 
cadence of the instrument. In addition, the T and EM values used for predicting 



the lightcurves were obtained from Figure 5.5 It can be seen that the EBTEL 
temperature is not in good agreement with the Ov data point, predicting that 
the loop will reach this temperature a short time later. 

The instrument response function used in the forward modelling of TRACE 
were those calculated in Chapter |4j The first small peak at 15:13 UT was repro- 
duced using the flare response function. However, in that case, the main peak 
was not predicted and so we find it inappropriate to use the fiare response to 
predict the temperature of TRACE during the decay phase. The active region 



response curve (Figure 4.7 panel c) however, predicts very nicely the main late 
phase intensity peak along with a smaller peak at 15:28 UT (which also corre- 
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sponds to the broad peak seen in Mgx formed at a similar temperature) as shown 



in the fifth panel of Figure 5.3 In 4.3.2 it was established that when imaging 



an active region, TRACE is most sensitive to plasma at 1 MK. For this event, 
by the time the TRACE intensity begins to increase, the main phase of the fiare 
has passed and the system is similar to a bright active region loop. Therefore, 
assigning a temperature of 1 MK to the time of the TRACE lightcurve maximum 
is justifiable. 



The emission measure was calculated by inverting Equation 5.1 and assuming 
an isothermal plasma at any given time. We have not included TRACE in the 
analysis of EM due to its broad and complex response function (see Chapter |4] 
for further details). For the uncertainty in the CDS EM, there are a number of 
factors to consider. These include uncertainties in the intensity of the line, the 



contribution function, the volume element dV assumed for Equation 5.1 and the 



CDS calibration. While the uncertainty in measuring the line intensity is small 



for strong lines such as those used in our study, (typically ~10%; Del Zanna et al. 



2001), a consideration of the contribution functions FWHM yields an uncertainty 



in the EM of up to 30%. In addition, the CDS calibration is known to be good 



to within 15-20% (Brekke et al., 2000). Considering these factors, the combined 



photometric error was taken to be 50%. 
5.2.2 Velocity 

The relative Doppler shifts at both footpoints were calculated for the duration 
of the fiare using the five CDS emission lines. The centroids were calculated 



following ^3.1.3.1 and were corrected for both heliocentric angle and an average 
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Figure 5.4: The line profiles of the Fe xix emission line for the right footpoint dur- 
ing the impulsive phase (a) and the decay phase (b) . Note the increasing intensity 
of the Fexii blend at ~592.8 A (Del Zanna & Mason, 2005) during the decay phase. 



The solid vertical line represents the rest wavelength while the dashed vertical line 
is its centroid. 



inclination of 44° to obtain the rest wavelengths listed in Table [5T| As mentioned 
previously, it is not possible to obtain absolute velocities using CDS and so these 
should be interpreted as relative velocities. Figure |5.4 shows the Fexix line 
profiles for the right footpoint during (a) the impulsive phase and (b) the decay 
phase. 



5.2.3 Modelling 

The EBTEL model simulates the evolution of the average temperature, density, 
and pressure along a strand at any given time ( ^1.2.2.2 ). The flare was mod- 



elled as a single, monolithic loop. This was justifiable since the observations 
suggest that most of the strands are heated in approximately the same way and 
at the same time. The standard EBTEL pre-fiare conditions included a temper- 
ature of 0.3 MK, an initial density of 5x10^ cm~'^ and an emission measure of 
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4x10^^ cm ^. The observations described in this chapter were, where possible, 
used to constrain the model's input parameters. The ranges of acceptable parame- 



ters studied, along with the observed values are shown in Table |5.2[ The method 
presented here does not use any rigorous statistical fitting methods. However, 



Adamakis et al. (2009) utihsed a Bayesian fitting technique to better identify the 



best-fit model parameters. The results of their study are also shown in Table 



5.2 for comparison. In general, the statistically robust parameters are in good 



agreement with those published in Raftery et al. (2009) with the exception of the 
direct /non — thermal ratio which is discussed below. The loop length was de- 
termined from magnetic field extrapolations of the region (P. A. Conlon, private 
communication). Since it was assumed that the heating function is associated 
with the HXR burst, the majority of which was not observed, the shape of the 
heating function was inferred both from previous observations of HXR bursts and 
the slow rise of the GOES SXR lightcurve. Thus, the most appropriate heating 
function was deemed to be Gaussian in shape. This choice was later justified by 



Adamakis et al. (2009) who proved statistically that a Gaussian distribution is 



is in better agreement with observations than a half-Gaussian distribution. The 
amplitude of the non-thermal electron fiux was constrained by the lower limit cal- 
culated from RHESSl observations and the width was inferred from the derivative 
of the SXR fiux (g]2) |Neupert| |T968| [Zarro fc Leinenl pSS) ). While the direct 
heating rate was not constrained by observations, it was assumed to have the 
same width as the non-thermal heating fiux and to occur at the same time. The 



range of parameter values shown in Table 5.2 correspond to the maximum and 
minimum values that produce an acceptable fit to data. The ratio of the heating 
components (i.e. direct to non-thermal) is also shown for the best fit parameters. 
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Table 5.2: Input parameters used for EBTEL simulation for work shown in this 



Chapter and work carried out by Adamakis et al. (2009). The parameters were 



constrained by data when possible and the ranges of parameters investigated are 
shown. 



Parameter Observed EBTEL Adamakis 



Loop half-length [cm] 


3x10^ 


(3 ±0.2) X 10^ 


(3.5 ±0.4) X 10^ 


Non-thermal flux 








- Amplitude [ergs cm~^ s~^] 


7x10^ 


5 X 10^^^ 


(3.2 ±0.4) X 10^ 


- Width [sec] 


~100 


100 ± 50 


92 ±20 


- Total [ergs cm~^] 


~1. 7x10^2 


2.5 X 10^°^^ 


3 ± 4 X 10^ 


Direct heating rate 








- Amplitude [ergs cm~^ s^^] 




0.7 ±0.3 




- Width [sec] 




100 ± 50 




- Background [ergs cm~'^ s~^] 




< 1 X 10-6 




- Total [ergs cm~^] 




175 ± 150 




Direct /non-thermal heating 


(best fit) 


~ 4 


23 ± 16 



Note the discrepancy between the Adamakis et al. (2009) heating ratio and the 
value used in this chapter result from the exclusion of background heating levels 
in the ratio calculated by Adamakis and the inclusion of the background heating 
here. 

In order to highlight the importance of the effect of cooling by conduction 
and radiation simultaneously, the data were also compared to the Cargill model 



( ^1.2.2.1D . Following [Antiochos fc Sturrock] ([1976]) , |Cargill| p993l [1994| ) presented 
a model that considered a flare that is cooling purely by conduction for a time 



Tc, followed by purely radiative cooling for a time r^., as given in ^ 1.2.2.1 The 
cooling times, along with and T^,, the time and temperature at which the 
cooling mechanisms change were calculated for the flare and compared to the 
results of the 0-D EBTEL model. 
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5.3 Results 

Combining the observations from the different instruments used for this study 
with results from EBTEL, the heating and cooling phases of this flare can be 



comprehensively described. These results are presented in Figs. 5.5 and 5.6 
These figures show the evolution of the flare through the dependence of temper- 
ature, emission measure, energy losses and velocity. 

5.3.1 Comparison of model to data 

The EBTEL parameters were allowed to vary within the limits listed in Table 



5.2 until a good fit to the cooling phase data points was found. The cadence of 
CDS was such that it did not observe any significant intensity deviations during 
the impulsive phase and since RHESSl was eclipsed during the impulsive phase, 
the only reliable observations available were for the decay phase. Therefore, no 
observations could be used to constrain the model during the rise phase of the 
flare. 



The top two panels of Figure [5~5] describe the evolution of the flare temperature 
and emission measure from both observations (data points) and model (solid 
line). The conductive and radiative loss curves generated by EBTEL for the flare 
are shown in the third panel of Figure |5.5[ Conduction was found to dominate 
initially, with radiation becoming prevalent for the remainder of the decay phase. 



This is consistent with previous observations (e.g. Aschwanden & Alexander 



2001[ [Culhane et al] [1994j). Both the Cargill and EBTEL simulations found 
conduction to dominate for the first 200 - 400 s of the decay phase, with radiation 
dominating for the remaining ~4000 s, referring to Tc and Tr respectively. The 
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time r* at which Tc ~ is ~15:24 UT in both cases. However, the temperature at 
which this occurs, T* = T(r*), was found to be ~12 MK and ~8 MK according to 
Cargill and EBTEL respectively. This is not surprising since EBTEL uses both 
radiation and conduction to remove heat from the corona while Cargill uses only 
one of these mechanisms at any given time. Therefore the overall efficiency of 



heat removal will be greater for EBTEL than for Cargill. The fourth panel of 5.5 



shows the GOES 1-8 A lightcurve for context. The last panel shows the velocities 
at the loop footpoints, calculated following the analysis in §5.2.2 where negative 
velocities correspond to upfiowing (blueshifted) plasma. The flow velocity at both 
left and right footpoints are shown for the coolest and hottest lines - He I and 
Fexix respectively while, for clarity, only the right footpoints for the remaining 
three lines are shown. The Mgx velocity simulated by EBTEL is represented 
by the thick black line. The simulations are in reasonable agreement with the 
observations. Upfiows are of course predicted during the evaporation phase and 
downfiows are predicted during the draining phase. However, the magnitudes are 
generally larger than those observed in the Mgx line. This likely to be as a result 
of the treatment of the HXR beam by EBTEL. By assuming all of the beam 
energy is transferred into upfiowing plasma, the velocity attained will be greater 
than the velocity observed which, in reality, will be driven by only a portion of 
the beam energy. The small peaks seen at later times in the simulated velocity 
curve are a result of the piece-wise continuous form used for the radiative loss 
function. 



Figure 5.6 shows the evolution of the flare through the interdependence of 
emission measure and temperature. The data points obtained during the analysis 



described in ^5.2.1 were computed at the same time for any one emission line or 
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bandpass. Figures |5.5| and |5.6| show four phases of the flare evolution: the pre- 
fiare heating of the plasma (A) followed by evaporation of hot plasma (B), cooling 
(C) and draining (D). 

5.3.2 Flare phases 

A) 14:45-15:10; Pre-flare phase: For the majority of this phase, the EBTEL 



parameters remained at quiet Sun values, as phase (A) of Figure 5.5 shows (i.e. 



between 14:45 and 15:10 UT). At 15:07 UT the EBTEL temperature and emission 



measure began to rise. Figure 5.6 shows the steep temperature gradient and the 



initial gradual rise in emission measure. However, as the fourth panel in Figure 



5.5 shows, the GOES soft X-rays began to rise slowly before this. Since the 
impulsive phase of the flare was not constrained, this pre-flare heating was not 
modelled. At ~15:00 UT, a small amount of Fexix emission was seen in the loop 



(Figure 5.1, first Fexix image). In addition, velocities of 90±16 km s ^ observed 



in Fexix can be seen in the bottom panel of Figure 5^ while all of the cooler 
lines remain at rest. This could be evidence for pre-flare gentle chromospheric 
evaporation in a small number of strands before the HXR burst. 

B) 15:10-15:17; Impulsive phase: During the impulsive phase of a flare, the 
standard model predicts the propagation of non-thermal electrons to the chromo- 
sphere where they heat the ambient plasma, causing it to rise and fill the loop. 
Upflows of 81±16 km s~^ in Fexix and simultaneous cool downflows of 16 and 
33±16 km s~^ He I and O V respectively were observed and shown in the bottom 



panel of Figure 5.5 A non-thermal electron flux of ~7xl0^ ergs cm ^ s ^ was 



determined between 15:16:30 and 15:17:30 UT from RHESSI observations. This 
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(B) (C) 
RHESSI 




15:00 15:15 15:30 15:45 16:00 

Start Time (26-Mar-02 14:47:02) 



16:15 



16:30 



Figure 5.5: Comparison of observations (data points) to EBTEL model (solid 
lines, except panel 4). From top to bottom, the panels describe the temperature, 
emission measure, loss rate, GOES SXR flux and velocity profiles of spectral lines 
of the left and right footpoints of the flare as a function of time. The adopted 
heating function for EBTEL is shown overplotted on the first panel. The data 
points in the second panel correspond in time to the ones above in the temperature 
panel. The EBTEL modelled velocity (thick black line in the bottom panel) is for 
Mgx. The dotted vertical lines correspond to the flare phases (A)-(D) explained 
in ^5X21 
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Figure 5.6: This shows the dependence of emission measure on temperature for 
both model and data. The different phases of the flare are marked (A)-(D). Over- 
plotted are the emission measure data-points as a function of their temperature. 

is slightly lower than the 3x10^° ergs cm~^ s^^ required to drive explosive chro- 



mospheric evaporation (Fisher et al, 1985c Milligan et al. 2006a). This may 



be due to the uncertainty in estimating the area of the HXR beam, although it 
is most likely due to RHESSI being in eclipse during the main release of hard 
X-rays. It is assumed that 7x 10^ ergs cm^^ s^^ is a lower limit to the peak HXR 
flux for this event. 



C) 15:17-15:24; Soft X-ray peak: The top panel of Figure 5.5 shows the tem- 



perature has peaked and begun to fall and that the emission measure and SXRs 



were at a maximum in this phase. As Figure 5.1 shows, the non-thermal flux 
levels are signiflcantly reduced at a value of ~9xl0'' ergs cm""^ s^^. Conduc- 



tion however, was highly efficient at this time (third panel of Figure 5.5). The 
blueshifts observed in Fexix (25±16 km s~^) during this phase may be as a 
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result of conduction driven gentle chromospheric evaporation. The lack of any 
significant downflowing plasma at this time corroborates this theory. 

D) 15:25-16:30; Decay phase: This phase is dominated by radiative cooling, 



as seen in the third panel of Figure 5.5 Velocities in Fexix were returning to 



quiet Sun values. If the upfiows observed at 15:32 UT in Fexix is not as a 
result of continued conduction driven gentle evaporation, they may be explained 
by the line profile containing components from evaporating strands that were 
heated after the main loop bundle. Between approximately 15:45 and 16:20 UT 
Mgx, Ov and He I showed downfiows of up to 42±16 km s^^. This implies loop 



draining was occurring (Brosius 2003). By the end of the simulation, all of the 



parameters had returned to quiet sun values. 



5.4 Momentum Balance 



Following the approach taken by Teriaca et al. (2006), the momentum, p = mv 



of the upflowing and downflowing plasma was calculated. Mass density p can be 
expressed as the product of the mean molecular mass, the proton mass and the 
total particle density: 

p = firripntot (5.2) 

where = 1.67 x lO'^'^g. fi = 1.27 for a neutral plasma and 0.61 for a fully 
ionized plasma and ritot can be expressed in terms of the chromospheric and 
coronal densities. Since m = pV, the mass can now be expressed as 

m = finipntotV = pnipntotAAh (5-3) 
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where A is the cross sectional area of the loop and Ah is the height through which 
the plasma can flow. The footpoint area was estimated from Ov images taken 
at the peak of the flare and was found to be ~ 10^^ cm^. This is consistent with 



a typical footpoint area (e.g. Canfield et al, 1987). Following Equation 5.3, the 
momentum can be expressed as: 



p = finipntotAvAh 



(5.4) 



for average upflow or downflow velocity v. 

For chromospheric evaporation into the chromosphere, we assume the total 
density can be taken to be 1.1 times the typical chromospheric density where 
rich ~ 10^'^ cm^^ (see e.g. Figure 1.2). From Figure 1.15 it is clear that the 



high density found in the chromosphere restricts the depth to which downflowing 



plasma can reach. This leads to an estimate of Ah ^ 200 km (Abbett & Hawley 



1999 Teriaca et al. , 2006 ). With average downflow velocities of ~20 km s during 



the impulsive phase of this flare, the momentum of the downflowing plasma was 
calculated to be 9 x 10^°=^^ g cm s"^ 

Chromospheric evaporation upwards into the corona can occur along approx- 
imately the entire loop length. Therefore, in the coronal calculation, we take 
Ah = 3 X 10^ cm. The coronal density can be approximated as ritot = 1-91 x rig 
where Ue ~ ^^/{EM/A x Ah) = 1 x 10^" cm~^. Considering an average upflow 
velocity of ~80 km s~\ the momentum of upflowing plasma was calculated to 
be 5 X 10^*^^^ g cm s~^. These values are consistent with momentum calculations 



carried out for flares of similar sizes in e.g. Canfleld et al. (1987); Teriaca et al. 
(|2006|. 
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5.5 Conclusions and discussion 

This chapter compares a flare observed with CDS, TRACE, GOES and RHESSI. 
Early in the impulsive phase of the flare, evidence of 8 MK emission and 95±16 km s~^ 
upflows suggest the pre-flare heating of the loop. During the impulsive phase, 
hot upflowing plasma at velocities of 81±16 km s~^ and cool downflows of up 
to 33±16 km s~^ imply explosive chromospheric evaporation, in accordance with 



previous observations such as Milligan et al. (2006a); Teriaca et al. (2006). Up- 



flowing plasma at velocities of 24±16 km s~^ were observed around the time of 
the SXR peak in Fexix, along with negligible (6±16 km s~^) He I upflows. Since 
conduction was found to be highly efficient at this time, this may be evidence for 
conduction driven gentle chromospheric evaporation. 

As expected, once the temperatures within the loop had fallen and tempera- 
ture gradients were significantly reduced, the conductive loss rates fell consider- 
ably. Thus, the late decay phase of the fiare was dominated by radiative cooling. 
Late in the fiare (post 15:45 UT), the draining of plasma from the looptop was 
observed in He I, Ov and Mgx. 

During the explosive chromospheric evaporation phase, the momentum of the 
moving plasma was calculated. The upfiowing plasma momentum was found to be 
4 X 10^°^^ g cm s~^ while for downfiowing plasma it was 9 x lO^''^^ g cm s~^. Not 
only do these results show that the momentum was conserved within the loop, 
they are also consistent with previous studies of fiares of similar GOES class. 



Canfield et al. (1987) for example found the upfiow momentum of a small C2 



class fiare to be ~2xl0^° g cm s ^ and the corresponding downfiow momentum 



to be ~6xl0^° g cm s ^ while Teriaca et al. (2006) found the momentum of a 
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GOES C2.3 flare to be (1 ± 5) x 10^*^ g cm in both directions. The balance 
between these values means momentum is conserved within the system during 
chromospheric evaporation. Since the momentum of the material appears to be 
consistent (at least for C-class flares), the momentum of the plasma may lead 
to furthering the understanding of the force applied to the plasma due to the 
electron beam. 

The observations of the flare were compared to the 0-D hydrodynamic model, 
EBTEL. The comparison of the simulations and observations were confined to 
the decay phase of the event, due to data limitations. It is believed that for 
future analysis, the investigation into different direct heating models (e.g. the 
width of the function and the time at which it peaks) along with the inclusion of 
low level beam heating early in the flare, may have the capability to reproduce 
the pre-flare heating that is observed in the GOES lightcurve. For this work 
however, the behaviour of the flare was tracked as it cooled through the response 
functions of the many instruments and emission lines that were available. This 
provided details such as the ratio of the heating functions, which cannot easily 
be obtained from observations. It was found that the data were best reproduced 
when the plasma was heated approximately equally by direct and non-thermal 
mechanisms. This implies that both of these processes are vital during the flaring 
process and that flares may not be energised solely by non-thermal particles, as 



previously believed (e.g. Brown, 1971). This is in agreement with recent results 



found by Milligan (2008). There it was shown that a non-thermal electron beam 
is not necessarily required to obtain the high-temperature, high-density material 
we see in flares. However, it should be noted that the EBTEL value of the flux 
of non-thermal electrons required for equal heating is below the critical value 
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for explosive evaporation hypothesised by Fisher et al. (1985c). This can be 
explained by the over-simplified treatment of the non-thermal beam by EBTEL. 
This requires caution for a fiare of this nature, where it is evident that non-thermal 
particles play an important role. The model assumes that all non-thermal energy 
is used for evaporating plasma upward into the loop. However, it is well known 
that should the beam fiux be sufficient, it will drive plasma downwards into 



the dense chromosphere (AUred et al., 2005 Woods et al. . 2004). Despite the 
approximations made by EBTEL, such as the homogenous nature of the loop or 
the disregard for the location of energy deposition, the temperature and emission 
measure curves reproduce observations very well. It is computationally efficient, 
running a complete simulation in a matter of seconds. 



The cooling timescales modelled by Cargill ( 1993 , 1994 ) were calculated to be 
within a factor of 2 of EBTEL timescales. The critical time, was found to be 
~15:24 UT for both models, however the temperature T^,{t^) did not agree. The 
temperature at r* was found to be 8 MK for EBTEL and 12 MK for Cargill. This 
discrepancy is as a result of the way in which the cooling mechanisms are treated 
by the models. Since EBTEL takes account of the energy loss by conduction and 
by radiation at all times during the event, heat is removed more quickly than by 
the Cargill model, which considers only a single energy loss mechanism at any 
one time. This highlights the importance of considering the simultaneous cooling 
by both conduction and radiation. 

This chapter has established a method that will be applied to the analysis 
of future events. For this case, the data was manually compared to theoretical 
models. However, the fitting of the parameters together with model comparison 
techniques has been investigated using a Bayesian technique for simulating values 
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from the posterior distributions of the parameters (Adamakis et a/., 2009). The 
results of the statistically robust approach have been found to be within the er- 



rors of the parameters presented by Raftery et al. (2009). The main discrepancy 
is regarding the ratio of the thermal to non-thermal heating. The value found by 



Adamakis et al. (2009) is significantly larger than that found by the author (~ 23 



vs ~4). This can be explained since in the calculation of this parameter, Raftery 



et al. (2009) includes the background heating while [Adamakis et al. (2009) does 
not. However in each case it was found that Fthermai I Fnon-thermai > 1, mean- 
ing the thermal (direct) heating of the plasma must have at least as significant 
an impact on the heating of the plasma as beam heating. Statistically, the use 
of a Gaussian profile for the heating function was found to be in better agree- 
ment with observations than a half-Gaussian, which corroborates the assumption 
that the heating function can be obtained from the derivative of the soft X-ray 
flux. The Bayesian technique will be used when comparing theoretical models 
to future data sets. The authors intend to carry out an investigation of flare 
hydrodynamics using the improved cadence and extensive spectral range of the 
Extreme Ultraviolet Imaging Spectrometer (EIS) on board Hinode. Combining 
these data with RHESSI spectral fits will vastly improve observations and allow 
for even more accurate modelling. 
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Chapter 6 

The flare-CME connection 



The connection between solar flares and CMEs has long been a topic of debate. 
We investigate this relationship using a well observed flare-CME event that oc- 
curred on the east limb of the Sun on 31-December-2007. The event was analysed 
using X-ray observations from MESSENGER and RHESSI, EUV observations 
from STEREO B/EUVI and white light observations from STEREO B/Corl and 
Cor2. The hydrodynamic evolution of the flare was analysed along with the kine- 
matic evolution of both the flare and the CME. Evidence for pre-eruption magnetic 
reconnection was found to occur at least 2 minutes before the CME was launched. 
This was evidenced by a soft X-ray source at the top of the pre-flare arcade and 
temperature and emission measure values significantly higher than equilibrium. 
At the time of eruption, a HXR source (30-50 keV) was observed above the ther- 
mal sources and the RHESSI HXR flux reached a maximum. Evidence of breakout 
reconnection was not observed until 8 minutes later when a second set of HXR 
fluctuations were observed in conjunction with the restructuring of neighbouring 
magnetic field. Since magnetic reconnection occurred before the CME erupted and 
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breakout reconnection was not observed until minutes after, it is evident that this 
event was not triggered by an ideal MHD instability or by external tether- cutting 
mechanisms. The pre-flare SXR source and the HXR source observed at the 
time of the eruption were both found to be located above the central arcade. This 
provides evidence that this CME was in fact triggered by internal tether- cutting 
reconnection. This work is currently in review for publication in the Astrophysical 
Journal (^Raftery, Gallagher, McAteer, Lin & Delahunt, 2010, in review). Con- 
siderable intellectual conversation was conducted with Dr. Peter Gallagher, Dr. 
Chia-Hsien Lin and Dr. James McAteer. In addition. Dr. McAteer constructed 
Figure \6A and Gareth Delahunt conducted the SXR time-height analysis with the 
assistance of the author and Dr. Peter Gallagher. 
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6.1 Introduction 



The behaviour of CMEs and eruptive flares have been investigated on two fronts - 
theoretically and observationally. Unfortunately, the comparison between theory 
and observations is difficult and in most cases is a qualitative one. The mechanism 
involved in the initiation of CMEs is still a topic of hot debate with many com- 
peting theories, three of which are internal tether-cutting, external tether-cutting 



and ideal MHD instability. These models are discussed in detail in ^ 1.3.1 



In the case of internal tether-cutting (Figure 1.22), magnetic reconnection be- 
tween the legs of a sheared arcade (i.e. "internal" reconnection) trigger the erup- 
tion. Following this, breakout reconnection above the freed "plasmoid" removes 
flux from its path and allows it to erupt as a CME. For external tether-cutting 



(Figure 1.23) or breakout reconnection (Antiochos, 1998; Lynch et al, 2004), these 
processes are reversed. The eruption is triggered by breakout reconnection above 
the central arcade (i.e. "external" magnetic reconnection). As the central arcade 
rises, the magnetic fleld in the legs is stretched and drawn together beneath the 
reconnection region, generating a second current sheet. Internal tether-cutting 
along this second current sheet frees the plasmoid, enabling it to erupt upwards. 
The flnal scenario involves the CME being initiated by an ideal MHD process 



(Figure 1.24). Should the central arcade be forced out of equilibrium by e.g. flux 



emergence or continued shearing of the fleld, it rises upwards to seek a new equi- 
librium. As this happens, it creates current sheets between its legs and above it 
with the overlying fleld, as before. Magnetic reconnection then takes place along 
these current sheets, resulting in the eruption of an accelerated CME. 

Observationally, complete analyses of flare-CME systems have been hindered 
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by the lack of suitable data. Considering the broad range of temperatures 
(~8,000 K to ~20 MK), energies (few eV to MeV), and distances (1 to >30 
R0), it is clear that a multi-instrument approach is necessary. In the past, the 



kinematics of CMEs (e.g. velocity, acceleration) have been well studied. Gal- 



lagher et al. (2002a) presented the first observations of a rising SXR source that 



occurred in conjunction with a coronal mass ejection. The kinematics of the SXR 
loops were analysed in detail and the thermal emission was found to originate 
from successively higher altitudes as the flare progressed, agreeing with the stan- 
dard flaring picture. The CME associated with this flare was studied in detail 



by Gallagher et al. (2003). The CME's acceleration was found to be best rep- 



resented by a double exponential function. Temmer et al. (2008) furthered the 
investigation into the flare-CME connection and found that the CME accelera- 
tion occurs simultaneously with the hard X-ray (HXR) burst of the corresponding 
flares. This lends further support to the proposal by [Zhang et al. (2001 ) who sug- 
gested that CMEs and flares are driven by the same mechanism but do not have 
a cause and effect relationship. The importance of the flare-CME onset has been 



well documented by Harrison & Bewsher (2007). Using the evolution of the pre- 
flare arcades though a series of EUV spectroscopic observations, they established 
the importance of the pre-eruption activity. Often, the temperature and density 
(emission measure) evolution of a system is conducted using spectroscopic data. 



such as that used in Harrison & Bewsher (2007) and Raftery et al. (2009). How- 
ever, since these instruments are most effective close to disk centre, observations 
of limb flares are not as readily available. An alternative method of studying 
flare hydrodynamics is to use the spectroscopic capabilities of RHESSI. However, 
since RHESSI is primarily designed to study high energy emission, caution must 
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be taken when analysing the lower energy end of the spectrum, especially when 
the attenuators are in use. As such, while high temperature emission can be 
modelled accurately with RHESSI, lower temperature emission (<5 MK while 
in the Al state) is harder to observe. In comparison, CMEs originating on the 
limb are better observed than those originating on disk. For limb events, the 
components of the CME can be easily observed against the sky and avoid con- 
tamination by disk emissions. As a result of these obstacles, the hydrodynamic 
evolution of eruptive flares with well observed CMEs are rare. In this Chapter, 
the hydrodynamic evolution of a CME-associated solar flare is examined in a 
unique way, using MESSENGER/SAX. The temperature and emission measure 
of the flare is investigated in conjunction with the kinematic evolution of both 
the CME and the post-flare loop system. The availability of this unique data 
set has made this extensive study possible. The instruments used were located 
throughout the heliosphere and provided excellent coverage of the event, both 
temporally and spectrally. Along with MESSENGER/SAX the instruments used 
included RHESSI, GOES-12, EUVI, Corl and Cor2. The observations and data 



analysis techniques are discussed in ^6.2 ^6.3 describes the main results of this 



investigation which are discussed in light of current theory in ^6.4 



6.2 Observations and data analysis 

This event occurred on the east limb of the Sun on 31-December-2007. The 
SXR flux began to rise from approximately 00:30 UT. The CME was launched 
at 00:48 UT and the SXR flux was above background levels for more than 4 
hours. As observed from Earth, the footpoints in the low corona were occulted. 
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Figure 6.1: Position of the spacecraft used to observe the flare-CME event on 
2007 December 31. RHESSI and GOES are orbiting Earth. This is a top-down 
grid overlaid on a side-on image taken with STEREO B/EUVI (blue) and Corl 
(black and white). 

Therefore, Earth orbiting satellites (RHESSI and GOES) only observed looptop 
emission from the event. As a result, it is likely that that GOES classification 
of C8.3 is an underestimation of the total flux. STEREO B and MESSENGER 



however, had a clear view of the entire system, as Figure 6.1 shows 



6.2.1 X-ray spectroscopy 

The RHESSI Al attenuators were in place for the duration of this event. As 



such, the RHESSI spectrum was only analysed above 6 keV (see Table 3.1 and 



Figure 6.2). A sample spectrum taken at the peak of the flare is shown in Figure 
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Figure 6.2: Top left panel shows the MESSENGER/SAX spectrum from the peak 
of the flare between 00:47 UT and 00:52 UT (solid black line). The isothermal 
fits, shown as red asterisk (low temperature) and blue diamonds (high tempera- 
ture) were applied between 1-9 keV. Top right panel shows the RHESSI spectrum 
taken between 00:47 UT and 00:48 UT (solid black line) and its fit components 
constrained between 6-30 keV: an isothermal function (blue), a broken power-law 
(red) and a single emission line (green) . The normalised residuals are shown in the 
bottom panels. 



6.2 (right) showing the observed flux (black) along with the fit components. The 



thermal continuum was well modelled using a Maxwell-Boltzmann distribution. 



as discussed in ^2.3 (blue line) and the non-thermal emission was modelled us- 
ing a broken power-law (red line). A single emission line was also included at 
8 keV (green line). By fitting the thermal part of the RHESSI spectrum with a 
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Maxwell-Boltzmann distribution, the temperature and emission measure could be 
calculated during the early decay phase of the event. As the lower energy range 
of the spectrum (<6 keV) was not analysed (due to the attenuator state), any 
plasma with a temperature of < 5 MK present in the loop was not observed with 
this instrument. The MESSENGER/SAX instrument was used to account for any 
low temperature emission. This instrument is designed to measure characteristic 
X-ray emissions from the surface of Mercury. Incident solar flux is measured us- 
ing SAX for calibration purposes (when analysing Mercury). SAX is sensitive to 
the 1-10 keV range, overlapping with the thermal range of the RHESSl spectrum. 
However, significantly better spectral resolution and sensitivity to lower energies 
allows both high and low temperature emission to be observed and modelled ac- 



curately. Figure 6.2| (left) shows a sample SAX spectrum taken at the peak of 
the flare. The spectra were fit using the OSPEX software in the SolarSoft SPEX 
package. Two isothermal functions were used to approximate a differential emis- 
sion measure function to accommodate both high and low temperature plasma 
present in the flare. 

6.2.2 Imaging 

RHESSl images were reconstructed using the CLEAN algorithm and detectors 
3, 4, 5, and 6 integrated over 2 minute periods. The source was imaged in the 
3-6, 6-12, and 12-25 keV energy bands. The 30-50 keV coronal sources imaged by 



Krucker et al. (2009) have also been incorporated into this study. The evolution 
of the flare and the acceleration phase of the CME were imaged using 171 A 
EUVI passband on board STEREO B, while the propagation phase of the CME 
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was imaged using Corl and Cor2. 



Figure [673] shows the evolution of the event. The left hand panels are interlaced 
difference images from EUVI and Corl (a - d) and Cor2 (e). Overlayed on these 
are the RHESSI sources, 3-6 keV (red), 6-12 keV (yellow), and 30-50 keV (purple, 
panels b and c only). The flaring region of the event has been expanded on the 
right hand side to highlight the evolution of a post-fiare loop system and the 
motion of the soft X-ray sources. 

6.3 Results 

This CME-flare system was observed for more than four hours using a multitude 



of instruments. Figure 6.3 shows the evolution of the system and Figure 6.4 
explicitly shows the motion of both the CME and SXR loop top sources (inset). 
As RHESSI emerged from eclipse at ~00:46 UT, a SXR source was observed to lie 
between and above two bright EUV ribbons (Figure [673^ ). Between 00:46 UT and 
00:48 UT, the SXR sources appeared to move in a southerly direction, decreasing 



in altitude at the same time (Figure 6.4 inset). The decreasing altitude of the 



SXR source may correspond to the collapse of a magnetic X-point in X- to Y- 



type reconnection (Sui fc Holman, 2003). At 00:48 UT (Figure 6.3d), a 30- 



50 keV coronal source was observed above the now extended thermal sources. 
Simultaneously, the CME was launched. Figure 6.3| 3 also shows the location of 



the SXR source to be at the base of the CME, between the flare footpoints. 
Following the CME eruption, the SXR source began to rise beneath the CME 
while still moving south along what would become the post-flare arcade. This 



can be seen in the main frame of Figure 6.4: the direction of motion of the post 
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Figure 6.3: Evolution of the 31-Dec-2007 flare-CME as observed by EUVI, Corl, 
Cor2 (panel e only), RHESSI 3-6 keV (yellow), 6-12 keV (red), and 30-50 keV 
(purple). 
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Figure 6.4: Position of the CME front in the Corl field of view (black crosses) 
along with the centroid of the RHESSI sources: 3-6 (blue diamonds), 6-12 (red 
triangles), 12-25 (green squares), and 30-50 keV (black asterisk). The inset panel 
is a magnified view of the flaring region. Time stamps for both the flare and CME 
are also marked. 



flare loops follow the same path as the CME. Figure 6^:: shows that the fleld 
surrounding the erupting arcade has been disturbed by the launch of the CME, 
primarily to the south. This disturbance of the neighbouring fleld is also evident 



in Figure 6.3 i, along with the 3 part structure of the CME - the core, the cavity 
and the front. The evolving post-flare arcade can also been seen in this flgure. 



The bottom panel of Figure 6^ shows the system approximately 2 hours after 
eruption. The CME is observed in Cor2, the post-flare loop system has evolved 
and the X-ray sources are still located above the post-flare EUV loops, at the 
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base of the CME. 



Figure |6.5| shows the kinematics of the flare and the CME, along with the 
hydro dynamic evolution of the flare. The CME liftoff time (00:48 UT; defined as 
the time the CME is first observed in EUVl, marked in Figure 6.3)p ) is highlighted 



in each panel with a dashed vertical line. Figure 6.5 1 shows the displacement of 
the RHESSl SXR sources (3-6, 6-12, and 12-25 keV) relative to Sun centre. The 
30-50 keV sources are also shown here. The average velocity of the loops as they 
descend (pre 00:50 UT) is 18 ±4 km s~^. The height of the ascending loops (post 
00:50 UT) were fit with a constant acceleration model: 



h{t) = ho + vot + lat'^ (6.1) 



where /iq, fo, and a, averaged across the three energy ranges are 711.1 ± 0.5 Mm, 



15 ± 2 km s ^, and —1.1 ± 0.9 m s ^. These fits are shown in Figure 6.5 1. 



Figure |63p shows the displacement of the CME apex relative to Sun centre 
measured using STEREO B instruments (EUVl, Corl and Cor2). A detailed 
analysis and comparison to theoretical models of the CME velocity and accelera- 



tion can be found in Lin et al. (2009). Panel (c) of Figure 6.5 shows the RHESSl 
25-50 keV HXR lightcurve (green). The peak in the hard X-rays corresponds 
exactly to the time the CME is launched. Overplotted on this is the GOES 1-8 A 



lightcurv^ Figure 6^i shows the thermal evolution of the flare using results 



from RHESSl, SAX, and GOES. The GOES temperature (and later, emission 
measure) were obtained from the ratio the two flux channels (see ^3.4|). It is clear 



^The appearance of multiple peaks in the GOES lightcurve is as a result of occultation and 
not multiple events. This was verified by comparing EIT intensity (partially occulted) to that 
of EUVl on STEREO B (not occuhed). 
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Figure 6.5: Evolution of the flare and CME kinematics and flare hydrodynamics 
in time, (a) SXR displacement from Sun centre, together with the fits to the 
constant acceleration model of the 3-6 (blue diamonds), 6-12 (red triangles), and 
12-25 keV (green squares) sources. The displacement of the HXR 30-50 keV source 
is also shown (black asterisk), (b) CME displacement from Sun centre in the EUVI 
(black asterisk), Corl (red diamonds), and Cor2 (blue diamonds) fields of view, (c) 
GOES 1-8 A soft X-ray lightcurve (black dashed line) overplotted on the RHESSI 
25-50 keV lightcurve (green solid line), (d) Flare temperature measured by GOES 
(black solid line), RHESSI (green triangles), and the two components of the SAX 
fit: high temperature (blue diamonds) and low temperature (red asterisk), (e) 
Flare emission measure from GOES (black solid line), RHESSI (green triangles), 
and the two components of the SAX fit: high temperature (blue diamonds) and 
low temperature (red asterisk). 
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that the lower temperature component of the SAX spectrum remains approxi- 
mately constant during the initial and early decay phases of the flare while the 
higher temperature component rises quickly to a maximum value of 19±2 MK. 
The high and low temperature components converge as the flare decays. Figure 



6.53 shows the evolution of emission measure of the flare, again from RHESSI, 



SAX, and GOES. As expected, the level of high temperature emission increases 
as the flare approaches its peak. Surprisingly however, the level of cooler emission 
also increases during the flare impulsive phase. 



6.4 Conclusions and discussion 

This chapter presents observations of a CME and associated solar flare that were 
observed using a range of spacecraft positioned at various points in the helio- 
sphere. The eruption took place on the east limb of the Sun and the footpoints 
were occulted as viewed from Earth. Before the eruption, a SXR source was ob- 
served in conjunction with increasing flare temperature. At 00:48 UT, the time 
the CME was launched, a 30-50 keV HXR looptop source was observed. Follow- 
ing the liftoff of the CME, the SXR source was found to rise beneath it as the 
temperature and emission measure both began to fall. 

The pre-eruption increase in temperature and emission measure, along with 
the presence of the SXR source early in the event is strong evidence for pre- 



eruption reconnection. Following the interpretation of Sui & Holman (2003), we 



also believe that this event exhibits signs of X- to Y-type reconnection. |Sui Ez 



Holman (2003) observed a looptop source that appeared to shrink during the 



impulsive phase of a flare (Figure 6.6). Following the flare peak, the loops were 



170 



6.4 Conclusions and discussion 




Figure 6.6: Observations published in Sui &: Holman (2003) that display similar 



behaviour to the observations presented in this chapter. Top panel: RHESSI light 
curves in three energy bands (from top to bottom): 3-12, 12-25, and 25-50 keV, 
scaled by 2.0, 0.5 and 1.0 respectively. Middle panel: Time histories of the loop 
height (obtained from 10-12 keV images) and the coronal source height (obtained 
from 10-25 keV images). Bottom panel: Height of the loop and the coronal source 
at different energies at 23:11:00 UT. The horizontal bars represent the energy band- 
widths of the RHESSI images. 

observed to rise. In addition, they too observed a HXR coronal source that 



was located above the looptop sources (Figure 6.7). They hypothesise that a 



current sheet is located between the looptop sources (on disk in Figure 6.7) and 
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Figure 6.7: RHESSI images in different energy bands from 
The three contours on the solar disk indicate the top of loops in the energy bands 
(contour line shade from light to dark): 6-8, 10-12, and 16-20 keV. The contours 
above the limb are (from light to dark) 10-12, 12-14, and 14-16 keV bands. The 
asterisk marks the centroid of each source. The crosses mark the two footpoints of 
the X-ray loop. 

the coronal source (off limb). In this chapter, the coronal source was observed 
for only two instances in time so a direct comparison is difficult. However, the 
apparent downward motion of our looptop (SXR) source is in very good agreement 



with the motion of that of Sui & Holman (2003). They explain the stratification 



in temperature (energy sources) within the flare loops as the heating of material 
close to a current sheet formed between the looptop and the coronal source. 



According to Krucker et al. (2009), this coronal source is the location of particle 



acceleration in the Dec-31 event. This, therefore is a good explanation for why 

inset 



the 12-25 keV source is located above the 3-6 keV source (see Figure 6.4 



for clarification). The temperature distribution can be explained in terms of the 
energy supply by reconnection at the coronal source heating the plasma closest to 
the current sheet. This produces the higher energy sources. Subsequent cooling 
of the inner loops result in the lower energy sources observed at lower altitudes. 
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Figure 6.8: Schematic diagram showing the asymmetric eruption of a CME that 



fits well with the observations presented in this chapter. (Tripathi et al, 2006) 



In addition, the loop shrinkage is believed to be due to a change in the magnetic 
field configuration as an X-point collapses into a current sheet. As this happens, 
the looptop source moves downwards, as is observed in this chapter. Since the 
observations presented here are clearly in good agreement with |Sui &: Holman 



(2003), we can conclude that we are observing the changing in magnetic topology 



from an X-point (X-type) to current sheet (Y-type) reconnection region. 

The extension of the SXR source across the arcade top may be interpreted 
as successive reconnection events progressing along the arcade. These reconnec- 
tions are likely be the source of the heated flare plasma, either by chromospheric 



evaporation (Brown, 1971 Milligan et a/., 2006a Raftery et a/., 2009) or by a 



collapsing magnetic trap (Joshi et al, 2008 Karlicky & Kosugi, 2004 Somov & 



Kosugi, 1997; Veronig et al, 2006a). Pre-eruption magnetic reconnection would 



also explain the apparent motion and elongation of the SXR source: as the event 
progressed and neighbouring loops within the arcade reconnect, emission from 
more and more loops contribute to the SXR source, resulting in the apparent 
elongation of the source along the arcade. This is similar to the asymmetric 



eruption described in Tripathi et al. (2006) and shown in Figure 6.8 
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The timing in this event is crucial for the interpretation of the initiation 
mechanism. The three indicators are the initial rising of the CME, the evidence 
for magnetic reconnection in the central arcade (e.g. rising temperature, X-ray 
sources) and evidence for reconnection in neighbouring arcades. This event began 
at~00:46 UT with the appearance of a SXR source at a height of 22±1 Mm 
above the photosphere. This corresponds very closely to the height of the pre- 
emption arcade loops. At this time, the temperature of the region has risen from 
its equilibrium value of 5.3±2 MK to 17±2 MK (see Figure |6^). The rise m 
temperature and the SXR source are clear evidence that magnetic reconnection 
has taken place before the CME is launched. This rules out the ideal MHD 
trigger as a mechanism for CME initiation. At ~00:48 UT, approximately 2 
minutes later, we see the first evidence of a rising CME in EUV difference images 



(see Figure 6.3 d) with the SXR source at the base of the CME. It is not until 
~00:56 UT, that any restructuring of side lobes is observed. This would indicate 
that the breakout reconnection occurred after the internal reconnection between 
the legs of the arcade. It is possible that the limitations of the EUVI's 171 A 
passband did not observe the neighbouring loops until after they were heated 
and were cooling back into the 171 A passband. However, observations of the 
284 A passband, whose peak sensitivity is at 2 MK did not observe any change 
until this time either. As well as this, RHESSI observed an increase in HXR 
fluctuations from 00:57 UT. This can be interpreted as evidence that breakout 
reconnection occurs after the CME is initiated. Therefore, we conclude that this 
event is most likely triggered by an internal tether-cutting mechanism. This is 



in good agreement with Chifor et al. (2007) who also diagnosed tether-cutting 
CME initiation from the presence of HXR sources (or "precursors") during the 
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early phase of an eruption and observed a similar asymmetric evolution to their 
eruption. 

The unique multi-spacecraft observations used for this study, specifically the 
high spectral resolution of SAX, enabled a detailed analysis of the hydrody- 
namic evolution of an eruptive flare. High resolution data from the STEREO 
and RHESSI instruments allowed the kinematic evolution of both the flare and 
the CME to be analysed simultaneously. Our results clearly demonstrate there is 
a physical connection between the two phenomena. We can also disregard a num- 
ber of theoretical models that are not applicable to this event and conclude that 
this CME was initiated by an internal tether-cutting mechanism, with breakout 
reconnection occurring afterwards. 
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Conclusions and Future work 



The goal of my thesis work was to better understand the behaviour of the active 
Sun. This incorporated a study of a confined solar fiare, a study of an eruptive 
flare/CME system and an investigation into the sensitivity of EUV imaging tele- 
scopes. In this chapter, the primary findings of these investigations and their 
implications are presented. This is followed by a discussion of the problems en- 
countered in this thesis and how they they can be rectified in the future. 
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7.1 Thesis results and comments 



The goal of my thesis was to improve our understanding of the active Sun, specif- 
ically the behaviour of solar flares. It was found that the use of multi-wavelength 
analysis is essential when studying the hydrodynamic properties of solar flares. 
The broad temperature range that these events cover require at least observations 
from HXR through SXR and into the EUV regime. Ideally this should be ex- 



tended even further to include UV, white light and radio emissions (e.g. Alexander 



& Coyner, 2006 


Benz & Pianezzi 11997; 


Benz et al. 


2009; Coyner &; Alexander 



2009 Fletcher et al. . 2007). Attempting to understand the full effects of a solar 
flare using a single passband is not advisable. The temperature of a flare can 
change very quickly, especially during the impulsive phase and early decay phase 



(Raftery et al., 2009, 2010b). Therefore a single passband, sensitive to a narrow 



range of temperatures will only observe the flare for a short period of time. A 



clear example of this is shown in Chapter [sj the 171 A images in Figure 5.1 show 
that the TRACE instrument did not register the occurrence of a solar flare until 
approximately 30 minutes after it began. Even the use of multiple broadband 
telescopes for such studies must be done with care. For example, [Aschwanden 



& Alexander (2001) utilised integrated flux measurements from four imagers - 
TRACE, Yohkoh/HXT, SXT and GOES - to analyse the cooling curve of the 
Bastille day flare. The cooling curve was obtained by assigning the characteristic 
temperature of a particular instrument to the time of maximum flux, much like 



the approach taken by Raftery et al. (2009). Unlike Raftery et al. (2009) how 



ever, integrating the flux of a source observed with a broadband imager can lead 
to problems in obtaining the temperature of the emission, especially during the 



177 



7.1 Thesis results and comments 



impulsive phase of a flare. As Chapter |4] ( |Raftery et al. 2010a) clearly demon- 
strates, the temperature responses of broadband instruments vary significantly 
under atypical conditions, such as a flare. Therefore, the nominal temperature 
assigned to these instruments is not likely to be reliable during the impulsive 
phase and early decay phase of a flare. 

The question of plasma temperature can be alleviated through the use of spec- 
troscopy. Emission line spectroscopy, specifically the use of rasters, is extremely 



useful for investigations of solar flares (e.g. Chapter pi Brosius, 2003 Cirtain 



eTd] [20071 iMiUigan et olj |2006a||bl [Raftery et alj [20091 |Schmelz et al.\ [20071 ). 



The ability to image a solar flare simultaneously in multiple temperatures can 
provide a detailed insight into the evolution of the event along with the distribu- 
tion of plasma temperature and densities in the system. It can also be used for 
the analysis of plasma velocity through the use of Doppler shifts. In Chapter [5] of 
this thesis, we utilised the presence of upflowing Fexix plasma at 95 ± 16 km s~^ 
to diagnose pre-flare heating. A combination of hot (Fe xvi and Fe xix) upflows 
at up to 81 ± 16 km s~^ and cool (He I and Ov) downflows at approximately 
30 km s~^ were characteristic of explosive chromospheric evaporation. Other 
studies have utilised the versatility of the multi-temperature nature of spectro- 



scopic observations. For example, Milligan & Dennis (2009) used the broad range 
of emission lines observed by the EUV Imaging Spectrometer (EIS) on board Hin- 
ode to flnd the cut off temperature at which plasma is evaporated upwards into 
the loop rather than downwards towards the chromosphere. This was achieved by 
studying the Doppler shifts of 15 Fe lines from 0.05 - 16 MK. A clear cut-off was 
observed Fexiv (2 MK) whereby plasma at temperatures above that ion were 
evaporated upwards and plasma at temperatures below flowed downwards. 
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This thesis also showed that soft X-ray spectroscopy is also a very powerful 
tool when analysing the hydrodynamic properties of plasma (Chapter |6, Raftery 



et al, 2010b). The sensitivity of thermal continuum to both temperature and 



density makes it ideal for analysing their variation in time. A comparison of 
the RHESSl SXR continuum to that of SAX revealed that while RHESSI is not 
as reliable when analysing plasma at temperatures of less than ~5 MK (in the 
Al state), it is ideal for observing emission at high temperatures that cannot be 
observed with e.g. SAX. In Chapter [6] ( [Raftery et al. 2010b), the combination 
of these two instruments enabled the author to accurately analyse a solar flare 
across its entire evolution, from pre-fiare emission at <5 MK using SAX to 19 MK 
using RHESSl. The limitations of the SAX spectrum, sensitive to emission up to 
~ 9 keV constrained the upper fitted temperature to ~17 MK. Therefore, careful 
analysis of both instruments provides the most complete and accurate results. 



7.1.1 Principal results 

The primary aim of this thesis was to better understand the behaviour of solar 
flares. This was done by examining the differences and similarities between an 
eruptive and a conflned flare through quantitative and qualitative comparison of 
multi-wavelength observations to theoretical models. From this comparison, we 
draw the following conclusions regarding the behaviour of flares. 

• In the case of the conflned flare (CP) the HXR proflle consisted of a single 
burst during the impulsive phase that lasted approximately 3 minutes. The 
eruptive flare (EF) however was observed to have a single HXR burst during 
the impulsive phase, lasting 4 ± 2 minutes, followed by HXR bursts of lower 
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magnitude but similar duration during the decay phase. This difference 
is most likely due to the fact that the CF does not experience continuous 
reconnection to any great extent during the decay phase. However, as a 
CME erupts above an EF, it is believed that breakout reconnection takes 



place to free the CME (Antiochos, 1998 Moore & Sterhng 2006 Raftery 



et al. 


2010b; 


Zhang et al. 


2001 



cause of the continued HXR bursts. 

The continued HXR flux that occurs in conjunction with the hypothesised 
breakout reconnection, mentioned above, corresponds in time to the re- 
structuring of the magnetic field surrounding the erupting arcade. This 



manifests itself as an "EIT wave" (Thompson et al., 1998). The nature of 



these disturbances has been debated for many years, with the interpretation 



varying between various types of waves (e.g. Long et al., 2008; Patsourakos 



et al, 2009) to the chromospheric component of breakout reconnection (e.g. 



Attrill et al. , 2007). The association between the HXR flux and this "wave" 
lends support to the CME theories: that the wave is in fact observations 
of the coronal magnetic field being restructured to allow the release of the 
CME. 

The time to heat the flare to its maximum temperature took a comparable 
amount of time for both events. The CF took approximately 11 minutes to 
reach its maximum temperature of 15±2 MK while the EF took approxi- 
mately 20 minutes to reach a temperature of 19±2 MK. 

The time it took for the flares to cool varied considerably. In the case 
of the CF, the SXR flux took approximately 1 hour to reach background 
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levels following the peak. The EF flare however, took more than 7 hours to 
return to equilibrium. This could also be attributed to the extended period 
of reconnection hypothesised above, leading to the continued heating of the 
plasma within the flare for a longer time, thus keeping temperatures higher 
for longer than in the conflned case. 

• The topological conflgurations of each system were considerably different. 
The CF existed as an individual loop with discrete footpoints. The flare was 
an isolated event that did not affect the global magnetic fleld to any great 
extent. The EF however occurred in extended arcade whose footpoints were 
extended ribbons. It was a global event that occurred in conjunction with 
a CME eruption and an EIT wave that reconfigured the magnetic field on 
a scale much larger than the fiare. 

The comparisons and differences between the eruptive and confined events has 
led the author to the conclusion that while eruptive and confined fiares have very 
different properties and manifest themselves in different ways, the fundamental 
processes involved are the same in both cases. 

The analysis of the individual events had some interesting conclusions that 
are important in their own right. The main results and conclusions from each of 



the events studied in Chapters [5] and |6] are discussed in detail in §7.1. 2| and |7.1.3 
The use of EUV imagers was key to both the solar flare investigations and so a 
study into the sensitivity of EUV imagers in the 171 A passband was carried out 
to ensure accurate interpretation of EUV images in Chapter |4j The results and 



conclusions of this study are presented in ^ 7.1.4 
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7.1.2 Hydrodynamic modelling 

1. The use of a 0-D hydrodynamic model, EBTEL, proved to be extremely 
useful for the interpretation of observations. Despite averaging tempera- 
ture, density and pressure through the corona, it reproduced observations 
accurately and in a very efficient manner. Combining modelling results with 
observations revealed flare properties that would otherwise be unattainable 
such as the significance of various heating and cooling mechanisms. 

2. The approximate magnitude and duration of the HXR burst were deter- 
mined from this comparison and found to be 5 x lO*^®^^-* ergs cm~^ s~^ 
and 100 ± 50 seconds respectively. These are not exact values as the over- 
simplified consideration of the HXR heat fiux by EBTEL will affect the 
total flux. In order to accurately determine the total non-thermal flux, a 
more thorough treatment of the effects of beam heating must be considered. 
The non-thermal flux determined by EBTEL is most likely underestimat- 
ing of the actual non-thermal flux since it is approximately an order of 



magnitude less than the flux required by Fisher et al. (1984) for explosive 
chromospheric evaporation. This is an important result as it will effect how 
EBTEL should be treated in future investigations. 

3. Despite the issue of the treatment of the non-thermal heat flux, EBTEL 
enabled the author to get a reasonable approximation for the non-thermal 
heat flux and for the direct heating rate. Although the values obtained from 
the model were not absolute, it was however, possible to establish that solar 



flares are not heated purely by beam heating as hypothesised by e.g. Brown 



(1971). They also require some kind of direct heating mechanism such as 
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shock heating or particle trapping at the top of the loop (e.g. Joshi et al. 



2008 Karlicky & Kosugi, 2004; Veronig et al. 2006b) 



4. An investigation into the significance of fiare cooling mechanisms proved 
the importance of simultaneous cooling by both conduction and radiation 
throughout a solar fiare. Conduction can be considered independently at 
the beginning of a fiare decay phase when temperatures are high and con- 
duction is efficient. Likewise, radiation can be considered independently 
late in the decay phase when temperatures are low and conduction is not 
efficient. However, for the majority of the fiare cooling phase, both con- 
duction and radiation are important mechanisms for the removal of heat 



from the system. Therefore, following the techniques of e.g. Aschwanden 



& Alexander 


( 


2001 


); 


Cargill 


( 


1994 


); 



(1994), using a single loss mechanism at any one time will significantly 
reduce the efficiency of heat removal from the system, resulting in longer 
cooling timescales and higher predicted temperatures, as demonstrated in 
Chapter |5| 

5. Studying the CDS spectral data and the associated Doppler shifts revealed 



evidence for pre-fiare heating and chromospheric evaporation (Raftery et al. 



2009 ) . Analysis of the Fe xix images revealed the loop emitting at 8 MK up 
to 10 minutes before the peak of the hard X-rays. Analysis of the Doppler 
shifts of that line resulted in the fastest velocities for the entire event, at 95 ± 
16 km s~^. Following the precursor phase, explosive chromospheric evapo- 
ration was determined from a non-thermal fiux of 5 x lO*^^^^-* ergs cm~^ s~^ 
combined with 8 MK plasma rising into the loop at 81 ± 16 km s~^ and 
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0.25 MK plasma being forced towards the chromosphere at 33 ± 16 km s ^. 
These observations (with the exception of the low non-thermal flux, see 



above) are characteristic of explosive chromospheric evaporation (Fisher 



et ai, 1984 Milligan et ai, 2006a). The explosive evaporation is believed 



to have been followed by a period of conduction driven gentle chromospheric 



evaporation (Fisher et al. 1984 Milligan et al, 2006b). This was evidenced 



by plasma rising into the loop at much lower velocities: 24 ± 16 km s~^ for 
Fe XIX and 9±16 km s^^ for He I. This is believed to be driven by conduction 
as the flux of non-thermal electrons had almost returned to background lev- 
els by this time and, according to EBTEL simulations, conduction is highly 
efficient. Thus, the most likely explanation is conduction fronts from the 
hot (~20 MK) loop top propagate towards the chromosphere where they 
thermalise the plasma at a much slower rate than a non-thermal flux. The 
small pressure increase due to heating results in a minimal deviation from 
equilibrium. Thus plasma is forced upwards into the loop. The pressure 
is not sufficient to overcome the large densities through the chromosphere 
and so no downfiows are observed. 

6. Despite the lack of a rigorous statistical test when comparing the results 
of the hydrodynamic model to observations, the input parameters used by 
the author in the EBTEL modelling of a confined fiare and published in 



Raftery et al. (2009) agree very well with those found by Adamakis et al. 



(2009) using Bayesian statistics to optimise the model parameters. 
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7.1.3 CME initiation mechanism 

1. The initiation mechanism for the launch of a CME was investigated by qual- 
itative comparison of multi-wavelength observations and theoretical models. 
SXR loops are believed to be a signature of magnetic reconnection (^1.2| 



Dennis & Schwartz, 1989). Therefore, the presence of an SXR source at the 
top of the arcade before the CME erupts suggests magnetic reconnection 
took place within the arcade system (i.e. not breakout reconnection). 

2. The presence of an X-ray source observed before the CME began to rise 
excludes triggering by an ideal MHD instability. According to the catastro- 
phe model, a CME fiux-rope forced out of equilibrium should rise to find 
a new equilibrium. This is believed to take place initially through an ideal 



MHD process i.e. without magnetic reconnection (e.g. Forbes & Isenberg 



1991 Forbes & Priest, 1995 Isenberg et ai, 1993). Therefore, the possible 



evidence for magnetic reconnection observed before the loops began to rise 
(mentioned above) suggests that the catastrophe model is not applicable to 
this event. 

3. Approximately 8 minutes after the main CME lift off, restructuring of the 
surrounding magnetic field was found to correspond to secondary HXR 
fluctuations that occurred after the main impulsive burst at the time of 
launch. This may be representative of breakout reconnection occurring after 
lift-off. Breakout is expected to occur high above the photosphere and so the 
location is not discernible from current instruments. However, considering 
the chromospheric field motions along with secondary HXR peaks imply 
some large-scale field reorganisation which fits with the expectations of 
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Antiochos (1998). However, unlike the predictions of the breakout model, 



the CME is not triggered by breakout reconnections, rather they occur in 
order to allow for the propagation of the CME. 

4. Considering the observational evidence and a qualitative comparison to the 



models presented in Moore & Sterling (2006), we can conclude that the 



CME in question was triggered by internal tether-cutting reconnection, fol- 



lowed by breakout reconnection after the CME had been launched (Raftery 



et ai, 2010b). This is an interesting result as evidence for internal tether 



cutting is reasonably scarce (e.g. Chifor et al. . 2007 Sterling & Moore 



2001), compared to observations and simulation of the other two models 



see e.g. Chen & Shibata, 2000 Lynch et al. 2004 Roussev et al. 2003). 



5. The behaviour of the eruptive flare was found to be in good agreement with 



the CSHKP model (Carmichael 1964; Hirayama, 1974; Kopp & Pneuman 



1976 Sturrock 1966). The positioning of the HXR source above the pre- 



existing SXR source at the time of eruption agrees with the general topology 
of the CSHKP model, as does the motion of the SXR source between the 
footpoints of the CME. In addition, the connection between the SXR source 
and the CME topology and the relationship between the launch of the CME 
and the hydrodynamic behaviour of the flare reaffirms the postulation of 



Zhang et al. (2001): flares and CMEs are separate manifestations of the 



same driving mechanism, and therefore are not independent of each other. 
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7.1.4 EUV temperature response functions 

1. A new method for investigating the temperature response of EUV im- 
agers has been developed and presented in this thesis. This method ac- 
counts for the varying emission measure with temperature for four types 
of solar plasma: coronal holes, quiet sun, active region and flares. The 
method not only enables the interpretation of different plasmas, it can 
be applied to any EUV imager whose wavelength response is in units of 
DN phot~^ cm^ str pix~^. This will facilitate the comparison of results 
from different instruments. So far it has been applied to six instruments in 
the 171 A passband, however it can easily be adapted to accommodate any 
of the typical EUV passbands commonly used in solar physics. 

2. The interpretation of images taken with EUV imagers requires special at- 
tention. The typical response curves for the 171 A passband revealed stark 
differences between the response of EUV imagers to quiet sun/active re- 



gions and the extreme conditions of flares and coronal holes (Raftery et al. 



2010a). The maximum sensitivity of a 171 A passband imager when observ- 
ing coronal holes is 7.5±0.3 x 10^ K. This is cooler than the typical peak 
temperature of 10^ K. During a flare, the instrument will observe not only 
the "characteristic" temperature peak at 1 MK but a second peak above 
6 MK. 

3. The temperature response of 6 EUV imagers in the 171 A passband were 
studied and found to be in very good agreement. There was an offset 
between the instruments' responses, primarily as a result of their varying 
aperture size. However the shape of the AlA wavelength response was dra- 
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matically different to all the others and this had a large impact on the overall 
sensitivity of the instrument. It was found to be almost 300 times more sen- 
sitive than the next most sensitive instrument, EUVI on STEREO A. This 
suggests that developing instruments to study the spectrum around 171 A 
wavelength will produce similar results despite different optical systems or 
detector technologies. 



7.2 Future work 

This thesis incorporates many basic assumptions adopted in the solar physics 
community. For example, coronal abundances are assumed throughout this thesis. 
However, it has been shown that the ion abundance in the corona changes during 



non-relativistic electron events such as solar flares (Ramaty et a/., 1980 Reames 



et a/., 1985). The assumption that abundances remain constant throughout an 
energetic event is into question. While it is likely that this assumption is valid 
through the decay phase of flares, changing abundances during the impulsive 
phase will affect the hydrodynamic properties derived from observations and used 
to diagnose energetic phenomena. In order to fully understand the properties of 
solar plasmas, this effect should be taken into account. 

A second approximation is that of ionization equilibrium. For example, the 
hydrodynamic models used in Chapter [5] assumes that the heating rates involved 
in the solar flare are longer than those of the atomic processes. Again, since 



Equation 2.14 is based upon the assumption of ionization equilibrium and is used 
to calculate the emission measure and temperature of emission lines during solar 
flares, this calls into question the validity of our results. A very large quantity 
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of energy (~10 J; Emslie et al, 2005) is deposited in the chromosphere over a 
very short period of time (~minutes) and so it is quite possible that the heating 
timescales are at least comparable to the ionization timescales (~100 seconds; 



Mariska, 1993). In addition to this, when modelling the thermal continuum of 



RHESSI and MESSENGER, the assumption is made that the thermal electrons 
adhere to a Boltzmann distribution of velocities. This assumption must be made 
for RHESSI spectroscopy as there is no other way to differentiate the thermal from 
non-thermal components of the spectrum. Should there be a second population 
of electrons in, for example a bump-on-tail instability, which is known to occur 
during the production of type-Ill radio bursts frequently associated with solar 



flares (Benz & Saint-Hilaire , 2003), this second population would not be modelled 
accurately using a Maxwell-Boltzmann distribution. This would alter both the 
incoming flux of electrons and the temperature of the system. 

This thesis tested the robustness and applicability of the EBTEL hydrody- 
namic model. The model proved to be highly efficient in producing results and 
is an ideal candidate for use with statistical methods. It will be very useful for 
probing the broad parameter space such as the heating functions of solar flares. 
By narrowing the parameter space, this 0-D model can be used as a "flrst esti- 
mate" for more robust 1-D hydrodynamic models, such as those by Bradshaw fc 



Cargill (2005), Warren & Winebarger (2007) and Antiochos & Klimchuk (1991). 



In addition, the EBTEL model will be very useful in the investigation of the 
multi-thread model of solar flares. This theory involves the magnetic reconnec- 
tion of individual fleld lines within a larger loop (like strands within a rope). 
The number of fleld lines reconnected at a given time will combine to give you 



the time dependent heating of the loop (Klimchuk &; Cargill 2001 Klimchuk 
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et al. 2006; Patsourakos & Klimchuk, 2006). The short computation time of the 



EBTEL model will enable investigators to probe the frequency and magnitude of 
individual reconnection events, thus furthering the understanding of solar flares. 

The instruments available for analysis of the Sun have undergone significant 
improvement since the days of Skylab. The number of available spacecraft is at 
an all time high and the quality of the instruments by far surpasses anything that 
has gone before. For example, the innovative approach to stereoscopic imaging 
of energetic events used by STEREO has revolutionised our understanding of 
coronal activity, from the structure of coronal loops (Aschwanden et al, 2008a|[b 



2009) to understanding the propagation of CMEs (Byrne et a/., 2009 Maloney 



et al. 2009) 



The launch of Proba-2 has brought the commissioning of a new EUV imaging 
telescope, SWAP. The extra-wide field of view of this instrument will hopefully 
facilitate a more accurate analysis of the acceleration phase of CMEs. CMEs 
are believed to be accelerated close to the surface of the Sun, in the vicinity of 



the edge of EUV images (Chapter^ Gallagher et al. 2002a, 2003 Lin et al 



2009 Raftery et al. 2010b). SWAP however, will allow users to study evolving 



CMEs out to 3 Kq which will cover the entire acceleration and early propagation 
phase. In addition to this, the upcoming launch of SDO will inundate the solar 
community with more data than ever before - terabytes per day! The extremely 
high cadence and resolution of its instruments will give investigators a completely 
new insight into the dynamic nature of the corona. 

The CDS instrument has recently been accompanied by the Extreme ultra- 



violet Imaging Spectrometer (EIS; Culhane et al. 2007) on board the Hinode 
satellite. This instrument has higher cadence, spatial and spectral resolution 
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than CDS which has aheady advanced the understanding of coronal physics sig- 



nificantly in its three years since launch (e.g. Jin et al. , 2009 Mcintosh & De 



Pontieu 2009 Milligan k Dennis 2009). 



Despite these significant developments both in theory and observation, there 
remain many unanswered questions in the field of solar physics. With the pro- 
posed commissioning of new and innovative instruments such as Solar Probe, 
some of these questions may be answered. However, there remains significant 
work to be done. 



7.3 Outstanding questions 

In the past three decades, there have been many developments in the area of solar 



physics, specifically in the field of coronal disturbances. Figure |7.1| describes 
this idea well. Magnetic reconnection is believed to occur within the diffusion 
region (DR) below an erupting prominence (EP). Slow mode standing shocks 
(SMSS) occur at the boundary of outflow jets which may generate a fast mode 
standing shock (FMSS) if the downflow jet is supersonic. As the flux rope evolves 
into a CME, a heliospheric interplanetary (IP) shock can be driven infront of it, 
producing metric type-II bursts. Finally, large scale restructuring of the magnetic 
fleld can be observed to propagate radially from the eruption site. Despite this 
apparently concise picture, it is still unclear how these various phenomena are 
connected and often they are studied as independent events. Here I discuss some 
outstanding questions regarding the overall future direction of flare physics. 
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Figure 7.1: Schematic diagram from Warmuth (2007) describing the various coro- 



nal disturbances that are known to occur during a solar eruption (adapted from 



Aurass et al. 2002). 



7.3.1 What is the role of magnetic reconnection? 

It has become clear over the past few decades that magnetic reconnection plays 
a role in solar activity as a mechanism for converting magnetic energy into heat 



and kinetic energy (Priest & Forbes, 2000). The exact nature of this role, how- 



ever, remains elusive. The concept of magnetic X-points being the site of particle 



acceleration and heating, for example, was proposed over 50 years ago by Gio 



vanelli 


( 


1946 


) and 


Hoyle 


( 


1949 



the Masuda flare (Masuda et a/., 1994) that this could be confirmed observa- 



tionally. Currently, the standard thick target model for solar flares by Brown 
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( 1971 ) hypothesises magnetic reconnection at the top of the flare loop as the par- 
ticle acceleration mechanism for precipitating electrons. In the CSHKP model, 
magnetic reconnection is believed to be the method by which a rising plasmoid 
is untethered from the solar photosphere. Observational evidence for these pro- 
cesses so far consist of secondary effects such as particle acceleration, heating 
effects and rising loops. Confirmation of the actual reconnection process has not 
yet been seen. The extremely short length scales over which magnetic reconnec- 
tion occurs ('~meters; 



Petschek, 1964) is well below the resolution of even the 
most technologically advanced instruments. Resolving this process is essential 
for fully understanding the processes involved in energy release on the Sun. 



7.3.2 What is the nature of CME initiation? 

CMEs have been the focus of the solar community since Col. Sabine hypothesised 
a connection between the Sun and the Earth more than 150 years ago. Since 
then, significant progress has been made in understanding the nature of these 
eruptions. Their kinematics have been analysed from many viewpoints over their 



Asai et al. 


2006 


Chifor et al, 2006 


2007 



Feynman & Martin 1995; Gallagher et al., 2002a) to the heliosphere (e.g. Byrne 



et al. 2009 Gopalswamy et al, 2009a Maloney et al. 2009; Schwenn 2000). 



Despite this however, much remains unknown regarding CMEs. While many 



models exist to explain their initiation (e.g. Antiochos, 1998 Forbes & Isenberg 



1991 Moore & Sterling 2006), we are far from isolating one in particular that is 



accurate every time. Observations of X-ray sources at the base of these features. 



such as those observed by e.g. Raftery et al. (2010b) and Krucker et al. (2007 
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2008a||b| [2009| may help to shed some hght on the nature of their initiation 
and acceleration. Unfortunately, the dynamic range of RHESSI, the instrument 
used to image these X-ray sources, is such that any chromospheric emissions will 
dominate the spectrum, making it impossible to observe coronal emission. Thus, 
occulted flares are the most reliable. This however, severely limits the number of 
events suitable for study. 



7.3.3 The connection between EIT "waves" and CMEs? 

Since their discovery in 1998, so-called EIT "waves" have been a topic of hot 
debate. Initially, they were proclaimed as fast-mode magnetoacoustic waves: the 



coronal counterpart of chromospheric "Morton waves", following Uchida (1968). 



EIT waves were believed to be fast-mode magnetoacoustic waves as this would 



explain the EUV brightenings, slow expansion rates (Wang, 2000) and reflection 



and refraction at coronal hole boundaries ( Gopalswamy et al. , 2009b ; Long et al. 



2008 Veronig et al.. 2006a). Since then, two more theories have been suggested. 



Attrill et al. (2007); Chen et al. (2002); Delannee et al. (2007) have suggested 



various models by which EIT waves are not true waves, but are the propagation 
of perturbation sources due to a CME lift-off. Finally, a recent suggestion by 



Wills-Davey et al. (2007) states that "EIT waves" are MHD solitons. 



Recently, an event was studied by Patsourakos et al. (2009 ) using both STEREO 



spacecraft. The spacecraft were located 90° from each other at the time of the 
event, which occurred on the west limb from STEREO B's point of view and at 



disk centre from STEREO A. 3D modelling of the event led Patsourakos et al. 



(2009) to conclude that initially, the "wave" is associated with opening field, but 



194 



7.3 Outstanding questions 



eventually disassociates from the magnetic field and travels as an independent 
wave. The analysis however, remains controversial and inconclusive. There is 
still significant work to be done before a definitive understanding of these events 
can be established. Continued observations and analysis of quadrature events 
using the capabilities of RHESSI X-ray imaging is currently the most productive 
method for establishing this. 



7.3.4 How do we predict solar storms? 

The prediction of solar storms, both fiares and CMEs has not yet been tied down. 
The statistics for fiare prediction are such that if one predicts no fiares will be 
produced in a 24 hour period, more often than not, this will be correct. There- 
fore, with such poor statistics, it is difficult to make a quantitative prediction 
for these events. Currently work is being carried out using machine learning 
techniques that consider the complexity, fractal dimension and magnetic field 
strength, among other parameters, to identify regions on the cusp of energy re- 



lease (Colak & Qahwaji, 2009 Conlon et al. 2008; Gallagher et al. 2002b Leka 



& Barnes, 2007; Qahwaji & Colak 2007 Wheatland, 2005). The use of fiares 



as a mechanism for CME prediction is also important. Svestka (2001) pointed 



out that fiares are excellent indicators of CME activity. Using fiare monitoring 
resources such as www.solarmonitor.org to gather observations into one central 
database can help in making the best informed and complete predictions. The 
process of modelling these events from the Sun out to the Earth is as complex as 
the initial prediction. Where a CME occurs is an important aspect of its poten- 
tial impact on Earth. A CME on the west limb of the Sun will not pass close to 
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the Earth. However, a CME that erupts from the east hmb has a good chance 
of interacting with the Earth's magnetic field. This, of course is important for 
the maintenance of satellites and electricity grids, along with the protection of 



astronauts (Smith & Scalo, 2007). While incremental progress is being achieved 
in this field, there remains a long way to go before it is sufficiently reliable. 
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